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UNDERGROUND SHOCK SIMULATION SESSION 


A REVIEW OF PRESENT ESTIMATES OF NEAR SURFACE 
GROUND MOTIONS FROM NUCLEAR DETONATIONS 
AT OR NEAR THE EARTH’S SURFACE 


by 
Gilbert W. Ullrich, Capt., USAF 
VELA Seismological Center 


A Review of Present Estimates of Near 
Surface Ground Motions From Nuclear Detonations 
At or Near the Earth's Surface 


by Captain Gilbert W. Ullrich 
Abstract 


Estimates of the ground motion environment are important for design 
and cost trade studies of land-based weapons systems that may be sub- 
jected to the effects from both single and multiple nuclear detonations. 
In addition, these estimates provide a basis for evaluation of simula- 
tion techniques that are designed to replicate this environment. Assess- 
ment of the single burst environment estimates indicates that the airslap 
induced component has been previously characterized. Characterization of 
the low frequency component, which has been less successful, should be 
divided into three regions based on range from the detonation point. 
Characterization of low frequency motions in the close-in range (less 
than about 900m from a megaton burst) is presently uncertain because 
of possible significant differences between high explosive and nuclear 
effects. Characterization of low frequency motions in the far-field 
region (greater than 1500m from a megaton burst) is improving with 
understanding developing that the governing physics is Rayleigh Wave 
effects. Characterization of the transition region (between approxi- 
mately 900m and 1500m from a megaton burst) is incomplete but the magni- 
tude of motions may be estimated. For sufficiently separated bursts in 
~ a multiple burst environment, superposition provides only a baseline 
with at least enhanced airblast effects and spall contributing to sig- 
nificant variations. 


INTRODUCTION 


Estimates of the ground motion environment are important for design ~ 
and cost trade studies of land-based weapons systems that may be subjected 
to the effects from both single and multiple nuclear detonations. In 
addition, these estimates provide a basis for the evaluation of simula- 
tion techniques that are designed to replicate this environment. There- 
fore, as a background for the discussion of ground motion simulation, we 
shall review the present understanding of near-surface ground motions from 
nuclear detonations on which these estimates are based. In this review, 
we highlight areas of increased understanding that have developed during 
research in support of the development of Multiple Launch Point (Ref 1) 
missile basing systems. We shall first discuss the ground motions 
expected from a single nuclear detonation and then the effects of 
multiple, near-simultaneous detonations. We limit our discussion to 
expectations of near-surface motions resulting within a few seconds 
from megaton-size detonations within a few kilometers range. The 


discussion should, to first order, also be applicable to kiloton size 
detonations with time and dimensions scaled by cube-root of the yield. 
The limitations on time and range excludes some possibly important, 
late-time effects that are poorly understood. 


TYPES OF SINGLE BURST GROUND MOTION 


Ground motion is produced from the detonation of a nuclear device 
near the earth's surface through both the initial coupling of radiation 
from the device (Ref 2) and by the load resulting from the expanding air 
shock which results from the energy deposited in the air (Fig 1). If 
the height of burst is sufficient, no radiation will be coupled to the 
ground and all ground motion will result from the air shock. Thus, 
airblast-induced ground motion is applicable to a near or above-surface 
detonation, while radiation-coupled ground motion is applicable to near- 
surface and buried configurations. We define a buried detonation as one 
in which no significant air shock is produced. 


The expanding air shock can be viewed (Ref 3) as a series of point 
loads being applied to the surface of the ground as a function of time. 
At a point just below the ground surface and away from the region below 
the detonation point, three categories of signals from the series of 
point loads occurs. One category of signals results from the airblast 
load directly above the point being considered. The ground motion 
resulting from this category is designated airslap induced ground motion. 
A second category of signals results from the transmission through the 
ground of airblast-induced loads applied closer to, but not at, ground 
zero. The ground motions resulting from this category of signals will] 
be designated upstream-induced ground motions. The third category of 
signals, designated as direct-induced ground motions, results at least 
partially from the airblast loads applied in the ground zero region. 
The reason for separating this last category of signals from the 
upstream-induced motion is that the directly coupled energy may also 
contribute, or even dominate, the direct-induced ground motion. 


As a result of the deposition of energy into the ground by both 
direct coupling and the airblast loading, stress and momentum fields 
are established which are then either transmitted away by slower 
propagating waves or dissipated within the region of interest. Trans- 
mission mechanisms include shear waves and surface waves; dissipation 
mechnisms include plastic deformation and work against gravity. Thus, 
while initial motions are generated as the result of stress wave propa- 
gation, additional motion may be the result of either slower propagating 
waves or insufficient time to dissipate the momentum. This difference 
is important because cube-root-of-the-yield scaling does not apply to 
all the mechanisms that stop motion. 


"[RALUUe YSP{qULe JY} YFLM JUBAUNDUOD SUOLJOW de,suLe 
BU} YZLM SUOLJOW puNouH SNOLUPA JO SaWL} LBALWUe ay} SJoLdap weubeLp wamo, ayy, “UoLBba 
(4a}e49) Ouaz punoub ayz wWouy ACME HULAOW ZseLqule 343 s}o.dap weuberp doz ayy “(ph J$2y) 
31LS snoauabowoy e uo euawouayd yOYUS punoub a4} Jo UOLZeJUaSaUd|au OLJeWAaYyIS “| suNnbl4 


WY ‘JONVH 
| 4 
| ye < 
sh > 
| P . 
NP 4 
| aS S = 
| Ss S 4 = 
S <> . = 
| S Ss °” 
SS s > 2 
ONINNAULNO | eS Co = 
| a = 
Wz ™ 
I | = 
VA | 2 
NOILISNVHL JIWSITSHAd NS = 
” 


dJINGNI 
197410 


MIOHS 
annoys 


Because of the relative propagation speeds of the air shock and 
signals propagated through the ground, two major regions of near-sur- 
face ground motion occur. In the first, or close-in region, the air 
shock propagates faster than signals transmitted through the ground, 
and therefore, the airslap-induced ground motions occur first. In the 
second, or far-field region, the air blast propagates slower than 
Signals through the ground, with the result that the airslap-induced 
motions are not the first to occur. This latter condition is termed 
outrunning. In addition, substantial evidence exists that these two 
regions are separated by a third, distinct region of ground motion 
where transitional phenomena, between close-in and far field phenomena, 
occur. 


In all three regions the ground motion may be separated into two 
distinct components based on frequency (Fig 2). The high frequency 
component (greater than approximately 5 hertz) is dominated by the 
airslap induced ground motion, with the timing of this component 
associated with the arrival of the air shock. The vertical ground 
motion associated with this component depends on the air shock over- 
pressure, the soil compressibility, and the air shock impulse up to the 
arrival of the signal reflected from the first material interface with 
significant impedance increase. The physics of this phenomena has been 
understood with several successful prediction techniques available 
(Refs 3, 4, amd 5, for example) depending on the desired level of 
sophistication. While the high frequency component of horizontal 
ground motion has been somewhat less adequately characterized, important 
horizontal effects were usually dominated by the low frequency components. 


The low frequency components (approximately .5 to 5 hertz) show 
distinctly different characteristics depending on the ground motion 
region. In the close-in region, which generally extends to 600 to 900m 
(~ 2,000 to 3,000 ft) from ground zero for a one megaton event, the low 
frequency ground motion ie thought to be dominated by a single phase of 
upward and outward motion that was a combination of upstream induced and 
direct-induced signals followed by “crater related" motion that 
attenuated repaidly with increasing range (Refs 3, 4, 5). In the far- 
field region, beyond about 1,500m (~ 4,800 ft) from a one megaton event, 
the low frequency component is a nearly sinusoidal motion, termed "ground 
roll" (Ref 4), that contains at least a few cycles of motion. The 
maximum amplitude of the ground motion attenuates slower with range than 
in the close-in region. In addition, since the airblast is usually 
traveling slower than the seismic velocities of the earth materials, 
first motion is upward for points at or near the ground surface, and 
significant motions may occur prior to, or during, airblast arrival. 

In the intermediate region between approximately 900 and 1,500m (~ 3,000 
to 4,800 ft), there is no attenuation, and possibly increases, of low 
frequency amplitude with range. Waveforms in this range change from the 
close-in type to the far-field type in a presently uncharacterized manner. 
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ESTIMATES OF LOW FREQUENCY MOTIONS FROM A SINGLE BURST 


Two types of techniques have been used to estimate the magnitudes 
of low frequency ground motions. The first type is numerical simula- 
tions, using basic physical relations in continuum mechanics computer 
codes, of megaton yield nuclear events. This technique, however, can 
only be accurate if all the physical processes important to ground 
motions are adequately described in the codes. Since only limited 
ground motion data (particularly in the close-in region) exists from 
previous high yield nuclear detonations on the earth's surface (Ref 6), 
limited direct tests of those physical relations are possible. Instead, 
data measured during high explosive tests provide the prime test of the 
physical relations. Previous attempts have been notably unsuccessful 
in reproducing the close-in ground motion (Ref 7) but have had some 
success in reproducing the far-field motions (Ref 4). 


Thus, empirically based techniques (Refs 3, 4, 5), which place 
heavy reliance upon generalizing the ground motion data measured dur- 
ing high explosive experiments, have been used as an alternative to 
estimating the direct induced/"crater related" component and, to 
a lesser extent, the ground rol] component of ground motion. These 
empirical techniques, however, require extrapolations in (1) yield, 
(2) type of explosive, and (3) geologic conditions to provide esti- 
mates of ground motion conditions. We will examine the validity of 
those extrapolations in light of some recent results. 


CLOSE- IN 


Recent numerical simulations, with material properties based on 
in-situ test data, have been able to replicate the arrival and, to a 
lesser extent, the peak velocity of the low frequency component of 
motion in this range (Ref 8). However, major discrepancies between 
empirical and numerical predictions still -remain in the expected wave- 
forms and, therefore, maximum displacements (Fig 3). Tke present 
empirical procedures are dominated, in this range, by data from high 
explosive experiments. The data from high explosive tests does indeed 
consistently show the types of waveforms indicated by the "crater 
related" phase of the empirical predictions of Fig 3, but the unscaled 
durations are on the order of 0.1 to 0.3 seconds (Refs 4, 7, 9, 10, 11). 


Within the close-in range, hydrodynamic air shock calculations 
indicate essential differences between the air overpressures from high 
explosive and nuclear sources (Fig 4). In the numerical model of the 
high explosive event the positive phase of overpressure is terminated, 
in this close-in range, by the expanding detonation products which 
sweep the atmosphere away for a period of time (a phenomena supported 
by measured data that, however, may be questioned). Thus, the HE air- 
blast within this region not only has a much reduced positive phase 
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Figure 5. Comparisons of vertical motions at the 0.46m (1.5 ft) depth 
and 7.6m (25 ft) range from a single charge (MB I-2) with motion at the 
same depth and 8.0m(26.25 ft) range from the nearest charge interior to 
a multiple burst charge array (MB I-8). This comparison shows that the 
upward displacements are inhibited by later airblasts such as that causing’ 
the second downward velocity spike on I-8. 
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duration, but may also have a very pronounced negative phase during 
which the ground surface is almost free of atmospheric pressure. This 
difference may result in an essential difference in the mode of defor- 
mation from high explosive and nuclear sources. In the high explosive 
experiments tensile failure and pore air expansion (to be described in 
the multiple burst discussion) may dominate; while these effects may be 
much less significant for ground motions from high yield nuclear explo- 
sions. 


This suggestion of HE/NE difference in the close-in region should 
be considered only a hypothesis, to be used with caution. Some evidence 
supporting this hypothesis was obtained during the Phase I experiments 
of MISERS BLUFF (Ref 12). Comparisons between motion measurements at 
similar locations from the nearest charge showed that, where airblasts 
from other charges prematurely terminated the negative phase from the 
closer charge, the long duration of the vertical component of crater 
related motion was significantly reduced (Fig 5). We note, however, 
that similar differences did not occur in the horizontal motions. 
Further, inclusion of a crude pore air model in recent pretest predic- 
tions, that appear to be consistent with measured data, of the MISERS 
BLUFF Phase II single detonation, resulted in predictions of important 
contributions to crater related motions. The present result of this 
discussion is that significant uncertainty in close-in, low-frequency 
ground motions from nuclear detonations exists. 


FAR FIELD 


While recent analysis has only increased questions related to the 
close-in ground motions, analysis of the far-field ground roll is 
developing a physical understanding of this component of ground motion. 
In particular, numerical simulations (Ref 13) of the Pre Mine Throw IV-6 
and Pre Dice Throw I[I-1 TNT detonations have demonstrated a good capabil- 
ity of reproducing the low-frequency ground roll measured during both 
events (see Fig 6 for example). The two tests were conducted on signifi- 
cantly different sites (dry playa for Pre Mine Throw and essentially 
saturated clays and sands for Pre Dice Throw) and showed substantially 
different ground-rol] frequencies (8 hertz for PMT and 2 hertz for PDT) 
and durations. The numerical simulations reproduced these differences 
in spite of using simple site models with elastic properties only, not 
including crater-induced effects, and using a nuclear airblast loading 
(of twice the high explosive yield) that was much abbreviated in the 
ground zero region. 


These results led to a further modeling of ground roll using an 


analytic solution of Rayleigh Surface Waves on layered, elastic sites 
(Ref 14). The analytic results demonstrated the same capabilities to 
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Figure 6. Comparisons of calculated and measured vertical velocities for 
Pre Mine Throw IV-6 at the 1.5m (5 ft) depth and the 97.6m (320 ft) and 
357m (1170 ft) ranges. The top is the nearer range. 
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model the differences in ground roll between the two test events, again 
provided that a much abbreviated airblast model was used in the ground 
zero region. 


While some questions remain about the general applicability of these 
analytic and numerical results, the trend of the analysis indicates that 
the dominant component of ground roll in the far-field ground motions is 
a long wavelength Rayleigh Wave, with both amplitude and frequency that 
depend on the elastic properties of the site. We also suggest that, 
although there are the differences previously described in the close-in 
airblasts, these Rayleigh Waves are generated from portions of the air- 
blast loadings that may be sufficiently similar between high explosives 
and nuclear devices of twice the yield. We speculate that the loading 
characteristics that control are the gross energy, in the frequency ban 
of the natural Rayleigh Wave frequency for each site, coupled to a depth 
of 1/2 to 1/3 the depth comparable to the range of interest. If these 
hypotheses are correct, then (1) present ground rol] relations should be 
modified to express the dependency on elastic properties of the site, 
(2) ground roll in surface nuclear events should be similar to ground 
roll in surface high explosive events of one-half the yield, and (3) 
ground roll should be simulated by replicating only the controlling 
characteristics of the upstream Toad. 


TRANSITION 


Finally, both recent numerical simulations of nuclear events (Ref 15) 
and empirical predictions based on measured data (Ref 5), have indicated 
that a transition region should exist between approximately 900 and 1,500m 
(~ 3,000 to 4,800 ft), from a one megaton nuclear detonation. Numerical 
results that had a high density of output in this region showed that peak 
motion values were very sensitive to, and sometimes increased with, range. 
To first order, however, these results were consistent with no change in 
magnitude of motion with range. Because of complexity and limited atten- 
tion, waveform characteristics in this region have been inadequately 
described. 


GROUND MOTIONS FROM MULTIPLE DETONATIONS 


The ground motions from a single nuclear detonation provide a basic 
background for all land based systems designed to survive a nuclear 
attack, However, the emergence of attack vehicles with multiple war- 
heads has resulted in the need to describe the environment from multiple, 
near simultaneous detonations. As a result of this need, a two phase 
program of high explosive experiments, designated MISER'S BLUFF, has been 
conducted to measure the effects of multiple, simultaneous detonations on 
airblast and ground motion environments (Ref 16). 
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The principal objective of the Phase I tests was to obtain data from 
a series small-scale multiple-burst experiments to develop waveform 
synthesis procedures for predicting ground motions interior and exterior 
to a multiple charge array (Ref 17). The multiple-burst array design was 
developed by considering an attack on a hexagonal closely-spaced array of 
launch points where two-thirds of the launch points were targeted (Fig 7). 
The unit geometry of such an attack becomes a single ring of 6 charges at 
the corners of a hexagonal, equally spaced from the center point, with 
the distance between charges equal to the distance from any charge to the 
center. The first logical extension of this geometry, which requires 24 
charges, is 3 rings of charges, formed by locating charges at the corners 
of an inner hexagon surrounded by 6 outer hexagons. In this discussion, 
we shall use data froma single charge event (I-2), a 6 charge event (I-4), 
and a 24 charge event (I-8). 


Since an obvious method of predicting ground motion produced by the 
detonation of multiple charges is the direct addition (superposition) of 
the ground motions produced by each charge, the principal questions that 
were asked during the analysis of Phase I ground motion data were (1) 
"In what regions of the multiple burst test bed did superposition work?" 
and (2) “What were the mechanisms that influenced ground motions during 
the times that superposition failed?" Prior to the Phase I experiments, 
a review (Ref 17) had indicated that superposition did not always work; 
however, this conclusion was based on data from multiple charge experi- 
ments where there was thought to be strong inelastic interactions between 
bursts. The charge spacing during the Phase I tests was probably suffi- 
cient to preclude these inelastic interactions. Probably as a result, 
all motions exterior to the charge arrays and at depth interior to the 
charge arrays were adequately characterized by superposition (Ref 12). 


However, superposition was unsuccessful in modeling the motions 
from near-surface masurements interior to the array. Comparison of the 
vertical velocity measured at the .46m (1.5 ft) depth near the center of 
the I-4 array with the band of waveforms constructed by superposition of 
the single burst data (Fig 8) shows at least four examples of the failure 
of superposition. The first failure is the peak downward velocity at 
approximately 25ms (milliseconds) which is underestimated by superposi- 
tion because of the enhancements, beyond superposition, in air over- 
pressure above the measurement location. The second failure is in the 
first upward peak at 50ms which is somewhat overestimated by superposi- 
tion. The third failure is that the downward acceleration following the 
first upward peak is much greater in the superposed wave than in the 
multiple burst data. The multiple burst data indicates a gravitational 
acceleration downward; the single burst data, even before superposition, 
had a downward acceleration that was greater than could be produced by 
gravity alone. Apparently, the material in the multiple burst spalled 
after the rebound; however, the material in the single burst did not 
spall at this depth. Indeed, tensile failure, if it occurs, will result 
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Figure 7. The geometry produced by considering a hexagonal packing of 
launch points where every third launch point is not targeted. The launch 
points are represented by squares, and the targeted launch points are 
represented by squares in circles. 
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in a significant deviation from superposition. The fourth, and most 
impressive, failure of superposition occurs between 60 and 250ms, where 
a major upward and then downward component of velocity occurred in the 
multiple burst experiment, while no similar component was produced by 
superposition of single burst data. 


This fourth failure of superposition requires special consideration 
because the mechanism associated with this failure will become less 
important with increases of yield. With the recognition of the 
occurrence of this phenomena, additional measurements were made on 
event I-8 to further isolate the characteristics of this motion phase. 
In particular, comparison of data from measurements at .30 and .46m 
(1.0 and 1.5 ft) depths near the center of the I-8 test array (Fig 9) 
show that (1) the initiation of the upward motion propagated downward 
at a slow propagation velocity and (2) that the upward motion was 
effectively terminated (earlier than in the 1-4 observation) by an 
even slower downward propagating signal that was generated as a result 
of the airblast arrivals (with peak overpressure an order of magnitude 
less than the initial overpressure) from the second and third rings 
of charges in the I-8 array. Additional observations (Ref 12) showed 
that similar upward motions occurred, at .46m (1.5 ft) depth, elsewhere 
in the Phase I test beds and were observed 20 to 25 msec after the 
beginning of an airblast negative phase. Thus, the fourth failure 
of superposition appeared to be associated with a slowly propagating 
signal resulting from surface underpressures. 


These characteristics have been explained (Ref 12) by considering 
the expansion stress-strain relationship of the mixture of the soil 
and air initially entrained in the soil, using a simple piston analogy 
(Fig 10). By dividing the mass of soil into a rigid piston of mass 
Polo (where op, is the in-situ soil density, and % is a unit length) 
with length (1-n)2. (where n is the volume fraction of air-filled 
voids) and a volume of air with negligible mass and an initial pressure 
of P, for the volume n&.; for a constant unit cross section and assum- 
ing no air flow, we can derive (Ref 12) the stress-strain relation 


p=P, { —L~-}7’ (1) 


where the strain, e«, is defined as 


= . : -] (2) 


and y is the gas constant for adiabatic expansion. In addition, the 
sonic velocity is then 


cz = (He) (ee)? (Tt! (3) 
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Figure 9. Vertical velocities at the 0.3 and 0.46m (1.0 and 1.5 ft) depth near 
the center of MB I-8. The data indicate that the upward acceleration was 


initiated at the shallower gauge 0.01 sec before the deeper gauge and was 
terminated by a downward propagating wave. 
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Figure 10. A piston analogy, ignoring air flow, of the expansion 
of air entrained in soil. The top piston represents the reference 
soil-air system, and the bottom piston represents the expanded 
system, 
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providing an unstrained wave speed of 


When we use (Ref 12) a y for air of 1.4, an initial pressure of 100 kPa 
(~ 14.7 psi), an initial density of 1.7 Mg/m? and an n of .25, then from 
(4) we obtain the unstrained wave speed of 18.3 m/s (60 ft/s) which is 
similar to the wave speeds shown in Fig 9, 


Thus the mechanism, designated pore-air expansion, associated with 
the upward motion that produced the fourth failure of superposition is 
identified as the expansion, caused by a negative phase of overpressure, 
of the air originally entrained in the soil. This mechanism is time- 
dependent and is delicately balanced between overburden weight and magni- 
tude of underpressure. Calculations using models of this process show 
that it becomes much less significant in megaton size detonations, than 
in yields associated with high explosives, because of the increased 
influence of gravitational effects. However, because the process is a 
response to transient reductions in atmospheric pressure (negative phases 
of overpressure) and, as discussed earlier, the region where "crater 
related" motions were measured may be subjected to a pronounced negative 
phase of overpressure, this process may play a major role in producing 
motions that have been attributed to cratering effects. 


SUMMARY 


We thus provide the following assessment of the present understanding 
of ground motions from single, and multiple, nuclear detonations. In the 
single burst case the ground motion may be divided into a high frequency 
and low frequency component. The high frequency component is generated 
by the local airslap. The low frequency component, however, represents 
the transmission of, and response to, stresses transmitted through the 
ground from both local and non-local energy sources. 


The physics of the vertical component of high frequency ground motion 
has been understood, with the horizontal motion Jess adequately treated. 
Thus, prediction techniques exist and simulators (such as HEST and DABS) 
that model the local airblast also model the high frequency ground shock. 


The physics that control the low frequency ground motions has not 
been well understood. The phenomena, and understanding (which is still 
uncertain) associated with low frequency motions differ as a function of 
range, as indicated in Table I for a one-megaton nuclear detonation. In 
the close-in region there are significant differences between the airblast 
characteristics from nuclear and surface high-explosive detonations. 
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These differences may result in an essential unsimilarity in the dura- 
tion of motion, with the surface high-explosive source enhancing this 
duration. Thus, while the maximum velocity associated with the low 
frequency ground motion can be approximated, the displacement and wave- 
forms are in question. Therefore, simulators (such as MINE THROW and 
centrifuge tests) should be evaluated for their phenomena definition, 
and not phenomena replication, potential. In the intermediate region 
the phenomena is ill-defined, possibly as a result of inadequate study, 
but the magnitudes of motion can be estimated. Therefore, one may be 
‘able to define requirements for motion magnitudes; and stimulation, 
(possibly by multiple DIHEST lines) as opposed to simulation, may be 
sufficient. In the far-field region, developing understanding is 
indicating the governing physics should be similar in high explosive and 
nuclear detonations. In addition, apparently only limited explosion- 
characteristics need to be replicated permitting the use of BLEST, 
DIHEST, or buried charge simulators, with simulation verification and 
effects replication possible for this region. 


The single burst environment definition provides a basic building 
block for multiple burst environment definition. For sufficiently 
separated bursts in a multiple burst environment, superposition provides 
a baseline with at -least enhanced airblast effects, tensile failure, 
and pore air effects (described in the multiple burst discussion) 
contributing to significant deviations. Of these latter effects, the 
airblast enhancements and tensile failures prior to airblast arrival 
are expected to remain important at large yields. 
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INTRODUCTION 


Since the advent of the Nuclear Test Ban Treaty, the Defense 
Nuclear Agency and the Armed Service Laboratories have relied to 
a large extent on chemical explosive sources for the proof testing 
of structures and military equipment against the nuclear airblast 
and ground shock environment. The phenomenological investigation 
of nuclear effects, especially investigations concerning crater- 
ing and ground motion, have relied upon chemical explosive to 
simulate the nuclear source. In both usages, simulation implies 
replication of only the nuclear effects of interest. In this 
study we have limited the simulation objective to the direct- and 
cratering-induced ground shock, the airblast-~induced ground shock 
and the airblast. This is graphically illustrated in Figure l. 
How well high explosive simulation sources replicate the airblast 
from nuclear weapons can be readily evaluated since considerable 
data exists for atmospheric nuclear tests made prior to the 
Nuclear Test Ban Treaty. In addition, the computational aspects 
of both high explosive and nuclear effects have been studied in 
considerable detail. 


In May, 1970, Physics International Company (PI) proposed a 
method of using chemical explosives to reproduce the cratering- 
and direct-induced ground motions of a nuclear surface burst. 
This method has become known as the MINE THROW technique. 


The first test of this technique was the MINE THROW I event. 
This event had as its objective the reproduction of the directly 
coupled ground motion produced by the JOHNIE BOY nuclear event 
which was a 500 ton nuclear explosion buried at 58.5 cm in 
Area 18 alluvium at the Nevada Test Site. No attempt was made to 
match the airblast overpressures on MINE THROW I. As the coupled 
airblast energy did not appear significant, it was postulated 
that the crater formation for this event was dominated by the 
directly-induced motion and would not be severely influenced by 
the difference in airblast between the JOHNIE BOY event and 
MINE THROW I. 


The specific technique for designing the MINE THROW I experi- 
ment is as follows. The contour of constant peak stress was 
obtained from the finite difference calculations of the JOHNIE BOY 
event, corresponding to the detonation pressure of the explosive 
used. At each point along that contour, both the pressure as a 
function of time, P(t), and the time integral of P(t), or 
specific impulse, were also determined from these calculations. 

An explosive charge was then shaped in such a way that it repro- 
duced the pressure history and the total specific impulse along 
this contour. In practice, an iterative series of finite diff- 
erence calculations were performed, tailoring the high explosive 
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Figure 1 High explosive simulation of nuclear detonation. 
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to produce the same boundary and initial conditions along the 
above described contour. 


The final charge configuration is shown in Figure 2. A 
comparison of the final craters for MINE THROW I and JOHNIE BOY 
showed an agreement in volume and shape within 11 percent. 
Although JOHNIE BOY had only a few ground motion gages, and 
these were placed at a large radius from the source, the corres- 
ponding gages on MINE THROW I showed very similar displacements 
and ground motion waveforms. A complete description of the MINE 
THROW I event is given in Reference l. 


The results show that the MINE THROW technique is a valid 
technique for reproducing the directly-induced ground motion of 
a nuclear surface burst of a known degree of coupling. 


Subsequent to the MINE THROW I experiment, PI performed 
calculations (Reference 2), on an explosive configuration which 
would simulate the directly-induced ground motion of the CACTUS 
event, a 17 Kt surface nuclear explosion at the Pacific Proving 
Ground. In this case, the relative airblast-induced motions were 
much larger than for JOHNIE BOY. It became clear from these 
calculations that the airblast-induced motions (see Figure 3) on 
the horizontal plane added a significant impulse that should be 
included in the simulation technique. Thus, in addition to the 
target response effects of airblast, and airblast-induced ground 
motions, there are important phenomenological reasons why a 
standard cratering simulation technique should include the proper 
pressure profiles and timing of the airblast. The requirements 
for this airblast addition is shown in Figure 4. 


In March, 1976, Physics International Company proposed a 
method for applying the required close-in airblast loading to the 
ground surface in conjunction with the MINE THROW technique to 
better simulate the cratering-induced and airblast-induced ground 
motions. This technique is schematically shown in Figure 5. 

The MINE THROW charge is shown coupled at its edges to a surface 
charge which extends in a thin sheet above the ground surface to 

a range of approximately 20 meters. The subsurface charge design 
is accomplished by using the technology developed for MINE THROW I. 


The surface charge design required development of a new 
technology. This effort has been underway for about two years. 
During that time a series of one-and two-dimensional calculations 
were performed to establish the design elements such as the 
explosive thickness, the standoff distance and initiation pattern 
for an ANFO surface explosive charge. As a result of that work, 
a preliminary surface charge design was developed. 
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CAST PENTOLITE 


CONSTRUC- 
TED FROM 


6 METER 50 Ib. BAGS 


+ 7.3 METER 


STATUS: FULL SCALE FIELD SIMULATION OF 0.5 kt JOHNIE 
BOY NUCLEAR EVENT SUCCESSFULLY CONDUCTED 
ON DECEMBER 15, 1971 


RESULTS: JOHNIE BOY DIRECT AND CRATERING-INDUCED 
GROUND MOTION AND FINAL CRATER WELL SIM- 
ULATED. NO ATTEMPT TO SIMULATE AIRBLAST- 
INDUCED GROUND MOTION 
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Figure 2 MINE THROW 1-- Direct and cratering-induced 
ground motion, cratering, and ejecta. 
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TIME: t = ty 
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GROUND SHOCK 


STATUS: CONCEPT DEVELOPMENT AND VALIDATION 
PROGRAM IN PROGRESS. INITIAL CHARGE 
GEOMETRIES DESIGNED AND INITIAL COMPUTER 
CALCULATIONS IN PROGRESS. 
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Figure 3 Airblast-induced ground motion. 
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@ CORRECT TIME OF ARRIVAL OF “AIRBLAST” 
@ CORRECT PEAK PRESSURE AS FUNCTION OF RANGE 


@ CORRECT IMPULSE AS FUNCTION OF RANGE 


DESIRABLE 


REASONABLE APPROXIMATION TO NUCLEAR SURFACE 
BURST AIRBLAST FOR P < 100 psi {i.e., AT A RANGE > 
SOURCE DIMENSIONS} 


Figure 4 Requirements for airblast. 
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Figure 5 Charge design concept. 
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© DIRECT AND CRATERING INDUCED. i«*™” 


The procedure for performing the surface charge design is 
shown schematically in Figure 6. The nuclear airblast conditions 
of peak pressure, P, total impulse, I, and time of arrival, TOA, 
along the ground surface are well known. In order to simulate 
these conditions data from both experiments and calculations must 
be obtained for thin slabs of high explosive of total thickness, 
A, elevated a distance, V, above the ground surface. The required 
high explosive results are shown in Figure 6. In general, they 
relate the peak pressure at the ground surface, P', the total 
impulse at that point. I', and the corresponding time of arrival, 
TOA', to both V and A. Once these high explosive results are 
known, it is a straight forward procedure to develop the basic 
charge design to accomplish the simulation of the baseline nuclear 
conditions. The procedure is shown at the bottom of Figure 6. 
Sufficient information is available to obtain both V and A as a 
function of range and also the lighting times, LT, of the explo- 
sive charge vs range. 


An important assumption made in this preliminary design 
process, was that ANFO in thin sheets acted nearly as an ideal 
explosive. There was, at that time, some data in existence which 
indicated the detonation velocities for thin slabs of ANFO might 
be considerably reduced from the quoted Chapman-Jouget (C-J) 
value as shown in Figure 7. A number of 1D calculations were 
performed using the simple geometry shown in Figure 8. First 
order corrections for small thicknesses of ANFO were used as 
shown in Figure 9. The assumptions being recognized, a series of 
experiments were begun at this time, investigating the detonation 
characteristics of thin slabs of ANFO. 


Using this model of an ideal explosive, calculations were 
performed to determine the thickness of ANFO required as a function 
of range which is shown in Figure 10, and the times of initiation 
as a function of range for these thin slabs as shown in Figure ll. 
These resulted in a preliminary charge design designated 1A 
which is shown in Figure 12. Recognizing the assumptions used in 
Charge Design 1A, experiments investigating the detonation 
characteristics of thin ANFO slabs were performed to check its val- 
idity. These experiments consisted of detonating thin slabs of 
ANFO (17.5 cm to 7.5 cm) thick, 0.6m wide X 1.0m long and with a 
50 cm void. Time of arrival pins were installed along two lines 
to determine velocity of detonation as shown in Figure 13. These 
results show the presence of a “valley” in the detonation velocity 
vs range curve for distances within 50 cm of the booster. It 
appears that ANFO, which is a composite explosive, does not attain 
its Chapman-Jouget. or steady state detonation characteristics 
immediately upon igniting, but because of the multiplicity of 
chemical reactions taking place, takes some time to reach the 
classical C-J conditions. During the transition period between 
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Figure 6 Standard source surface charge design. 
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Detonation trajectory for undiluted 94/6 ANFO. 
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Figure 8 Computational geometry. 
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FOR ANFO THICKNESSES ~ 100 cm 
Poy = 55kb 
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Figure 9 First order correction for small ANFO thickness. 
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Figure 10 Explosive thickness versus range, charge design lA. 
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Initiation time versus range, charge design 1A. 
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Figure 12 Charge design 1A. 
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Thin ANFO experiments. 
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boosting and the C-J conditions, the detonation velocity, D, is 
substantially lower than the C-J detonation velocity. The 
velocity slowly increases so that at the end of the valley, 

D_. for bulk ANFO, is reached. Very little previous work has been 
pérformed to examine this phenomena of a non-ideal explosive, yet, 
this data is crucial to this technique since thin sheets of ANFO 
must be used to simulate the nuclear airblast-induced ground 

- motion. : 


Lately, a grain burn model has been implemented in the PISCES 
computer code to attempt to account for the non-ideal behavior 
of ANFO. Some initial 1D and 2D calculations have been performed. 
These calculations have been based on an ANFO grain burn model, 
shown in Figure 14. The results of these calculations are shown 
in Figure 15. Much more needs to be done to verify the adequacy 
of the model. It is felt that the parameters needed for the 
ANFO grain burn model, must depend only upon the properties of 
the ANFO itself, thus, 1D spherical and cylindrical calculations, 
uSing the grain burn model, should show the presence of the detona- 
tion valley if the parameters are to be picked correctly. 


A full-scale test, based on charge design 1A, could be con- 
structed by using four foot by eight foot trays as shown in 
Figure 16. These trays would be made of low mass paper honeycomb 
Supported by low mass styrofoam pillars. The depth of the trays 
would be adjusted for the correct depth of explosive. The styro- 
foam pillars would be cut to the proper height above the soil 
surface. Several hundred of these trays would be required in an 
arrangement as shown in Figure 17. This shows an artist's con- 
cept layout of how the below ground MINE THROW charge is inte- 
grated with the above-ground airblast charge. A suitable protec- 
tive housing must be provided for this entire arrangement, since 
it would be very susceptible to weather conditions, especially 
wind loadings and precipitation. Each of the trays could be 
detonated at the proper time in concentric arrays as shown in 
Figure 18. Designs have been completed and tested which show 
the feasibility of these initiation geometries using both 
electrical timing and pyrotechnic delays. 
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NON-IDEAL EXPLOSIVES: 

e CHEMICAL ENERGY IS RELEASED OVER A FINITE 
PERIOD OF TIME AFTER THE ARRIVAL OF THE 
DETONATION WAVE. 


e AN ACTIVATION ENERGY IS REQUIRED TO INITIATE 
THE ENERGY RELEASE. 


THUS GRAIN BURN MODEL ASSUMES: 


dF 


de G(E,, 7), where 


F = FRACTION OF HE CONSUMED AT A POINT 


E, = ACTIVATION ENERGY = aE, 
E, = TOTAL ENERGY RELEASE/UNIT WEIGHT 


7 = TIME OVER WHICH E, IS COMPLETELY RE- 
LEASED AT A POINT 


Figure 14 ANFO grain burn model. 
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DETONATE 1D CALCULATIONS 
| | | | ; | | | (PLANAR SYMMETRY) 
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Figure 15 Thin ANFO calculations, grain burn model. 
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Figure 16 Typical panel. 
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1. Introduction 


Cooper (1975) asks the question "Cannot small scale experiments 
credibly be used to study phenomena (e.g., height-of-burst and layering effects) 
and scaling rules?" This question is answered and implications on the modeling 
of ground shock as well as the response of buried structures is addressed. 

The approach is based upon performing subscale experiments which are 
dynamically similar to their full scale prototypes. Two events of different 
sizes are called similar if each of the pertinent thermomechanical field 
variables (e.g., stress, velocity, etc.) are equal to within a simple scale 
factor for the two events. To use this concept, it is necessary to correctly 
identify the pertinent fields and to control all the test parameters to ensure 
similarity. The more common practice of extrapolating the results of small 
scale experiments without due regard to similarity requirements will also be 
briefly discussed to highlight pitfalls of such a procedure. | 

Section 2 is a complete similarity analysis for the thermomechanical 
response of a continuous media. Here the theory is developed and similarity. 
requirements common to all cratering and ground motion applications are derived. 

Section 3 gives results demonstrating the applicability to the problem 
of cratering. In particular, a soil strength model is presented which explains 
the dependence of apparent crater volume upon charge size and material strength, 
where the latter is given by a Mohr-Coulomb representation. 

Section 4 gives the results of using this same soil strength model to 
develop a set of non-dimensional parametric curves for the prediction of ground 
motion. Comparison with limited field data is given and a discussion suggesting 
refinement of the theory and further verification of its applicability is 
included. | 

Section 5 is a discussion of the feasibility of various types of 
subscale experiments which must be conducted on a centrifuge to achieve the 
elevated gravity field strength required by similarity. These include ground 


motion as well as structure-soil interaction. 
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Cs Similarity in Ground Motion Experiments 


Ground shock events are of a general class of phenomena involving 
detonation sources and resulting dynamic stress and deformation fields. It will 
be shown that a similarity analysis derived for cratering due to near surface 
bursts is also applicable for ground motion. Schmidt and Holsapple 
(1978a,1978b) present a general analysis of the requirements of similarity with 
specific application to cratering. A condensation of that analysis appropriate 
to the ground shock problem is repeated here. 

The parameters and fields used to define and characterize ground shock 


phenomena include the five fields: 


x(X,t) the spatial position vector (1.1) 
Q (Xt) the mass density per unit volume Clee) 
A(X, t) the stress tensor field (123) 
B(X,t) the body force per unit mass (1.4) 
e(X,t) the internal energy per unit mass (1.5) 


where X is the initial position and t js the time. These fields determine the 


auxiliary fields 


aX,t) = R(X, t)-X the displacement vector field ; (1.6) 

¥(X,t) = d x(X,t) the velocity field (1.7) 
dt . 

a(X,t) = d ¥(X,t) the acceleration field (1.8) 
dt 

E(X,t = Grad x(X,t) the deformation gradient tensor (1.10) 

LOGE = id 6%, yILE(X,t)T | the velocity gradient tensor Gbaae) 


The operator d/dt denotes the time derivative at fixed initial position 
X, while Grad denotes the gradient with respect to X at fixed time t. A 
Superposed arrow denotes a vector and a subposed tilde denotes a tensor. The 
superscript -1 denotes a tensor inverse, and the superscript T denotes the 
transpose. The tensor Lis the unit tensor. The symetric part of the tensor 
LE=1) is the infinitesimal strain tensor and the symmetric part of Lis the rate 


of strain tensor used in fluid mechanics. 
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The fields are not independent but must satisfy the balance equations 
of mechanics. In the absence of heat conduction and heat supply, these equations 
have the form 
balance of mass 


Qdet (F) = @, (2) 


balance of angular momentum 


" 
eo 
—_ 
Ww 
— 


es 


balance of linear momentum 


and balance of energy 
Qe= Tr (TL). (5) 


nee QR 9 is the initial density, div refers to the divergence with respect to 
x, and Tr denotes the trace of a tensor. In addition, these fields are related 
by whichever consitutitive equations describe the response of the media in 
question and by the initial and boundary conditons. 

At surfaces of discontinuity; e.g., detonation fronts, these equations 
are augmented by jump conditions. For a shock or a detonation wave with local 
normal T moving at speed U into an undeformed and unstressed material, the 


following equations apply 


mass: Ro U = Qu - Q (ven) (6) 


momentum : Th + Rev = 0 (7) 
m 

a> = ss =” 
energy: Ro Ule-e, ue ve vV-Q) +n s Tv = 0 (8) 
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and a compatibility condition 


U(F-l)i + V = 0 (9) 
where the variablesQ, V, J, e, E refer to the properties behind the front, Q is 
*he energy per unit mass added at the jump for a detonation wave, eo andy acne 
the energy and the density in the undeformed material ahead of the wave. 

Two different solutions corresponding to two different experiments, 
each satisfying the entire set of eqs. 1 through 9 are to be compared. All of 
the quantities associated with the second solution will be denoted by primes. 


These two solutions are said to be similar if the following relationships: 
K(X! t') 2 oC, KK t) (10.1) 


Q(X t') = HQ (Kt) (10.2) 


TO tty = %, 1H t) (10.3) 
Bit!) = &X, B(X,t) (10.4) 
e'(X',t') =%, e(K,t) (10.5) 
Q'(K',t") =X) a(%,t) (10.6) 


Koy X (10.7) 
for homologous time 
i =X, t. (10.8) 


The derived fields are related as a consequence of their definitions by 
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"Ox, th) =o, UK, t) (11.1) 


x 
v'(xtt') =%& W(X, t) (11.2) 
r=. 
t 
a (x't') = % “BGK, t) (11.3) 
m,)? 
F'(k) 0!) Xx F(X, t) (11.4) 
ro4 
X 
Le te Lt) (11.5) 
* 


Here the eight various factors %, are constants called scale factors. Note that 

a =X if in both solutions, the reference position X is the initial position: 
“x(t=0) =X. This is henceforth assumed. leaving seven independent. scale factors. 
Thus, in particular, eq.11.4 states that the deformation gradient must be the 
same at homologous points and therefore the strains are identical also. In more 
detailed terminology, the two solutions are said to be geometrically, 
kinematically and dynamically similar (Langhaar, 1951). 

The question of the existence of two different experiments related by 
these similarity requirements is to be investigated. Both sets of fields must 
satisfy the balance equations given above, regardless of the form of the 
constitutive equations. This requirement will give certain restrictions 
relating the seven scale factors. Additional restrictions from specific 
constitutive equations are derived separately below. 

Assume that a primed solution exists, substitute the primed fields into 
the balance equations 2 through 5 and use the similarity relationships (eqs. 10) 
to get 


Repeat (12.2) 
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cere Sh Teen) (14) 
sar aa 


The unprimed fields must also satisfy the balance equations. Hence equations 
12.1 and 12.2 are satisfied identically. 

For the balance of linear momentum eq. 4 can be used for div T in eq. 
13, giving 


At Ag &Xx loT=[%r xy xX,lOb. (15) 
x eye IN x, . ys 


This must hold for all] X and t. Thus, unless the acceleration field a is itself 
a scalar multiple of the body force b or either a or b is identically zero, it 


is necessary that 


T_ _ XQ x (16) 
Ky, OM) 
and 
XT = OX" (17) 
X x 


The most common body force per unit mass is a constant vector proportional to a 
constant gravity, 9g. In this case, the acceleration field which is a scalar 
multiple of b is at most a rigid body motion. If b= 0, then eq. 17 is not 
required. If a= 0, as in statics, then eq. 16 is not required. 
Now substitute Re from the balance of energy eq. 5 into eq. 14 to get 
X_eXe -%r | Tr (T 


ar ee 


L) = 0 (18) 


for all x and t. Again, assuming the stress and velocity gradient fields are 
neither identically zero, it is necessary that 
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RoK Ay . (19) 


Lastly, requirements brought about by detonation waves in the explosive 
Can now be included. For this case, both solutions must also satisfy the jump 
conditions across the wave (eas. 6-9). An analysis identical to that just given 
for the balance equations produces only one additional restriction which 
involves the scale factor Kp for the energy density of the explosive Q 


x =Ar 
TAS 
Altogether then, the balance equations plus the jump conditions provide 
four similarity restrictions among the seven scale factors relating the variable 
fields as defined by eqs. 10.1 through 10.8. Eqs. 16 and 17 are a consequence 
of the balance of momentum, eq. 19 results from the balance of energy and the 


. (20) 


energy jump condition supplies the remaining constraint, eq. 20. The balance of 
mass and the balance of angular momentum provide no restrictions. 

It is convenient for the present application to consider &,, XK, and Xe 
as independent. The four restrictive conditions (eqs. 16-20) can then be used 
to solve for the remaining scale factors: 


X= pal Xp (2.1) 
a 
AX, =X /%)° (21.2) 
x, = oi | (21.3) 
org 
tl 
xX, E 2 (21.4) 
LN 


Therefore it can be seen that, considering only the balance equations and the 
jump conditions, nontrivial similar solutions with arbitrary scaling of size, 


time and density are allowed. The four remaining scale factors must then 
Satisfy the four equations, 21.1 through 21.6. 
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Analysis of Specific Constitutive Equations 


The complete solution to the deformation and flow of any continuous 
medium depends on, in addition to the balance equations and the jump equations 
given above, the constitutive equations that describe the behavior of that 
medium. Therefore, while the above similarity requirements (eqs. 21.1 through 
21.4) are necessary, the question of their sufficiency must be answered. Note, 
however, that they are general and apply to all materials. 

The nature and type of equations that describe the material behavior 
for a given medium, and their role in determining similarity requirements is now 
considered. There are three media in the present problem: the soil, the 
explosive and the overlying air, at least for near surface bursts. Each of 
these media is extremely complex. The complete characterization of the soil 
requires concepts of compressibility, nonlinearity, yield, fracture, porosity, 
cohesion and others. As disussed by Schmidt and Holsapple (1978a, 1978b), it is 
not necessary to pick definite constitutive equations for the purposes of a 
Similarity analysis. The complexity of the soil behavior leads one to believe 
that any scaling or modeling of the soil itself is hopeless, except perhaps in 
extremely simple problems (Pokrovsky and Fyokorov, 1969}. As a consequence, one 
expects that the same soil medium must be used in two experiments if similarity 
is to be achieved. For the purposes of this discussion it will be assumed that 
the same explosive is used in both, and only the size is varied. For this case 
it is necessary that Xe= 1,%q = 1 and from eq. 21.4, OX, =AX,. Hence there 
remains only one independent scale factor A, and all others are determined, as 
follows 


Xx, =] (22.1) 
AKpe= Why (22.2). 
42 (22,3) 
X= 1. (22.4) 


In this case, similarity is possible only insofar as the constitutive 
equations are consistent with this scaling. The stress, the strain, the 
density, the internal energy and the heat of detonation must be the same at 
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homologous points. The body force must scale as the reciprocal of the size, and 
the scale facors for time and distance are equal. Consequently, if&X, = 1/19 so 
that a 1/10 size scale experiment is to be performed, the body force must be 10 
times larger. All events will occur in 1/10 the time over 1/10 the distance, 
and all velocities will be the same. 

Schmidt and Holsapple (1978a, 1978b) have addressed the question of the 
consistency of material constitutive behavior with the set of eqs. 22.1-22.4. 
They have proved the statement, "Complete and exact dynamic similarity can be 
achieved between two different experiments of arbitrary size scale in the same 
material as long as the constitutive behavior is 1) rate independent_and 2) has 
no inherent size properites." If either of these requirements is not satisfied 


then the material must also be modeled using the concept of equivalent 
materials. 

Summarizing the above development, the design of a subscale experiment 
that directly simulates a large full scale event is indeed possible. All 
dimensions are to be reduced by a common scale factor GQ,» which implies the 
yield of the device is reduced by &,) - The same soil and explosive type is 
used in both experiments. The body force, i.e., the gravity, must be increased 
by 1. The way that this can be achieved in practice is through the use of a 


centrifuge. 
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3. Subscale Cratering Experiments 


While the primary interest of this analysis is ground motion phenomena, 
a closely related problem is that of cratering due to near surface bursts. 
Recent cratering experiments by Schmidt and Holsapple (1978a, 1978b, 1979) have 
provided valuable experience in the conduct of explosive testing in the high-G 
environment of a centrifuge. Furthermore, as will be discussed in the next 
section, the ground motion due to near surface explosives in various experiments 
was correlated with the crater parameters in these same experiments (Cooper and 
Sauer, 1977). Consequently, it is appropriate to summarize recent centrifuge 
cratering experiments and results here. These subscale explosive cratering 
experiments were performed in Ottawa sand, desert alluvium, clay, and saturated 
sand. Two different explosives, PETN and lead azide, in sizes up to four grams, 
were used. These experiments were conducted using the Boeing geotechnic 
centrifuge at gravity fields up to 520 G's. 

As shown in the previous section, if the body force is increased by a 
factor of 520, then for this subscale experiment X = 1/520. The similar full 
scale experiment is 520 times larger in all linear dimensions. The corresponding 
charge size is (520)°9 times 4 grams, which is a charge size of over 0.6 kilotons 
of PETN, or an equivalent energy of anost 1 KT of TNT. Therefore the centrifuge 
experiments have given a direct simulation of relatively large scale field 
events. This direct simulation technique was demonstrated very successfully 
with a subscale replication of the JOHNIE BOY 500-ton shallow buried nuclear 
event conducted at 345 G's (Schmidt, 1978). 

The results of centrifuge experiments have been invaluable in 
understanding both size effects and material strength effects in the various 
media. The base line experiments in Ottawa sand are shown in Fig. 1 as a plot 


of the cratering efficiencies 
4¥, = Vo/W (23) 


versus the gravity-scaled size parameter 


We eg lap 2 (24) 


°° 
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where 


= crater volume 
soil density 
= charge mass 

= gravity 


Dau =f =< 
" 


= charge energy per unit mass 
& = charge density 


The high-G (equivalent large scale) events done on the centrifuge correlate 
extremely well with 1-G small scale experiments (Piekutowski, 1974, 1975) and 
all experiments fall on the same straight Jine with negligible scatter. The 
slope of this line, on log-log paper, is -0.472 which gives a crater volume 
dependence upon yield 


Additionally for Ottawa sand, the crater linear dimensions scale with the energy 
to the 0.281 power, over an energy range of almost 12 decades. The form of %, 
as given by eq. 24 also includes an adequate representation of the energy 
density term, Q. This was demonstrated over the entire range of Q encountered 
for common chemical explosives (Schmidt and Holsapple, 1978a, 1978b) and over an 
order of magnitude in range for impact experiments (Schmidt and Holsapple, 
1978c). 

The results for the other materials are superposed on the Ottawa sand 
curve and shown in Fig. 2. The features of these results can be comprehended by 
considering a given material such as the alluvium, for a fixed explosive. At 
smal] values of 17>, j.e, small yield, the data is essentially horizontal, which 
implies that the cratering efficiency is constant, and the crater volume is 
proportional to the charge mass. This is a regime where the familiar "cube- 
root" scaling holds. 

For larger values of charge mass for gravity), the results become 
asymptotic to lines that are parallel to the sand curve. For alluvium, this 
Tine is almost the same as that for sand. The clay results and saturated sand 
data (Schmidt and Holsapple, 1979, not shown) both become asymptotic to lines 
above the sand line, but have approdmately the same slope. Thus, in each case 
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there appears to be a transition with charge size and, for sufficiently large 
experiments, cube root scaling is replaced by a crater dimension dependence on 
charge mass with an exponent of about 0.28 for all materials. However, the 
exact value for this asymptotic yield exponent is tentative from the number of 
experiments conducted to date. 

A theory that accounts for this behavior and accounts for the 
differences between the various materials has been constructed and can be 
summarized as follows. Typically there are three measures of Stress or pressure 
that can be identified as contributing to this problem. These are 

1) the magnitude of the shock front 

2) the soil strength 

3) the lithostatic pressure. 

A measure of the strength of the outgoing shock is given by the 
magnitude of the Chapman-Jouget pressure in the detonation wave of the 
explosive. This pressure is given in terms of the explosive properties by the 
relation 


Pog = 2( X-1) 90 (26) 
where & is the perfect gas constant of the explosive products, and is 
approximately 2.8 for common explosives. 


The Mohr-Coloumb strength model for soils relates the yield strength to 
the pressure by 


Y =ctP tan @ (27) 


where c is the cohesion of the material and 6 is the angle of internal friction. 
The lithostatic pressure depends on the depth in a linear manner. A comparison 
of this quantity can be made for different size experiments using the 
lithostatic pressure at the base of the spherical explosive. For hal f-buried 
charges this quantity is 


P = Qga = g{ Ww) M3 (28) 
see 


where, a, iS the radius of the spherical explosive. 
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The relative size and importance of these three pressures can be 
measured by their ratios. The ratio of lithostatic to C-J pressure is given by 


1/3 (29) 


which, to within the factor of 2(%-1) (87@) is the gravity-size parameter TY, 
defined earlier: 


Mg (PPeg) 2 URAC Sle (30) 


SAIN 


The ratio of the yield strength to the C-J pressure is given by 


VP es = ct+P tan # = aS [ 1c re (31) 
Pes 2(6-1) eo + WN, tan As 
cal 
The term in square brackets is used to define the dimensionless parameter 
Sh ee ec +17, tan ¢ (32) 


Xa 


which depends on the material strength parameters c and 6, the explosive 
property Q and the previous gravity-size parameter TT,,, For a cohesionless 
material such as dry Ottawa sand, this material strength parameters is only a 
constant factor of tan ¢ different from TT. 

Various measures of cratering, and in particular the resulting crater 
volume can be expected to depend on these stress ratios. Thus 

TY 8h a Ny) (33) 

The explicit dependence on IY, is a measure of the importance of the gravity 
determined lithostatic pressure relative to the shock strengths. The explicit 
dependence on TT, is a measure of the effects of finite material strength 
compared to the shock strength. This strength term incorporates an additional 
implicit gravity-size dependence because of the tan 6 tem. 
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Examination of existing data for granular materials indicates that the 
explicit dependence on TY, is not as important as the implicit dependence via 
the Ty dependence. Furthermore, Vy depends linearly on f¥. Consequently, it 
seems plausible to combine these two dependences, introducing a single parameter 
incorporating, both iV, and Tr, in the form 


Tr, =, + kT (34) 
where k is a universal constant. Written out this gives 


17 » ae. +TY, (tan p+ k). (35) 
QS 
The constant k is expected to be relatively smal] compared to the tan 4 term for 
materials such as sand. For sand, tan ¢ is about 0.7 so that k should be small 
compared to 0.7. Only when tan 6 = O will the exact value of k be significant 
or determinable. 

A plot of TT, versus this combined strength-gravity-size parameter TY, 
gives an interesting result, as shown in Fig. 3. All results now superimpose on 
a single straight line, with equation as shown. If it is assumed that TQ, 
versus IT, is indeed the same for all materials, then the result when plotted as 
47, versus the original gravity-size parameter TY is as shown in Fig. 4 for the 
various combinations of soil media and explosive type. It is seen that these 
curves reproduce and explain the actual data points extremely well. 

From these results, an explanation of the experimental results becomes 
possible. For small events, TY, is small, and the cohesion term c/Qe dominates 
in TY. There is little gravity or size effect and cube root scaling holds. 
Only for a cohesionless material is this plateau not observed. 

The other extreme is for large size events. In this case IT, is large 
and the cohesion term can be ignored. The strength of the medium is dominated 
by the tan 6 term and by lithostatic pressure. In this extreme, and for a 
given material and explosive, the parameter TY, is only a constant factor 
different from Ole: and a straight line asymptote results on log-log paper. The 
level of this asymptote, compared to the base line sand data is determined by 
the magnitude of the tan 6 of the material compared to that for sand. 

Fig. 5 shows a schematic of the nature of the resulting curve for a 
given material compared to the sand curve. The predicted response is a curve 
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with the two straight line asymptotes as shown. The transition between the 
cohesion-dominated regime and the lithostatic-pressure dominated regime 7S shown 
on the figure as the intersection between the two asymptotes. It is given by 
the value 


Wee gfyi C (36) 
PS) 


—Ols QA(k + tan B) 


corresponding to the critical value of charge mass, : 
Ws eal? | (37) 


¢ Rsk ten BH 


For a given material, the value of ths transitional charge size can be 
calculated or read from the plot given in Fig. 6. Representative ranges for the 
Strength parameters of common materials are also shown on this figure. For 
example, desert alluvium at 1 G requires a charge mass of about 100 tons to 
reach the transition. The behavior below this critical size is cohesion 
dominated. Data from tests smaller than this size at 1 G would not extrapolate 
to the correct large yield behavior. 

Prior to the development of this strength theory, it was puzzling to 
compare the results of cratering in clay as compared to sand. It is usually 
thought that sand is strengthless and that clay has definite strength. This 
would lead one to expect that a given explosive would always produce a smaller 
crater in clay than in sand. However, while the results confirmed this at small 
1, values, large v, showed the opposite behavior. 

The reason for this is apparent by identifying the two components of 
the strength parameter model. Sand is essentially cohesionless, but it is not 
without strength. It has an angle of internal friction of about 25° and thus tan 
6 is about 0.7. Clay, on the other hand has a non-zero cohesion but has very 
smal] 6. Consequently, under large confining pressure, sand becomes much 
stronger than clay. For large scale events, a given explosive will produce a 
larger crater in clay than it will in sand. 

Alluvium has almost the same angle of internal friction as does Ottawa 
sand. In addition, it has non-zero cohesion. This cohesion depends strongly on 
moisture content. Therefore, it shows the cube-root plateau for small yields, 
but becomes asymptotic to the sand curve for large yields. This explains why 
small scale field events cannot be extrapolated to large scale events since they 
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do not fit the same scaling law as do the large event data. This fact takes on 
particular significance when extrapolating to megaton-sized nuclear events. 

This soil strength model for crater volume has been applied to various 
field data using a combined statistical/dimensional-analysis technique (Schmidt, 
et. al., 1978d). This on-going pilot study has provided a significant 
improvement in the data correlation for different size events conducted in 
various soil media over a range of burial depths. 
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4, Ground Motion Predictions 


As discussed above, direct experimental simulation of large scale 
explosive events is possible, on a centrifuge at greatly reduced scale. This 
approach has proved to be extremely fruitful for the cratering problem and may 
also prove to be viable for the near-field ground motion problem. Preliminary 
feasibility studies verify this hypothesis. In lieu of a direct experimental 
program, existing field data furnishes the basis for ground shock predictions. 
Cooper and Sauer (1977) have showed that the apparent crater volume can be used 
as a measure of the energy imparted to the soi] in correlating near-field ground 
‘motion measurements. This approach suffers from the fact that the crater volume 
is itself one of the dependent variables of an exposive event and is probably no 
easier to predict from existing field data than the ground motion itself. The 
correlation of ground motion to crater volume is not useful as a predictive 
tool. However, if one incorporates the crater volume prediction discussed in 
the previous section with the approach of Cooper and Sauer, then a prediction of 
the ground motion emanating from a given explosive event can be made. 

The crater predictions discussed previously are valid for hal f-buried 
charges. In this case the apparent crater volume is predicted by the formula 


Ty = We guiaety, O44 (38) 


where, expanding eq. 35 from the previous section 


ne =c + (tan 6 +k) Ig E Ie. (39) 
QR iQ} Ss | | 


and a value for k = 0.1 has been chosen as consistent with existing data. 
Cooper and Sauer (1977) give, for near surface bursts, a correlation 
for peak horizontal displacement di, as a function of range, r, as follows 


= 0.25| vi |3 (40) 


aT: Cer 


supposedly for all materials and explosives. 
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Equations 38 and 39 can be inserted into eq. 40 to obtain the peak 
displacement as a function of soil and explosive properties. However, the data 
used by Cooper and Sauer to obtain the above correlation consisted only of 
events using TNT. Consequently, using values for the energy density Q and the 
mass density of TNT, the dependence on material strength and explosive mass in 
egs units is 

d, = 0.0244 wlt/? [2.39 107! c+ 1.99 198 wl’? (tan 6 + 0.1) 10°92? (41) 
r3 g aS 


It is convenient to rewrite eq. 41 as 


1 -0.629 


4/3 
dy, me = 0.0244 [2.39 107 (42) 
W 


L 


c+ 1.99 1079 wl’? (tan g + 0.197 
Q 


where the term on the left_is a non-dimensional quantity that will henceforth be 
denoted by the symbol Kin: 
Wig ead ra . (43) 
W 
The value of 1s is a measure of the peak displacement d, at a given range and 
for a given charge size. It is a function of the explosive type and size and 


the soil material properties as given by eq. 42 for cgs units. A completely 
non-dimensional form of eq. 4? can be written as follows 


aN = 0.869 £ (0.1 + tan 6)41T, ‘s qe? (44) 
R 


Figure 7 shows the dependence of ee or charge size W for representative values 
of material properties. The corresponding general 1T, dependence is also shown. 
For a given -soAl such as desert alluvium, the character of the curve 
reflects the differing roles of material strength in the different size regimes. 
At small yields, the cohesion of the alluvium dominates. Consequently the crater 
volume scales linearly with the charge energy, and the ground shock parameter 
Tah is independent of charge size. At large yields the confining pressure is 
much larger and the dominant strength term is due to the angle of internal 
friction # in the usual Mohr-Coulomb strength criteria. In the latter regime, 


the relative cratering efficiency, the volume per unit charge energy, decreases 
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with increasing size and, as a consequence, the predicted ground shock parameter 
decreases with increasing explosive yield. 

For a given material, the transition from the cohesion dominated regime 
to the confining pressure dominated regime depends on the relative magnitudes of 
the cohesion and the angle of internal friction. Since the cratering 
efficiency, and hence presumably also the ground shock, differs markedly in 
these two regimes, it is important to note where the transition is for various 
materials. It is clearly inappropriate to measure data in one regime and 
extrapolate to the other. This observation is of special significance to a 
material with a transition at relatively large yields, where field data is 
likely to be in the cohesion dominated range, but results are desired for much 
larger yields. 

Typical values of TT, from several field shots are shown on Fig. 7. 
They correlate well with the curves. The reason is, of course, because the 
results of Cooper and Sauer which were an empirical fit to these very events, 
were used in developing the theory. However, al? shots shown for zero d.o.b., 
with only one exception (Mole 206), were 20 tons of TNT in more-or-less like 
soil media. However, two events with only small height of burst (tangent above) 
in a very different media (granite) are shown on this curve also. If it is 
assumed that the crater volume would have been a factor of 2 larger at zero 
d.o.b., then these two events also correlate well with the curves shown. Thus 
there is some verification of the correctness of the strength effects. 

Ground shock prediction as determined from Fig. 7 was, as the above 
development showed, based upon two major ideas. The first that of Cooper and 
Sauer (1977) relating peak ground displacement to crater volume based upon 
existing field data. These results, incorporated with the crater volume 
prediction results of Schmidt and Holsaple (1979) led to the ground motion 
results. However, at the resulting level of generality, certain limitations 
should be noted. 

The crater volumes predicted by centrifuge results of Schmidt and 
Holsapple depend explicitly on the following parameters of the event 


charge mass W 

charge energy density Q 

soil cohesion c 

soil internal friction angle § 


Soil density te 
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Their results correlated crater volumes over large ranges of these variables, 
for the case of half buried charges. 

On the other hand, the field data available to Cooper and Sauer for the 
zero depth of burial case were limited to the range of 20 to 500 tons. Of the 
eight field events shown in Fig. 7 all were with TNT and six were 20 ton shots, 
and only relatively small variations in material properties were present. 
Although the 500-ton data is not shown in Fig. 7, it is only a factor of 3 
larger in NV, and would be just to the right of the 100-ton points shown for — 
Mineral Rock and Mine Ore. These latter ‘two points which represent a limiting 
soil type had to be adjusted for an h.o.b. effect to be compared with the half 
buried data. Consequently, the hypothesis of Cooper and Sauer cannot be 
considered to be verified in large ranges in explosive type, charge mass and 
media properties. However, some obvious experiments can now he suggested to 
further verify this type of relationship. These will be discussed in the next 


section. 
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5. Subscale Testing of Ground Motion 


Pokrovsky and Fyodorov (1969) of the Soviet Union report on ground 
shock experiments conducted in the centrifuge. These experiments used No. 8 
electric detonators buried at 9.5 cm depth in dense sand. At a modeling scale 
of 65 (@, = 65), this corresponds to 362 kg of TNT buried at 6.2 meters. They 
sought to measure peak stress jisolines using passive identation gauges. These 
"dynamometers" were 2.5 cm in diameter and 1.8 cm thick giving them a 
corresponding large size when scaled to the prototype dimensions (1.6 meter x 
1.2 meter). Hence they acted as small structures rather than free-field stress 
gauges. In addition, this type of gauge is not noted for being able to record a 
peak transient stress amplitude, to say nothing of the reported calibration 
procedure which was based upon a falling weight scheme. 

Nevertheless such a procedure does suggest the possibility of the 
direct measurement of soil-structure response to the ground shock emanating from 
an explosive in proximity to the soil. A more sophisticated structure-gauge 
needs to be devised and could take on various configurations. One of the 
simplest would be to use a metal tube with either active dynamic instrumentation 
or to just record the final deformation using a yielding type of structure. The 
important point is that the similarity relationships derived for the general 
continuum in Section 2 also apply to the subscale structure. This has not been 
exploited, but does offer a very attractive approach to subscale testing of 
structural response to ground motion. 

To test the feasibility of such a concept, a simple experiment was 
devised to examine ground motion scaling on a centrifuge. Using the concept of 
"modeling of models", two shots were designed to be scaled replicas. Two sizes 
of explosive were used, 0.49 gns and 4.08 gms of PETN, giving an order of 
magnitude yield variation. The small charge was fired at a very high value of 
gravity (504 Gs) corresponding to a very large prototype yield equal to 40 tons 
of TNT. The larger charge was fired at 294 Gs to give the same value of IT, as 
required by the similarity constraint (eq. 24). 

Both shots were half buried spheres placed in "permaplast" oi! based 
modeling clay. This material was chosen for convenience in handling and allowed 
direct scribing of a grid network on the surface of the sample. Based upon the 
expected common value for cratering efficiency, TT, (eq. 23), the scale factor 
for linear dimensions relating the two shots was 
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1/3 


1.08 | = 2.03. (45) 
0.49 


Therefore the two samples were sized accordingly. Two sets of holes were 
drilled into each sample at increasing scaled radius. The large sample had both 
a 1/2 inch diameter series and 1/4 inch diameter series. The small one 
correspondingly had a 1/4 inch diameter series and an 1/8 inch diameter series. 
The preshot sample configuration can be seen in Figs. 8 and 9. In addition to 
the holes, various rectangular and radial grids were placed on the samples. A 
row of nine sand columns was also placed along the radial bisecting the 
rectangular grid on both specimens. The charges were fired and the resulting 
crater and residual permanent ground deformation compared. As can be seen in 
Figs. 10 and 11, the expPected similarity in crater size and in the various 
measures of deformation was achieved. Preliminary analysis of the data for 
horizontal ground displacement is shown plotted in Fig. 7. The actual measured 
values at various ranges were multiplied by a factor of two approximating peak 
values as suggested by Cooper (1971). 

The results of this experiment were very encouraging indicating that 
ground motion experiments can be conducted on a centrifuge. The requirements 
imposed by boundary conditions for such response will most likely dictate larger 
sample sizes or the use of a circular segment type sample. Piekutowski (1977) 
has demonstrated, using a plexiglass boundary, about 20% of the energy is lost 
based on resulting crater size when comparing a half space results with those 
from a full space. Using a 90° segment (quarter space) with steel or tungsten 
alloy walls to minimize the symmetry boundary losses and distortion, sample size 
could be increased two or three-fold for existing centrifuges (see Fig. 12). 

Subscale model testing for ground motion response on a centrifuge also 
requires the further development of instrumentation techniques. Preliminary 
tests using carbon piezoresistive gauges to measure free-field stress history 
were performed (Schmidt and Holsapple, 1979). The results from tests conducted 
in desert alluvium and in clay demonstrated the feasibility of using miniature 
gauges on-board a centrifuge. The minor problems encountered with high 
frequency low level signal transmission across slip rings can be overcome with 
adequate signal processing. If on-board analog amplification is not sufficient, 
an analog to digital conversion can be made. A typical scheme would be to use 
an eight bit parallel transmission employing nine slip rings. A highly noise 
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Fig. 10 


Post-shot configuration of smaller ground motion specimen, 504 Gs. 
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Fig. 11 Post-shot configuration of larger ground motion specimen, 294 Gs. 
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resistant 20 MHz bandwidth could be achieved with an amplitude sensitivity of 
one part in 256. The main source of electrial noise is due to the charge firing 
circuit, which could be circumvented by using on-board capacitors to eliminate 
the high current signal through the slip rings. 

Gauges to measure structural response directly include standard 
resistive strain gauge applications, variable inductance measurements and others 
depending on specific structural configuration to be tested. In additon to 
using standard resistive type strain gauges on structural elements, a recently 
developed technique using fibre optics could be employed. This is a dynamic 
displacement gauge using two fibre light paths which can measure surface motions 
by monitoring the reflected light into the second fibre using a silicon diode or 
a photo-multiplier tube detector. This scheme is particularly attractive for 
centrifuge application in that it could all be mounted on-board the rotor and 
would not be affected by high G forces. 

In general, model building techniques need further development. There 
is first the question of modeling a given site geology. If structural response 
is to be directly measured, the structure needs to be modeled consistent with 
the scaling discussed above. Lastly, miniature instrumentation compatible with 
centrifuge environment and slip ring transmission must be developed. All of 
these elements are within state-of-the-art techniques and it is ony the special 
application that needs attention. 
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6. Summary 


The successful application of centrifuge methods to explosive crater 
formation suggests further application to ground motion studies. A similarity 
analysis confirms the applicability of previously derived scaling requirements 
to this latter problem. A non-dimensional horizontal displacement parameter was 
developed utilizing the earlier results of Cooper and Sauer (1977) in 
conjunction with a recently developed soil strength theory for cratering 
(Schmidt and Holsapple, 1979). A plot of this prediction parameter (Fig. 7) 
compares favorably with limited field data and a simple subscale test performed 
on a centrifuge. Due to the limited range of explosive yield for existing field 
data, further experimental studies are suggested using the centrifuge. This 
gravity-scaled size concept embodies equivalent yields from 1 gm to as large as 
1 Mton of TNT. The feasibility of such ground motion studies has been 
demonstrated and various schemes for instrumentation for free-field as well as 
structural response are suggested. 


90 


7. References 


Cooper, H. F. (1975) A review of ground motion from nuclear and high-explosive 
experiments, DNA Strategic Structures Review Meeting, SRI, Menlo Park, CA. 


Cooper, H. F. and Sauer, F. M. (1977) Crater-related ground motions and 
implications for crater scaling, Impact and Explosion Cratering, (D. Jd. 
Roddy, R. O. Pepin and R. B. Merrill, eds.) Pergamon Press, New York. 

Langhaar, H. L. (1951) Dimensional analysis and theory of models, p. 68, John 
Wiley & Sons, Inc., New York. 


Piekutowski, A. J., (1974) Laboratory-scale high-explosive cratering and ejecta 
phenomenology studies, AFWL-TR-72-155, Air Force Weapons Laboratory, 
Albuquerque, NM. 


Piekutowski, A. J. (1975) A Comparison of crater effects for lead azide and PETN 
explosive charges, AFWL-TR-74-182, Air Force Weapons Laboratory, Albuquerque, 


Piekutowski, A.J. (1977) Cratering mechanisms observed in laboratory-scale 
high-explosive experiments, Impact and Explosion Cratering (D. J. Roddy, R. 
Q. Pepin and R. B. Merrill, eds.) Pergamon Press, New York. 


Pokrovsky, G. I. and I. S. Fyokorov, (1969) Centrifugal Model Testing in the 
Construction Industry, Volumes I and II, "Niedra" Publishing House, Moscow, 
1969. (Draft TransTation prepared by Building Research Establisnment Library 
Translation Service Great Britain, August 1975). 


Schmidt, R. M., (1978) Centrifuge simulation of the JOHNIE BOY 500 ton cratering 
event, Proceedings of the Ninth Lunar and Planetary Science Conference, 


Geochimica et Cosmochimica Acta, Supplement 10. 


Schmidt, R. M. and Holsapple, K. A. (1978a) Centrifuge crater scaling 
experiments I: Dry granular soils, Defense Nuclear Agency Report 4568F, 
Washington D. C. 


Schmidt, R. M. and Holsapple, K. A. (1978b) Theory and experiments on centrifuge 
cratering, submitted to J. Geophys. Res. 


Schmidt, R. M. and Holsapple, K. A. (1978c) A gravity-scaled energy parameter 
relating impact and explosive crater size, (Abstract) Trans. Am. Geophys. 
Union, Vol. 59, No. 12. 


Schmidt, R. M., Holsapple, K. A. and Fisher, L. D. (1978d) 
Statistical-dimensional analysis applied to the assessment of crater 
configuration, Final Report, DNA Contract 001-78-C-0326 (in preparation). 


Schmidt, R. M. and Holsapple, K. A. (1979) Centrifuge crater scaling experiments 


II: Material strength effects, Final Report, DNA Contract N1-78-C-0149 (in 
preparation). 


91 


‘ 
a 
a 
“ 
- 
ea 
Peat 
be le TE ee eT te sete one See Pe ees awn are si ‘ 
aie cot TEER TE hey TEE Rei OEE Gt Rise SAY Mane 2 RS EE 


92 


MULTIPLE BURST HE GROUND MOTION 
GENERATION TECHNIQUES 


by 
R.G. DeRaad, Capt., USAF 
DNA Field Command 


1p 


ABSTRACT 


This paper includes a summary of the background issues generating the 
MISERS BLUFF ground motion test program and a conceptual description of the 
ground motion experiment plan. Emphasis is placed on describing those methods 
and support systems which were newly developed in fielding a multiple burst 
high explosives test program. A description is given of (1) the explosives 
quality control effort, (2) the charge placement and stacking operations, (3) the 
firing/initiation system development, and (4) the charge performance and 
simultaneity measurement program. 

The author gratefully acknowledges the support received from the MISERS 
BLUFF Test Staff and Tech. Reps., Inc. in preparing this paper. 
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MULTIPLE BURST HE GROUND MOTION GENERATION TECHNIQUES 
by R. G. DeRaad, Capt. USAF 


SECTION 1] 
INTRODUCTION 


MISERS BLUFF was a Defense Nuclear Agency sponsored series of high explo- 
sives (HE) test events to investigate ground motions generated by single and 
multiple burst detonations. The experimental program was designed to support 
the USAF Missile X (MX) system development program. The basing mode investigations 
portion of the development program is vitally concerned with the soil response or 
ground motion generated by a potential nuclear attack on MX siting bases. The MX 
system concept incorporates a multiple-aim-point (MAP) basing design which relies 
on a combination of proliferation of target points and system hardening to achieve 
a desired survivability. Figure 1-1 presents, schematically, two potential MAP 
basing concepts, the trench and the shelter, along with the potential targeting 
concepts. The MISERS BLUFF experimental program consisted of a series of HE events 
to generate a ground motion data base applicable to the MX nuclear induced ground 
motion problem emphasized in Figure 1-1, i.¢e., the problem of multiple burst 
loading of an unattacked target point from simultaneous attack on nearby target 
points. 

The ground motion experimental program, while designed to support the MX 
system issues, was built upon previous DNA sponsored ground motion studies; it used 
charges of explosive weight and type which were well represented in the existing 
data base. Important DNA tests emphasizing ground motions relatable to MISERS BLUFF 
single burst events include programs such as MIDDLE GUST, DISTANT PLAIN, Pre-MINE 
THROW IV, and Pre-DICE THROW II. The Pre-DICE THROW series with the objective of 
developing comparable trinitrotoluene (TNT) and ammonium nitrate and fuel oi] (ANFQO) 
charge explosive response, were conducted at the same site as a portion of the 
MISERS BLUFF test. However, the uniqueness of the MISERS BLUFF test series is 
that its primary objective was to develop a data base for multiple detonation 
induced ground motions. 

The following section of this report describes the general MISERS BLUFF 
program and objectives. Section 3 details the ground motion experiments. Section 
4 is concerned witn the activities and the problem areas which are felt to be 
unique to the test program, in particular the multiburst test execution. 
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SECTION 2 
PROGRAM DESCRIPTION 


The MISERS BLUFF test program was conducted in two phases and consisted 
of 10 events. Phase I, a series of eight events using small (256 and 1000 1b) 
TNT spheres, was fielded at the White Sands Missile Range, New Mexico, from 
August to December 1977. Phase II, a series of two events using 120-ton ANFO 
charges, was fielded at the Planet Ranch in Western Arizona, from April to 
August 1978. The Phase I tests, including three multiburst events of hexagonal 
array pattern, provided a large quantity of baseline data from which an analytic 
multiburst model was developed. The Phase II events, consisting of one single 
and one hexagonal 6-charge array, provided data from which analysts can refine 
and prove the model both for the larger yield and in a different (MX typical) 
geology. 

Throughout the test program, single burst events in appropriate height- 
of-burst configurations, preceded the multiburst events. The objective of these 
single burst events was to serve as a calibration test for the follow-on multi- 
burst events, and to provide a baseline measure for evaluating the adequacy of 
superposition of single burst ground motion waveforms to represent the measured 
multiburst waveforms. 

Ground motion data measurements for each event included accelerations, 
particle velocities, soil stress, and soil strain (displacement). Measurements 
were taken both in the strong motion regions and in the far field or seismic 
regions surrounding the test bed. Extensive airblast data to support the ground 
motion data were taken, particularly on the multiple burst events. Approximately 
850 channels of the ground motion and the airblast data were obtained in the 
Phase I testing and similarly 550 data channels were obtained in the Phase I] 
testing. 
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SECTION 3 
GROUND MOTION EXPERIMENTATION 


The Phase I test program, conducted at the Queen 15 site on White Sands 
Missile Range, consisted of eight separate events. Table 3-] contains a listing 
of these events along with brief objectives for each. There were essentially 
two test series within the Phase I program. One, having a half-buried charge con- 
figuration, consisted of single burst Events 1 and 3 and multiburst Event 6. The 
other, having a surface-tangent-above configuration, consisted of the single 
burst Event 2 and the multiburst Events 4 and 8. 

Varying height-of-burst configurations were implemented, including the 
fully buried Event 5, for the purpose of differentiating and amplifying specific 
aspects of the ground motion response. The surface tangent events were intended 
to emphasize the airblast induced ground motions while the half-buried and the 
fully buried events were to emphasize the crater related motions. Event 7 was 
included for yield scaling verification in the half-buried configuration. 

Explosive charges in the multiple burst test events throughout MISERS 
BLUFF were configured in a hexagonal pattern(s). This pattern was derived from 
the MX basing concept arrangements shown in Figure 1-1. The Phase I charge spacing 
distances for multiburst events were based upon the crater volume scaling 
(1/v!/3) and the anticipated MX spacing distances. The objective being to subject 
the soil within the hexagonal pattern to environment levels relatable to the 
problem of interest. 

The ground motion measurement plan for the single burst events emphasized 
the crater volume derived spacing distances and for the multiburst events 
emphasized the charge array center point which was at a spacing distance from 
any charge. The center point was emphasized; this point is the location of inter- 
est for the MX system application, i.e., the unattacked launch point which could 
be subjected to multiple, enhanced ground motions. Examples of the ground motion 
measurement plan are given in Figure 3-1, a single burst event, and Figure 3-2, 

a multiburst event. 


The Phase II test program, conducted at the Planet Ranch site in Western 
Arizona, consisted of two events. Table 3-2 contains a description, with objec- 
tives, for each. These tests were conducted in a miniature or "scaled" alluvial 
valley typical of western land areas where the MX system may be sited. In addi- 
tion to providing an opportunity to test the multiburst ground motion predictive 
model developed from Phase I data, the second event of the Phase II testing 
generated data which will be used to evaluate the "valley reverberation" problem. 
Valley reverberation is a term applied to relatively late-time low frequency 
ground motions which are anticipated to develop under the massive loading of a 
nuclear attack on a MAP system. In contrast to Phase I, the Phase II multiburst 
event spacing distance was scaled from the anticipated airblast loading at MX 
spacing ranges. This resulted in increased environment levels at the center 
of the charge array relative to either the half-buried or the surface tangent 
events of Phase I. 

The following section of this paper will treat the unique systems and 
activities essential to the multiburst testing. Included are details of the 
explosives procurement and the placement, the multicharge initiation system, 
and the means used to evaluate charge performance. 


100 


*sabueyds $0 Bulu uazno 

SnsuaA Sabueyd JSauRau adUaN{ JUL |PENLeAZ 
“p-[ JUBA WOU} e2eP YZLM UOLZELaUWU0) 

"eVep Ysungiapnu yo AZLELQLONpoOudau ssauppy 


*sabueyd qi 942 pue ql OOOL 
u9emzaq BuLl[ eds paonpuL-vazeud UO eVeP |BPLAOUd 


"UOLZeUNBLZUuod paLung-Jfey *;apow uoLrqdoLpaud 
SLSBYLUAS WAOJAARM ZSUNGLZ[NW UOJ PLP APLAOUd 


“suolLzow punoub paonpul-yseiquie azeu 
-{WL[8 PUR [PLuNg-Jo-yydap Jo yDaJya ayy yey) 


"uOLZeUNBbLJUOd 
abueyd zuahue} adejuns *apow uoLzOLpaud 
SLSOYUZUAS WAOJSARM YSUNGLZ[LNW AOL eLEP AapLAoud 


*abueuw YLLM payenuaqze 9[F7ZLL SL UOLZOW 
punoub auaymM UOLbau aoue}SLp d}eLpawuaquL 
oud S*atL SUOLbau ,Neazel[d, PayuawWnuysuL ALL NJ 


auow e BuULpnjouL 4Nq *|-] JUdAZ Jo Zeaday 


; (, de [sare,, 

aunssauduaro Aq padnpul asoy} apadaud suoLjzow 
punoubh padnpul-vazeud ay} auayM 9SOU} Jue 
sobuew Butuunszno) sabuer punoub butuunuzno 7e 
yooyus punoub ut saduavajssip [eoL6o,ouawouayd 
auayds jUaHbueZ sNsuaA paLung-J ley AJL YUapT 


‘Juaudo|laAap [apow sLsayyUAS WAoJaAeM ySuNnq 
-L}[Nw vop eTeP auLp~aseq ysung atbuis apLAoug 


uoLtydLuosaq 


"/-] JUaAZ ydadxa Sauauds IN! PuNnod-QgO, sabueyd | Ly 


Juabuey a0esuns 
patang-4 eH 
paranq-} eH 


patungq-ALLn4 


quaBue, adejsuns 


palang-J eH 


quabuel, adejsuns 


Patunq-5 eH 


qsung-J0- 346 lay 


“SLMPay IS YUSAWZ J 8SeUdq 44NIG SYASIW 


Aeuse abueyd-p2 
(41 962) 
ysung aL buls 


Aeuse Xay abueyd-9 


ysung abuts 


Aeuse Xay abueyd-9 


ysung abuts 


ysung abuts 


ysunq atbuts 
ada} 


“Le 9LqeL 


90g / 
190 92 
490 EL 


das 22 


das / 


Bny €2 


Bny SL 


Bny 2 


a7eq 


>JLON 


[=I 


oA 


S71 


v-I 


€-I 


671 


L-I 
UdAA 


101 


SS3ULS T10S Y 


(HLdjd Ly) 
NOILOW GNNOWD @ 


isvisuly OF 


39V9 NOILOW GNAOYS- + 
ON39391 


+—A-—+—+— 2) a-a-2 Bt 


INL 


“SQL OOOL 


“(T-IGW) ynoAe, abeb 4Sunq abuts ‘*T-¢ aunBbLy 


JNI7 39V9 NIWW 


NOIL03S 
i] 
0°0b - e 7 
t | 
0°0€ ——@-@ 6 
oO 
™ 
= | 
a iy 
= feo“ e @ — 
= 
. | 
& 
° 


Ove OclL 06 GLO€ 09 06 Oct OSL 
JNI1 49V9 NIVWW 
MATA NW 1d 1344 NI 41V9S 


lem 
ool ~—s 0S 0 


+ 


102 


“(p-1W) ynoke, a6e6 ysunq adiainy *z-¢ aunbl4 


206 
EMENTISVEN cy 
INBW 
isviauiy = © 
LNIWIUNSVAN gy 009 0021 
NOILOW GNNOYD 
FuvHD 
ONa931 


Oc Gl OL GS O 
SYALAW-31VIS 


o08L 
00 INIT 49V9 NIVW 
an 
eDL2 
o00E o0hZ 
W 


103 


Table 3-2. MISERS BLUFF Phase II Program 


Event Configuration Objective Date 

II-1] 120 Ton Single burst 28 Jun 78 
ANFO Stack Baseline 

II-2 Six, 120 Ton Multiburst 30 Aug 78 


ANFO Stacks 


Data Acquisition 
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SECTION 4 
MULTIPLE BURST TEST EXECUTION 
4-1] EXPLOSIVES AND CHARGE CHARACTERISTICS 


The explosives used in MISERS BLUFF Phase I testing were all TNT spheres 
while the Phase II testing employed ANFO right circular cylinders with hemisphe- 
rical caps. TNT was specified for the Phase I 1000 and 256 1b charges due to 
the existing ground motion data base for TNT spheres of that size. The large 
quantity of explosives expended in Phase II dictated use of ANFO for cost reasons. 
Also, the Pre-DICE THROW testing had established the ground motion comparability 
for ANFO and the existing data base derived from TNT and nitro-methane explosives. 

The TNT spheres were purchased from the U.S. Army Ammunition Depot, 
Hawthorne, Nevada. A total of 40 nominal 1000 1b spheres and one 256 1b sphere 
were expended in Phase I testing. Each was cast with a center located pentolite 
booster which was in direct contact with the firing detonator. Detonators were 
placed at the sphere center through an access port approximately 3/4 inch in 
diameter. Tne detonator and the det cavity were back packed with one and one 
quarter (1-1/4) lbs of C-4 explosive. Charges were cast from Grade A TNT, using 
the pelleting method. The central pentolite booster was cast from a 50/50 mixture 
of TNT and PETN also using the pelleting method. The final charge density was 
1.55 g/cc. The specification allowed a +30 1b variation for the 1000 Ib charges. 
The average charge weight for the nominal 1000 1b charges was 101] lbs with maxi- 
mum extremes of +12 1bs and -16 lbs. 

The ANFO was obtained by bid purchase from Papago Chemicals, Inc., of 
Casa Grande, Arizona. The ANFO quality control specifications were the same as 
those of the Pre-DICE THROW and DICE THROW Events. The specifications included 
allowable fuel oi] content and constraints on ammonium nitrate prill size and on 
bulk density. ANFO quality checks were performed at the test site by the U.S. 
Naval Surface Weapons Center under direction of Mr. M. Swisdak. 

Samples were taken from the bulk ANFO spread in the yoids on each layer 
of each charge while under construction (see Figure 4-1). Samples were analyzed 
for both fuel oj] content and prill size distribution. Also, the weight of every 
20th bag was checked as the charges were constructed. Table 4-1, to quantify and 
qualify variations, contains a layer-by-layer ANFO parameter record, including the 
fuel oil content for Stack II-1. Table 4-2 provides a summary of the ANFO stack 
characteristics for all of the charges of the Phase II testing. Table 4-3 provides 
a comparison of ANFO parameters from Pre-DICE THROW, DICE THROW, and MISERS BLUFF. 

The data presented in Tables 4-1, 2, and 3 indicate some of the problems 
encountered in the ANFO quality control process. A fluctuating and generally 
low bag weight necessitated continual assessment: and adjustment of the number of 
bags and the bulk ANFO to be included to obtain the proper charge dimensions and 
weight. This topic will be treated in the paragraphs on charge stacking. Four 
of the seven stacks (Stacks II-1, 2, 4, 5) had an average fuel oil content below 
the specified allowable range of 5.75 to 6.75 percent. In addition, the fuel oil 
content varied widely within a given stack as can be seen from the oj1 content 
variance given in Table 4-2. (Note the comparison to DICE THROW and Pre-DICE THROW 
in Table 4-3.) The bulk density was generally on the high side of the specifica- 
tion and in some cases exceeded the allowable range of 0.84 to 0.90 g/cc. The prill 
sizes were adequate although they contained the largest quantity of fines allowed 
in the specifications. As it became apparent from field checks that the fuel 
oil content was low, the percentage of oi] was boosted at the plant. However, 
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Table 4-1. Stack II-1] parameters. 


No. No. 
of of Total Average Total Average 
Layer Layer Whole Bags No. Bag Layer Fuel 
No. Radius Bags Bulk Bags Weight Weight 0i] 
(ft) (1b) (tons) (%) 

] 8.542 95 8 103 46.8 2.408 4.78 
2 94 8 102 46.9 2.392 4.89 
3 47.4 2.415 3.93 
4 46.6 2.379 5.59 
5 47.) 2.402 5.25 
6 47.4 2.415 4.54 
7 48.4 2.47) 4,3) 
8 48.6 2.481 5.44 
9 520 2.598 4.09 
10 47.0 2.400 6.90 
11 46.0 2.369 5.64 
12 47.0 2.397 5.48 
13 46.4 2.369 4.80 
14 45.9 2.341 5.69 
15 46.6 2.374 4.02 
16 46.6 2.379 3.68 
17 48.8 2.488 4.83 
18 47.7 2.410 5.25 
19 47.4 2.417 3.88 
20 47.7 2.433 5.99 
21 48.2 2.456 4.59 
22 51.0 2.598 5.06 
23 49.2 2.512 B7e 
24 48.2 2.458 5.45 
25 48.4 2.466 a.42 
26 47.5 2.422 5.21 
27 47 .2 2.407 5.20 
28 49.] 2.504 5.99 
29 49.0 2.501 6.03 
30 47.2 2.410 8.88 
31 46.2 2.359 6.35 
32 46.2 2.356 8.37 
33 45.3 2.310 5.54 
34 45.4 2.315 5.68 
35 44.8 2.338 6.87 
TOTAL FOR 329) 280 3571 84.752 
CYLINDER 
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Table 4-1. Stack II-1 parameters (Continued). 
No. No. 
of of Total Average 

Layer Layer Whole Bags No. Bag 
No. Radius Bags Bulk Bags Weight 

Cat) (1b) 
36 8.542 92 9 101 46.2 
37 8.53! 92 9 101 46.4 
38 8.500 9] 8 99 45.8 
ag 8.458 90 8 98 45.6 
40 8.312 88 8 96 45.5 
4] 8.208 82 9 9] 45.4 
42 8.104 82 9 91 45.1 
43 7.948 82 6 88 45.8 
44 Ye a Vs 79 7 86 45.6 
45 7.604 73 7 80 46.6 
46 7.375 68 5 73 46.1 
47 7.135 64 5 69 44.8 
48 6.875 63 3) 68 44.4 
49 6.573 60 5 65 45.3 
50 6.219 48 o 53 44.1 
5] 5.802 42 3 45 44.2 
52 5.375 36 3 39 45.2 
53 4.875 neo 2 37 44.8 
54 4.250 26 Z 28 46.3 
55 3.479 15 2 17 45.3 
56 2.417 8 - 8 a5<3 
57 1.208 4 4 45.3 
TOTAL FOR CAP 1315 117 1432 
TOTAL FOR 
CHARGE 4606 3 5003 
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Table 4-2. Charge characteristics summary. 


Average 
Total Fuel 
No. of Total 07] 
Charge Bags Weight Volume Density Content Height 
(tons) (fis) (g/cm?) (%) (0) (ft) 

II-1-1 CYLINDER 5390 83.645 2936 0.913 6.5 12.8] 

CAP 1413 34.780 128] 0.870 7.4 

CHARGE 4803 118.425 4217 0.900 6.7+1.3 21.5 
II-2-] CYLINDER 357] 84.962 2985 0.912 5.424+1.12 13.02 

CAP 1432 32.715 1140 0.919 5.99+0.98 

CHARGE 5003 117.677 4125 0.914 5 5541.11 20.06 
II-2-2 CYLINDER 3468 85.598 2936 0.934 4 .0720.95 13.02 

CAP 1358 33.610 1204 0.894 5.3521.93 

CHARGE 4826 119.208 4140 0.923 4.30+1.26 20.80 
II-2-3 CYLINDER 3366 82.223 2985 0.882 6.92+1.5] 13.02 

CAP 1463 35.897 1288 0.893 4.5040.73 

CHARGE 4829 118.120 4273 0.886 6.0441.73 21.75 
II-2-4 CYLINDER 3367 82.120 2985 0.88] 5.4120.74 13.02 

CAP 1507 35.278 1177 0.960 5.48+0.19 

CHARGE 4874 117.398 4162 0.904 5.42+0.70 19.98 
II-2-5 CYLINDER 3468 84.071 2985 0.902 5.22+1.32 13.02 

CAP 1373 33.914 1176 0.924 5.50+0.79 

CHARGE 484] 117.985 4161 0.908 5.30+1.18 20.32 
II-2-6 CYLINDER 3468 83.936 2985 0.901 5.95+0.92 13.02 

CAP 1445 33.734 1157 0.934 6.04+1.63 

CHARGE 4913 117.670 4142 0.910 5.98+1.2) 21.08 
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Table 4-3. Comparison of Pre-DICE THROW, DICE THROW, MISERS BLUFF II-1 
and MISERS BLUFF II-2 bag weights and fuel oil content. 


Fuel 
Bag No. of 011 No. of 
Weight Samples Content Samples Density 
(1b) (%) (pe) (g/cm? ) 
Pre-DICE THROW 50.8+0.8 203 6.0+0.4 5] 
DICE THROW 50.4+1.0 480 6.140.4 89 0.914 
MISERS BLUFF II-1 49.341. 198 6.72153 50 0.900 
MISERS BLUFF II-2 
STACK 1 46.9+1.9 214 5.641.1 8] 0.914 
STACK 2 49. 342.1 266 AsSt lcd 82 0.923 
STACK 3 48 842.3 25] 6.041.7 76 0.886 
STACK 4 47 .8+2.2 247 5.440.7 73 0.904 
STACK 5 48 .841.2 253 5.31.2 71 0.908 
STACK 6 47 .822.0 248 6.041.2 7] 0.910 
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a discrepancy of from one to two percent of the fuel oil content seemed to exist 
between plant and field measurements. It is presumed that the high daytime tem- 
peratures (approximately 115°F) and the two-day delay, including transit to the 
test site, caused the apparent loss of oil. In general, the ANFO quality control 
was much less precise than that achieved in previous tests such as the DICE THROW 
Events. 


4-2 CHARGE PLACEMENT AND STACKING 


In Phase I testing, a charge placement technique was developed which 
allowed intimate contact of the half-buried and the fully buried spheres with an 
undisturbed soil cavity. A template was used to bore a-hole of the proper diameter 
(see Figure 4-2). The sphere was lowered into place with a lifting device which 
gripped the interior of the detonator access port. A small amount of soil matching 
grout was placed in the soil cavity to provide intimate contact at the soil/charge 
interface. A half-buried TNT charge is shown in place in Figure 4-3. 

The ANFO charge stacking process did not involve the placement problems 
presented by the TNT spheres but was complicated by variations in ANFO bag weight 
as mentioned earlier. The Event II-1] stack, as the result of an average bag 
weight of 49.3 lbs, was completed at a weight of 118.2 tons versus the design 
weight of 119.9 tons. However, it was very close to the design physical size of 
21.5 feet high with a radius of 8.5 feet and contained the design 4800 bags. A 
photo of the completed stack is given in Figure 4-4, For Event I1-2, the supplier 
was requested to provide a similar bag weight in order to obtain six stacks as 
identical as possible to the first event stack. However, bag weights tended to 
be lighter and fluctuated. It was necessary to compensate by adjusting the number 
of bag layers and the amount of bulk ANFO filling the voids. The adjustment 
approach used was: (1) determine the average bag weight as the cylinder was 
nearing completion and, (2) adjust the amount of bulk in the cylinder cap to pro- 
vide the needed compensation. Extra bag layers were required to obtain the desired 
cylinder height for Stack ] (2 layers), and Stacks 2, 4, 5, and 6 (1 layer each). 
This too resulted in weight variations which required compensation by varying 
the amount of bulk ANFO in the cap. The adjustment process was a dynamic effort 
which resulted in Event II-2 stacks which deviated from Event II-1] weight by less 
than 1 percent each. 

The stacking process for the six stacks of Event II]-2 required 11 days. 
Stacking was begun on 14 August and was completed 24 August, six days before event 
day. Each stack was constructed in a wooden shelter as seen in Figure 4-1, which 
gave protection from the elements and provided a form for the stacking process. Two 
stacking crews were employed to erect the charges. Charges 2 and 5 were completed 
on 15 August; charges 1 and 4 on 21 August, and charges 3 and 6 on 24 August. Note 
that opposite charges in the array were constructed at a given time so that one sec- 
tion of the array was not “aged" more than another. Figure 4-5 shows the completed 
charge array for Event II-2. 


4-3 ARMING AND FIRING SYSTEM 


The timing and firing system requirements specified for the MISERS BLUFF 
program were quite stringent. In the initial phases of program planning, the 
Data Analysis Working Group (DAWG) sponsored by Hq DNA, determined that simul- 
taneous detonations in a hexagonal pattern would provide the necessary data for 
the time related ground motion analysis at symmetry points involving two, four, 
and six charges. Further, nonsymmetric measurement points would allow analyses 
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Figure 4-3, Half-buried TNT charge. 
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Figure 4-4, Complete ANFO stack. 
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simulating non-simultaneous detonations. For these reasons and for simplicity 
the simultaneous detonation was specified. 

For practical, fielding purposes it was necessary to specify an allowable 
detonator firing time spread for the multiple burst initiation. Studies, reported 
by the DAWG, indicated that the time-on-target sensitivities for ground motion 
enhancement in the MX environment could be specified and that a standard deviation 
in burst timing of 0.1 seconds is essentially simultaneous. Time scaling to the 
smal] yield of Phase I would indicate essential simultaneity can be attained from 
detonation variations on the order of 1.0 to 0.1 ms. 

The Sandia Laboratories, Albuquerque (SLA) was selected to provide the 
arming and firing system for all MISERS BLUFF Events. The system was procured 
under specification that the maximum allowable time spread for multidetonator 
firing was not to exceed 10 us. A schematic of the firing system is shown in 
Figure 4-6. The schematic is representative of all Phase I Events and, with 
variations, of Event I-8 (the 24-charge event) and the tests of Phase IJ. The 
firing unit, a TC-132 capactive discharge unit designed for nuclear field test 
applications, is used to fire TC-246 exploding bridgewire detonators. A TC-369 
transverter charges the x-unit to 3 kV, which when triggered, discharges through 
a load ring into the detonator firing cables. 

Variations of the basic firing system were used in Event I-8 and in Phase 
II testing. In Event I-8, the TC-132 x-unit was used as an intermediate trigger 
to operate four TC-406 x-units. Each TC-406 unit then fired six TC-246 detonators, 
j.e., six of the TNT spheres. The TC-406 x-unit is identical to the TC-132 except 
for the triggering circuitry necessary to accommodate the output of the TC-132 
when used as a trigger. In Phase II, Event II-1, the TC-132 x-unit was used to 
fire seven TC-234 detonators which were located on the centerline of the ANFO 
cylinder. The TC-234 is the same detonator as the TC-246 except that there is an 
additional three and one-half (3-1/2) gram pellet of RDX explosive on the end of 
the detonator. The additional boost was desired, as an added safety margin, to 
initiate the pentolite boosters of the Naval Surface Weapons Center designed 
booster initiation system. For Event II-2, the TC-132 x-unit was again used as a 
trigger unit to operate six TC-406 x-units, one for each ANFO stack. Each TC-406 
unit, in turn, fired seven TC-234 detonators. 

Certification tests were performed by SLA to prove the system design for 
each configuration prior to fielding. The maximum time spread obtained for any 
of the system configurations was 49 ns with a simultaneity sigma of 27.1 ns. These 
results are well within the prescribed allowable 10 us spread time. 


4-4 SYSTEMS TO MONITOR CHARGE PERFORMANCE/SIMULTANEITY 


Optical systems were used as the primary means to measure the charge 
performance and the simultaneity in both phases of the MISERS BLUFF testing. The 
photographic and photometric diagnostic instrumentation consisted of (1) high-speed 
Dynafax cameras, (2) Unit Light Radiation (ULR) devices, and (3) streak cameras. 

Light breakout at the surface of the charge was used as the measure of 
simultaneity of detonations in multiburst events. The ULR devices, silicone solar 
cells with rise time of I us, were used to monitor light breakout at a selected 
point on each charge in the multiburst events. Data from the ULR devices were 
recorded by oscilloscope mounted cameras. These measurements were supported 
with coverage of all detonations by Dynafax cameras operating at 26,000 frames 
per second, i.e., an interframe time of approximately 40 us. ULR radiation 
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device data from multiburst Events I-4, ]-6, and 1-8 show a maximum charge breakout 
time spread of 6.6 us, occurring in Event I-4. Due to a limited supply of ULR 
devices (7), only seven charges were monitored on the 24-charge Event, 1-8. Remain- 
ing charges were monitored by the Dynafax cameras which could give a maximum time 
resolution of 40 us. Final timing data from multiburst Event II-2 are not 

available at this writing. 

For the Phase II] ANFO events, an additional diagnostic tool was added to 
those described above. A light transmission pipe system, developed by the Denver 
Research Institute, was installed in each ANFO stack to monitor the performance 
of the main booster assembly system. The light transmission system consisted of 
metal pipes, with polished interior, which were abutted to each of the seven octol 
boosters located on the centerline of the charge cylinder (See Figure 4-7). The 
light pipes provided a means of transmitting light generated by the booster detona- 
tion to the exterior of the stack where streak cameras (at a distance) were employed 
to record the time related appearance of light from the seven boosters. It is 
anticipated that these records along with the ULR data will provide an average ANFO 
detonation velocity for each stack in addition to simultaneity data. 

Much data have been gathered to document the MISERS BLUFF explosive 
charge preparation and the execution. The final evaluation of the charge perfor- 
mance characteristics must await the compilation and review of the data and 
preparation of the several reports. 
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SECTION 5 
CONCLUSIONS 


This paper has presented in brief summary the background issues giving 
rise to the MISERS BLUFF multiple burst ground motion test program. The ground 
motion measurement program was described in concept. Emphasis has been placed 
on describing those methods and support systems which were newly developed in the 
fielding of a multiburst test program. Unique problems were encountered and solu- 
tions were developed in the areas of (1) explosives quality control, (2) placement 
and stacking of multiple TNT charges and ANFO stacks, (3) firing/initiation systems 
for a variety of multiple charge explosives and, (4) methods to monitor the charge 
performance and the simultaneity of multicharge events. 

The MISERS BLUFF test program was highly successful in its primary ob- 
jective: that of developing a multiple burst ground motion data base for high 
explosive detonations. The data are currently undergoing evaluation in an attempt 
to understand the phenomena generated by the multiple energy sources interacting 
in the nonlinear regime. Observations thus far indicate that additional investi- 
gation may prove useful. They may include specialized small scale testing, i.e., 
1000 1b tests, designed to emphasize particular aspects of the ground motion wave- 
. forms through the variations in height-of-burst coupled with the variations in 
charge spacing distances. Laboratory testing may provide useful data on spall 
related phenomena and the convergence of horizontal waves at the center of multi- 
Source arrays. In addition, the information from MISERS BLUFF has emphasized the 
need to review and to reevaluate results from the previous testing. 

It is clear that HE multiple charge ground motion data are but a first 
step toward evaluating the MX nuclear ground motion problem. The need to establish 
the nigh explosive nuclear ground motion link is a very significant problem re- 
quiring further evaluation. Analytic evaluations and the comparison of existing 
data will be a part of the continuing effort. Additional testing and tne develop- 
ment of a high explosive charge more nearly duplicating the energy partitioning 
and the energy coupling characteristics of the nuclear event are likely. Much 
work remains to be accomplished on the ground motion issues. The MISERS BLUFF 
test program has provided a baseline which may be built upon through both analysis 
and test. 


120 


ACKNOWLEDGEMENT 

This paper has drawn upon information contained in proposals and pre- 
liminary data reports from participants in the MISERS BLUFF Test Program, 
particularly those of the Denver Research Institute, the Sandia Laboratories, 
Albuquerque, and the U.S. Naval Surface Weapons Center. 
REFERENCE 


Data Analysis Working Group, "A Review of High Explosive Testing to Investigate 
Ground Motions Pertinent to the MX Multiple Aim Point Systems," May 13, 1977. 


121 


toa tee 


Spot a see 


122 


BLEST — A GROUND SHOCK SIMULATOR 


by 
Frank E. Seusy 
Air Force Weapons Laboratory 
Harry E. Auld 
Civil Engineering Research Facility 


lah 


BLEST - A Ground Shock Simulator 


DEVELOPMENT 


The Berm-Loaded Explosive Simulation Technique (BLEST) was developed during 
the HARD PAN Test Program to provide a low-cost method of simulating airblast 
loadings over a large area. The requirement to load a large area was generated 
because the HARD PAN Test Site geology consisted of rock layers at relatively 
Shallow depths and it was predicted that the ground shock would outrun the air- 
blast environment. The BLEST design approach assumes that at large distances 
from a target, it is relatively unimportant if the correct pressure history is 


simulated so long as the correct total impulse is applied with the proper timing. 


The BLEST was initially conceived as an explosive weave or sheet overlaid 
with soil (fig. 1). Feasibility experiments were conducted using a weave of 
detonating cord overlaid with a soil overburden. The second concept was a matrix 
of charges placed in drill holes which approximated a plane of explosives (fig. 
2). There were substantial cost advantages with this approach. and it was selected 
for further development. In order to standardize on a limited number of charge 
weights, it was decided to vary charge spacing as a means of varying the effective 
areal charge density. A standard charge depth of 1.8 yes was selected where 


Wis the charge weight in pounds. 


Calibration tests are required to determine impulse as a function of charge 
spacing for a given explosive and for the soil at a particular test site. The 
impulse delivered to a test plane can be determined by using photopoles to map 
the vertical velocity of the overburden around charges and by knowing the over- 
burden density. Active instruments, including soil stress gages and acceler- 


ometers, are more expensive and are thought to be less reliable than photopoles. 
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Figures 3, 4, and 5 are BLEST calibration data presented in three different 
ways. Figure 3 is data from the HARD PAN Test Site with scaled impulse plotted 
versus scaled charge spacing. Figure 4 is impulse plotted versus charge spacing. 
The solid lines are derived from the straight line fit to the data in figure 3, 
with individual data points from HARD PAN overplotted. The dashed lines on 
figure 4 show the limits of charge spacing in terms of charge depth that have 
been tested. Figure 5 is impulse plotted versus the square root of the average 
charge weight per ft? times the overburden weight per #2. Figure 5 includes 
data from experiments configured as shown in figure 1 at the McCormick Ranch 
test area as well as data from experiments configured as shown in figure 2. 

No explosive energy corrections were made to the data in figure 5. The explo- 
Sive used on the HARD PAN experiments was an aluminized ammonium nitrate slurry. 
Detonating cord (PETN) was used on the experiments at McCormick Ranch. The 
important point to recognize in figure 4 is that very Jow impulses (<7 psi-s) 
require discrete charges of less than 25 lbs, which is near the lower bound for 


many low-cost slurries. 


The need for loading large surface areas was identified by studies of 
ground shock arrival times on vertical buried structures in layered media. The 
design of load areas and the design of impulse levels within the load area re- 
quires data on the seismic profile of the test bed, the location of target points 
of interest and the required simulation times at those target points, and the 
ground zero, height of burst and yield of the weapon to be simulated. The load 


area is bounded by the locus of points described by: 


o 
th 


Time of first arrival at target point 
Ty = Desired simulation time 
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AS 


GS 


= Time of arrival of the air shock at a specified point in the 


test bed 


= Travel time of the ground shock from the specified point in the 


test bed to the specified target point. 


The time of arrival of the air shock from nuclear detonations has been documented 


in a variety of forms. One which has been used extensively is given by reference 


1 


where, 
W 
R 


Tas 


if 


tt 


_ .5429W - 21,185RNe/? + 361.9R° Wl’? 
173 


we/3 + 2.04gRW'/3 + 2.687R° 


3 
+ 2383R 
Tas (2) 


2x weapon yield in KT for surface burst 


Range in K ft 


= Arrival time in ms 


The minimum travel time of the ground shock (compression wave) between two points 


on the surface of a layered geology is determined by evaluation of the following 


relation (reference 2): 


N 
N 1 ] X 
ie 3 2 Hi Va - one |e (3) 
1 i oN 


C 
wa N+1 


H. = Layer thickness, ith layer 

C. = Compression wave velocity, ith layer 
x = Range between two points 

N = Number of possible paths 


1. Brode, H.L., Height of Burst Effects at High Overpressures, RM6301, The Rand 
Corps., Santa Monica, California, July 1970. 


2. Crawford, R.E., Higgins, C.J., Bultmann, E.H., The Air Force Manual for 


Design and Analysis of Hardened Structures, AFWL-TR-74-102, Air Force Weapons 
Laboratory, Kirtland AFB, New Mexico, October 1974. 
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With definitions of the time of arrival of the air shock and of the travel 
time of the ground shock, the minimum value of the sum of the two for a given 
target range is the value for Ty> the time of the first arrival at the target 
point. The load area for a hypothetical case might appear as shown in figure 6 


for load points A, B, and C. 


At the boundary of the load area, the load is applied for zero time; i.e., 
there is no excess time for a portion of the pressure time history to reach the 
target point. Interior to the load boundary there is excess time, Ty > for the 
load to be transmitted to the target point. Interior to the load area boundary, 


equation (1) becomes: 


oe 


mie hag. © Tes (4) 
The interpretation of T used in the HARD PAN series is illustrated in figure 
7. If the pressure and impulse time histories shown in figure 7 are the ideal 
nuclear overpressure waveforms, and if the ground shock traveled at a constant 
velocity C, the target point could only receive the portion of the signal sent 
during Te On this basis, I, would be the design impulse at that point. Pro- 
grams have been written for the HP9820 to draw contours of constant impulse for 
design of a BLEST. The inputs to these programs are polynomial or segmented 


fits to the air shock time of arrival, ground shock transit time, and nuclear 


impulse as a function of range and time. 


There are non-simulation effects ieeoctatee with overburden fallback from 
the BLEST. Essentially the same impulse is delivered to the test bed by the 
fallback as was initially imparted to the test bed and overburden. The use of 
discrete charges and relatively high overburden velocities tends to scatter the 
overburden in space, and delay the return in time. Therefore, overburden fallback 


cause little concern in most test applications. 
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Figure 6. Load Area for Hypothetical Case 
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Figure 7. Interpretation of Ty 
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For small areas of relatively low impulse, it iS less expensive to con- 
struct a BLEST as shown in figure 1. The construction technique shown in figure 
2 was used in the HARD PAN series, and appears to be of more general interest 
for future applications, and therefore cost estimates for that technique will be 
presented here. The basic data used for estimating costs are shown in table 1. 


Estimated cost per square ft versus impulse are shown in figure 8. 


Table 1 
BASIC DATA FOR BLEST COST ESTIMATES 


: Drilling Hole Chg Cost Over 0.B. Cast 


Chg Wght Hole Diam Cost/ft Survey 1b Burden D yd? 
25 1.0 ft $5.00 3.60 .50 2 ft $2.00 
50 1.0 5.00 " n i i 
100 1.5 5.40 " ul i ‘ 
200 2.0 5.70 " " h ‘ 
400 3.0 6.20 H u 7 ‘ 
800 4.0 6.69 " " ; ' 


APPLICATION AND ANALYSIS 


The HARD PAN I Test Series was conducted in a layered geology near Trading 
Post, Kansas. Combined High Explosive Simulation Technique/Berm Loaded Explo- 
sive Simulation Technique (HEST/BLEST) test configurations were utilized in three 
of the tests (Events 2A, 2B, and 3) to simulate an environment of nuclear ground 
motions in a layered geology. The HEST was designed to simulate the effects of 
the propagating local airblast and the BLEST was designed to simulate the effect 
of upstream airblast. A scaled test structure was buried at a specified location 
under each HEST cavity and this structure was the target of the simulated nuclear 
blast. Bath the scaled structure and the free-field beneath the HEST cavity 


were heavily instrumented to measure their motion and stress response to the 
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simulated nuclear environment. Each of the three events was designed to simulate 
a peak overpressure of 1200 psi at an equivalent nuclear range of 725 ft from 

a 125-kt yield detonated at a 210-ft height-of-burst. Figure 9 shows the size 
and shape of the three HEST/BLEST events, and table 2 contains some additional 
facts about the details of the experiments. For a more complete description, 

see reference 3 and 4. 


Table 2 
HARD PAN I Events 


Dimensions, ft Simulation 


Tinie. ** 
Event | HEST Array | BLEST Array* ms 


64 x 96 158 x 222 19 Dec 74 
88 x 88 256 x 316 24 Aug 75 
120 x 200 450 x 550 4 Dec 75 


Notes: * The BLEST arrays were oval. The dimensions were the maximum spatial 
extent of the discrete charges in directions perpendicular and parallel 
to the long axis of the HEST area. 

**x Simulation time refers to the time to which the experimental pressure 
history is to match the design pressure history. 
The simulations of the HARD PAN I Test Series did not produce crater--or 
direct-induced effects. This permitted the exclusive study of airblast-induced 
effects. However, even with this simplification, two types of airblast-induced 


ground motions are possible in layered media, viz., those resulting from local 


airblast effects and those resulting from upstream airblast effects or outrunning. 


In conventional definitions of outrunning, such as that given in reference 


3. Doran, J., HARD PAN I Test Series and Instrumentation Plans, AFWL-TR-75-249, 
Vol. 1, Air Force Weapons Laboratory, Kirtland AFB, New Mexico, December 1975. 


4. HARD PAN I Test Series Design and Construction of Test Facilities, AFWL-TR-/6- 
60, Air Force Weapons Laboratory, Kirtland AFB, New Mexico, June 1976. 
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5, the outrunning region is considered to commence at a range where a headwave 
from a deeper refracting layer emerges at the ground surface simultaneously with 
the arrival of the airblast loading (fig. 10). Reference 2 presents a technique 
for calculating this range (designated R2 in figure 10). It should be noted that 
outrunning may commence while the airblast shock front velocity, U, still exceeds 
the seismic velocity of the upper material. In this paper, the definition of out- 
running is extended to include any ground motion resulting from a source other 
than the local airblast loading. This expanded definition of outrunning permits 
the contributions from upstream effects to be included regardless of their time- 


of-arrival with respect to the local airblast-induced ground motion. 


Until fairly recently, outrunning ground motions were considered to be prim- 
arily associated with ranges corresponding to fairly low overpressures. This 
misconception may have resulted from the unique geology present at the Nevada Test 
Site, where most of the nuclear test ground motions data have been obtained. 

It is now known that outrunning ground motions can occur at high overpressure 
levels, so long as the geology consists of high seismic velocity layers located 
near the ground surface. HARD PAN I produced significant outrunning ground motions 


at the 1200 psi overpressure level. 


Appendix A to reference 6, documents computer calculations of Event 3 which 
were performed to estimate the magnitude of the separate effects from the BLEST 
loading alone (designated BLEST only) and from the HEST loading alone (designated 


HEST only). The problems were calculated with TOODY II-A (ref. 7), a finite 


5. Sauer, F.M., et al, "Empirical Analysis of Ground Motions and Cratering (Part 
Four)," Nuclear Geoplosics, DASA-1285(IV), Defense Atomic Support Agency, 
Washington, D.C., May 1964. 


6. Auld, H.E., et al, ANALYSIS OF HARD PAN I FREE-FIELD DATA, Vol. 1: Summary of 
Test Data, AFWL-TR-78-66, Air Force Weapons Laboratory, Kirtland AFB, New 


Mexico, in publication. 

7. Bertholf, L.D., Benzley, S.E., TOODY II-A Computer Program for Two-Dimensional 
Wave Propagation, SC-RR-68-41, Sandia Laboratories, Albuquerque, New Mexico, 
November 1968. 
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difference code that uses Lagrangian equations of motion to calculate two-dimensional 
wave propagation in layered media. The typical layered geology of the HARD PAN 

I test site that was utilized in the calculations is indicated in figure 11. 

The water table is located at a depth of 8 ft’ and high seismic velocity refracting 


layers (C > U)* are located at depths of approximately 17, 32, and 52 ft. 


Results of the calculations are shown in figure 12 for the structure range 
of 725 ft. The depths of 14, 30, and 56 ft were chosen to show the extent of 
the variation of the BLEST contribution. The outward component of the BLEST had 
a more significant relative effect than the vertical. In regions other than the 
outrunning limestone layer (56 ft), the BLEST had little effect on the first out- 
ward peak due to the HEST, but it increased the outward velocity significantly 
after this peak. In the outrunning limestone layer the BLEST-induced first arrival 
and the peak outward velocity associated with the BLEST were approximately four 


times the peak associated with the HEST. 


Comparisons of these numerical simulations with the experimental data from 
the HARD PAN I Test Series reveals that they did not accurately predict the 
magnitude of the measured ground motions, particularly the horizontal component. 


However, the calculational trends were useful in interpreting the experimental data. 


Typical outrunning ground motion records for Event 2B are shown in figure 
13 for a depth of 48 ft. At this depth and range the initial motion results from 
outrunning. The initial vertical motion may be either up or down, depending on 
the local geology and whether the depth of interest is above or below the refracting 
layer generating the first-arriving headwaves. The initial horizontal motion 


will always be outward regardless of the location with respect to the layer. 


* C = seismic velocity of refracting layer (compression wave). 
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Dry Clay 125 C235 


Wet Clay 128 5,000 0.40 


Limestone 160 10,400 0.30 
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Limestone 160 10,400 0.30 -1000 
160 9,100 0.40 -36 
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Limestone 90 13,300 0.30 -1000 


Shale 150 10,400 0.30 -1000 


Typical Geologic Profile, HARD PAN I Test Site 
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Figure 12. HEST and BLEST Contribution to Velocity Waveforms at 
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Figure 13. Typical Qutrunning Ground Motions in Event 2B 
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Figure 14 illustrates the theoretical situation in which a layer with a high 
seismic velocity is bounded on either side by materials with lower seismic veloci- 
ties. Three types of headwaves will be generated--a compression wave (P-wave), 

a von Schmidt wave (PS-wave), and a shear wave (S-wave). The equations for cal- 
culating the propagation angles and thereby the direction of initial motions for 


each wave are indicated in figure 14, 


Upper and lower bounds to the measured peak particle velocities for each 
event are shown in figure 15. The upper bound. to the peak vertical velocity 
(fig. 15a) at the surface is nearly coincident for all events, although the lower 
bounds at the surface are more divergent. In near surface layers, the local air- 
blast (HEST) dominated the vertical motion histories. At greater depths, the HEST- 
induced particle velocities attenuated relative to the surface and outrunning 
(BLEST or upstream HEST-induced) motions became more important. The lower bounds 
at depth are nearly equal in all events, but the upper bound for Event 2A is sig- 


nificantly greater than the upper bounds for Events 2B and 3. 


Bounds to the peak horizontal velocities are Shown in figure 15b. The 
‘most obvious feature is the great difference between Event 2A and Events 2B and 
3. For the geology at this test site, the peak horizontal velocity due to the HEST 
should increase to a depth of about 10 ft and then should decay rapidly until it 
becomes essentially zero below a depth of about 60 ft. The upper bound follows 
this pattern in the near-surface material, but it does not decay to zero. This in- 
dicates that outrunning had a large effect on the magnitude of the horizontal 


velocity below the near-surface material. 


There is no generally accepted technique available for accurately predicting 
outrunning ground motions. However, empirical prediction techniques based on the 


available high-explosive ground motion data for the outrunning region are given 
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Definitions: 
C = propagation velocity Cp, > Choe Co, 


P = compression wave Cs, >C. , Cy 
PS = von Schmidt wave 

S = shear wave 
— > indicates direction 


—p of first motion a ee Cs, 1€. 
sinw. = Cs fen, 
sing, = €, /C. 


Figure 14. Headwaves in Layered Medium 
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Upper and Lower Bounds to Peak Velocity Data (2 of 2) 


_ Figure 15. 


in reference 8 and for the available nuclear data are given in reference 5. 


On the basis of a thorough analysis of the observed experimental free-field 
data, estimates of outrunning motions based on the available prediction techniques, 
and the results of the finite difference calculations, it can be concluded that 
the outrunning ground motions in Events 2B and 3 were adequately simulated whereas 
in Event 2A they were over simulated; perhaps by as much as a factor of two. 

It can further be concluded that the BLEST generated an impulse in excess of the 
required amount in Event 2A. This thesis is supported by an analysis of the limited 
photopole and acceleration data that were available in the BLEST region of the ex- 
periment (ref. 9). After Event 2A, a BLEST calibration series was conducted at 

the HARD PAN Test Site using the same slurry used for Event 2A and subsequent events. 
The data from a previous test series at McCormick Ranch, which was used to design 
Event 2A, was suspect because of differences in the soil properties and minor 
differences in the explosive slurry. A consistent method for mapping overburden 
velocities was not developed during the McCormick Ranch series, which may have 

also been a portion of the problems. All HARD PAN BLEST data following the on-site 
calibration series was consistent. It has been concluded that BLEST arrays should 


be calibrated on site for future field tests or experiments. 


POTENTIAL APPLICATIONS 


The BLEST is the most energy efficient simulation technique currently avail- 
able for application of surface impulsive loads. It cannot be used in the immed- 


jate vicinity of a test article because of the high stress levels which are produced. 


8. Higgins, C.J., and Schreyer, H.L., Analysis of Qutrunning Ground Motions, 
AFWL-TR-74-220, Air Force Weapons Laboratory, Kirtland AFB, New Mexico, 


May 1975, 


9. Webster, J.B., I]1, Maj., and Suesy, F.E., “Simulation Performance for HARD 
PAN I, Event 2, Memorandum for the record, 22 January 1975. 


However, valid simulation can be provided in cases where sufficient distances are 


available between the explosive plane and the test article. 


Several applications of BLEST have been proposed and are currently under 
study. The current basing concept for the MX missile system envisions a large 
number of vertical shelters, a fraction of which are occuppied by missiles. The 
shelters surviving an attack may need to survive multi-burst effects, i.e., the 
effects from nearly simultaneous detondtions on neighboring shelters. Because of 
the large spacing between shelters, it has been postulated that outrunning ground 
motions may be an important effect. OQutrunning ground motions are not well] under- 
stood from single bursts, and less well understood for multiple bursts. The BLEST 
offers potential for experiments at sufficiently large scale to increase under- 
standing of both the single and multi-burst case. A major reason for large scale 
experiments is to permit experiments to be run on realistic geologies. The 
BLEST may not be technically or economically feasible at full scale, and other 


testing techniques may be required. 


Deep underground missile basing concepts are also under study. The BLEST 
may also have potential for ground shock environment studies and small scale 
testing. Again, full scale may be technically and economically impractical be- 


cause of the extremely large amounts of energy required. 


POTENTIAL IMPROVEMENTS 


Earlier in the paper it was stated that photopoles were the current tool 
for measurement of the impulse delivered by a BLEST array. While photopoles 
have performed better than accelerometers and stress gages, they are not with- 
out shortcomings. One of the basic difficulties with photopoles is that con- 
siderable care must be taken to avoid obscuring the field of view with over- 


burden from other portions of the test bed. The photopole itself is simple 
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and inexpensive, but the data reduction process is fairly expensive and time 
consuming. Some effort will be expended in the next year in an attempt to iden- 


tify potential alternative measurement techniques. 


The calculations for loading area and loading time have been accomplished 
in the past with seismic velocities. This may be overly conservative, in that 
the important effects are probably better defined by peak propagation velocities. 
Peak propagation velocities have not been used in the past because of the obvious 
increase in design complexity. A useful study of load area and time would most 
likely include three-dimensional calculations and supporting field experiments. 
A substantial effort would be required if significant geologic variations were 


included. 


The current recommendation that the BLEST be calibrated on site with the 
intended explosive may be excessively conservative. A theoretical model verified 
experimentally may allow adequate prediction from the soil density and the 


available explosive energy. 
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I. INTRODUCTION 


The Direct-Induced High Explosive Simulation Technique (DIHEST) is a 
method of simulating ground shocks which result from the direct coupling 
of energy from large explosions. The original weapons simulation work was 
accomplished by the Air Force Weapons Laboratory (AFWL) and the University 
of New Mexico Civil Engineering Research Facility (UNM/CERF) with funding 
from the Defense Nuclear Agency (DNA). Recent earthquake simulation work 
has been accomplished at UNM/CERF with funding from the Electric Power 
Research Institute (EPRI) and the National Science Foundation. 


A typical DIHEST array consists of a row of vertical dril? holes in 
the ground which are filled to some height with explosives and then stemmed 
with sand or grout (Fig. 1). Multiple detonators are placed in each-hole 
and fired simultaneously to create a shock wave which is approximately 
planar. DIHEST arrays have been fired with total explosive charges as 
large as 40 tons. Array dimensions have been as large as about 300 ft deep 
and over 1100 ft long. 
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II]. DIHEST CHARACTERISTICS 


The peak ground motions resulting from a DIHEST array are expected to 
attenuate with range approximately as shown in Figure 2. At ranges near the 
array, the relief effects do not affect peak parameters of acceleration (a), 
velocity (v), and displacement (d); and the attenuation is planar and due only 
to inelastic effects. Beyond some range, relief waves from the top and 
bottom of the array begin to overtake the loading waves and cause the 
attenuation rates to approach cylindrical conditions. Finally, at some 
greater range, relief waves from the ends of the array overtake the loading 
waves; and a transition to spherical attenuation might be expected (although 
few DIHEST data have been recorded in this region). The transition regions 
for the parameters of acceleration, velocity, and displacement are expected 
to be different. The accelerations transition over a short range, the dis- 
placements over a longer range, and the velocities somewhere in between. 

The transition region between two zones has been taken as a point, called 
the transition-point range, Ro 
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Log-Scaled a, v or d 


Planar 
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Cylindrical 
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Spherical 
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Log-Scaled Range 


Figure 2. Effect of Finite Array Dimensions on Peak 
Ground Shock Attenuation 
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III. GROUND MOTION PREDICTION 


The major ground motion quantities of interest are peak horizontal 
acceleration, particle velocity, and displacement. In the absence of 
material property variations, the parameters which control DIHEST ground 
motions are: 


Areal Charge Density, a 
Array Height, H 
Array Length, L 
Surcharge Above the Explosive, S 
Range to the Point of Interest, R 
The surcharge term, in effect, is a measure of the location of the 
free surface. Gravity may also be influential, but there are insufficient 


data available at present to understand its effect. 


Dimensionally consistent scaling relations between the ground motions 
and the array parameters are: 


= Fy (R/a. H/o, L/Hy R/S) 


om 
Q 
I 


Fe (R/a, H/o, L/H, R/S) 


< 
{l 


d/a = f3 (R/a, H/o, L/H, R/S) 


It has been observed experimentally that if L is greater than about 2H and 
if R is less than 2L, the effect of L can be assumed to be small. Also, 
if observations are restricted to the mid-depth of the array, the effect 
of S can be assumed to be small. With these assumptions, the relations 

at the mid-depth of the array are: 


a-a = f; (R/a, H/a) 


v= fz (R/a, H/a) 


d/a = f3 (R/a, H/a) 


The scaled motion quantities are usually plotted versus scaled range, R/a. 
The effect of H/a is evident from the relative positions of the results 
from individual experiments with varying H/ca. 


The DIHEST data available prior to the recent SIMQUAKE I test series 
are found in References 1, 2, 3, 4, 5, and 6. These data are analyzed in 
Reference 7. The prediction relations for motion amplitudes at the mid- 
depth of DIHEST events, which were derived from these data and a large 
number of finite difference calculations (see Ref. 7), are as follows: 


157 


asa = 2900 (R/a) ~2-89 for R/a < 2(H/a)%*® 
aso = 11,680 (H/a)t*4(R/a)2°? ~— For Ryo > 2(H/a)?°® 
v= 49 (R/a) 72-18 for R/a < 2.6(H/a)9*>2 
v= 313 (H/o) (R/a) 7"! for R/a > 2.6(H/a)°*>¢ 
d/a = 10.5(R/a) 2°24 for R/a < 1.15(H/a)?°/? 
d/a = 13.3(H/a)'*®8(p/o)7t-8* for Ryo > 1.15(H/0)9°” 
where 


a = acceleration in g's 
v = particle velocity in ft/s 
d = displacement in in 
H = array height in ft 
R = range in ft 
a = areal charge density in 1b/ft® of TNT 
Three important points should be made with regard to these relations. 


First, the normalized transition range, R O/H is dependent upon H/a. 
The dependencies are: 


Acceleration: RA/H = 2.0(H/a) "9°" 
Velocity: R/H = 2.6(H/a) 0°48 
Displacement: Ro/H = 1.15(H/a)~?°2° 


These relations are not evident in the data alone because of scatter, but 
they are derived from finite difference calculations. 


The second important point jis that shallow events (array mid-depth 
within about 15 ft from the surface) produce greater scaled horizontal 
displacements than deep events. The reason for this behavior is not 
clear, but material property differences near the surface and/or gravity 
are possible causes. The displacement equations given above are for deep 
events. Displacement from shallow events will be about twice those 
estimated from the given equations. 


Third and last, most confidence can be given to the equations beyond 


the transition range. The behavior in the planar region remains more 
uncertain because very few data are available in this region. 
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The SIMQUAKE I tests (Ref. 8) were the most recent DIHEST events 
in soil. They were conducted at the McCormick Ranch test site of the 
University of New Mexico in May and June 1977. The main objective of 
the tests was to excite models of nuclear power plants with horizontal 
ground motions having amplitudes and frequency contents on the order of 
those expected in very strong earthquakes (Richter Magnitude of 7 to 8). 
Figures 3 and 4 show plan and elevation views of the test site. The test 
was originally to consist of two arrays 75 ft deep by 200 ft wide with 
25-ft surcharge. The back array, designed to fire first, was loaded with 
6.28 1b/ft? of ammonium nitrate/fuel oi1 (ANFO) (approximately 5.65 1b/ft? 
TNT equivalent). The front array was loaded with 3.84 1b/ft* of ANFO 
(3.45 lb/ft? TNT equivalent). The arrays were separated by 100 ft to in- 
sure the integrity of the front array since the time delay was set at 1.5 s. 


Due to a firing system malfunction, only the back array fired during 
the main test’. The front array remained intact and was fired a few 
weeks later. As a result, two complete sets of DIHEST data were obtained. 
The back array test has been designated SIMQUAKE IA (SQIA) and the front 
array test SIMQUAKE IB (SQIB). 


In Figures 5, 6, and 7, the horizontal motion data for SQIA and SQIB 
at the mid-depth of the arrays are compared with the prediction relations 
given earlier. The level of agreement is quite reasonable with the 
exception of the acceleration prediction for SQIB which jis low by about 
a factor of two. Since the prediction relations were developed indepen- 
dent of the data, the quality of agreement suggests that future DIHEST 
events in soil may be predictable within a reasonable level of confidence. 


*Two sequentially timed arrays were recently fired in SIMQUAKE IT at 
McCormick Ranch. 
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Figure 5. Comparison of SIMQUAKE I Acceleration Data and Predictions 
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Figure 6. Comparison of SIMQUAKE I Velocity Data and Predictions 
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Scaled Displacement, d/o [in/(1b/ft7)] 
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Figure 7. Comparison of SIMQUAKE I Displacement Data and Predictions 
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IV. GENERALIZED EMPIRICAL PERFORMANCE CHARTS 


In order to prescribe a DIHEST design procedure, it is desirable to 
have generalized performance charts of a relatively simple nature. The 
major complication of the empirical equations presented in the previous 
section is the dependence of the normalized transition range, Rg/H, on 
H/o. For DIHEST events with H/a in the range of 1 to about 50, average 
Ro/H values for acceleration, velocity, and displacement are 0.75, 0.82 
and 0.61. Figures 8 and 9 compare acceleration and velocity data from 
DIHEST Improvement Program Test IIA (DIP IIA) and Mini-SIMQUAKE-F (MSQ-F) 
(Ref. 7) with predictions using the equations with variable Ro/H, and 
with predictions using the constant values of Ro/H given above. DIP IIA 
and MSQ-F span over more than an order of magnitude in both areal charge 
density and scaled array height. It can be seen that constant values 
of Ro/H provide predictions which are low by about a factor of 2.5 for 
DIP IIA and high for MSQ-F. 


A similar comparison for displacements in DIP IIA and MSQ-B is given 
in Figure 10. MSQ-B is used because MSQ was a double-array event and 
displacements from MSQ-F cannot be estimated accurately because of their 
interaction with MSQ-B. The predictions for MSQ-B used the modification 
for shallow events discussed earlier (Section III). In Figure 10, the 
predictions for DIP IIA using constant Ro/H is about 60 percent low. The 
predictions for MSQ-B are about 10 percent high. 


If deep events are considered, the assumption of constant Ro/H 
results in errors ranging approximately from a factor of -2 at smal] 
H/o to a factor of +2 at large H/o. At intermediate values of H/a, 
the error will be less. For preliminary design purposes these results 
will usually be acceptable. 


Figures 11, 12 and 13 are generalized performance charts for DIHEST 
arrays in soil, based on the constant Ro/H values given previously. 
Figure 11 shows scaled accelerations (ata) versus scaled range (R/a), 
where the scaling factor, a, is the areal charge density of the DIHEST 
array in 1b/ft?. Figure 12 shows particle velocity (v) versus scaled 
range. Figure 13 shows scaled displacement d/a versus scaled range. 


Figures 11, 12 and 13 can be used to predict the approximate 
centerline, midheight ground motions from a given soil as follows: 


(a) For the particular parameter of interest, determine the 
transition-point range, Ro, by multiplying the array height, H, by 
the appropriate constant (0.75 for acceleration, 0.82 for velocity, 
and 0.61 for displacement). ; 


(b) Plat the planar attenuation zone down to the scaled transition- 
point range (Ro/a). Then plot the cylindrical attenuation zone (inter- 
polating where necessary). 


(c) The resultant plots can be “unscaled" by using the appropriate 


value of a to yield direct plots of a versus R, v versus R, and d versus 
R, 
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Figure 9. DIP IIA and MSQ-F Velocities Compared with Predictions 
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Figure 10. DIP IIA and MSQ-B Displacements Compared with Predictions 
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Figure 11. Soil DIHEST Accelerations 
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Figure 12. Soil DIHEST Velocities 
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Figure 13. Soil DIHEST Displacements 
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¥. DESIGN PROCEDURE 


The design of a DIHEST array to simulate a given weapon yield has 
not previously been a well-defined procedure. It has been limited by 
the lack of a consistent prediction procedure for determining ground 
motions resulting from nuclear explosions and by the lack of a definite 
design procedure for DIHEST. It is believed that prediction procedures 
for large single-point explosions are still lacking. However, a design 
procedure is now available for determining the DIHEST dimensions re- 
quired to produce a given set of ground motion parameters. 


In order to design a DIHEST array, it must be established that the 
required peak parameters are compatible (i.e., that all of the peaks can 
be derived through differentiation or integration from the other wave- 
forms). Generally speaking, any set of peak parameters obtained from 
data from a large explosion or an earthquake will be compatible. 


Given a compatible set of ground motion parameters, the current 
DIHEST design technique consists of the following procedures: 


(a) On Figure 12, plot the required value of velocity. Use the 
intercept of the velocity with the planar attenuation line 
as the initial guess for a value of the scaled velocity 
transition range [(Ro)y/a]. Use this value also as the ini- 
tial value of scaled target range (R,/a). 


(b) Plot the value of (Ro)y/a and Ry/a on Figures 11 and 13 and 
find the resultant value for scaled acceleration (aa) and 
scaled displacement (d/a). With the desired value of accel- 
eration, solve for a required a. Also solve for a using the 
desired value of displacement. If the two values of a do not 
agree, select a new value of (Rq),/a. If the a for accelera- 
tion is smaller than the a for displacement, decrease the 
value of (Ro) /a. 


(c) With the newly selected value of (Rg),/a, go to Figure 12 and 
find the required value of Rt/a [which now will be different 
from (Ro)y/al. Repeat step b until agreement is reached be- 
tween the values of a. 


(d) Using the values of (Ro)y/a and a found above, calculate the 
array height, H. Check the ground motions which result from 
the selected values of a and H using the equations of Section 
III. Assume a value of L (usually about 2 to 3 times H) to 
calculate the charge weight of the total array. 


It should be noted that the above design procedures are based on ground 
motions along the centerline of the array. This is adequate for deeply- 
buried target structures. However, for surface target structures a 
correction is necessary. The procedure for this correction is given 

in Reference 7 and can cause changes up to about a factor of two. In 
general, all parameters (acceleration, velocity, and displacement) 
increase as the surface is approached. Acceleration is approximately 
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the same; velocity is increased by about 50 percent; and displacement 
is increased approximately 100 percent. 


The relationships between the DIHEST array variables of H, a, 
(Ro)y/a, and total charge are shown in Figure 14. The array length-to- 
height ratio was assumed to be three. This plot shows that large values 
of (Ro)y/a rapidly lead to large array heights and very large total 
charge weights. It is therefore necessary to minimize the values of 
(Ro)y/a used in the above design procedure. Sometimes this may require 
some degree of compromise in the required ground motion parameters (such 
as reducing displacements or allowing larger accelerations). While us- 
ing Figure 14, it should be noted that H is usually limited by geology 
(bedrock or water table) or by drilling limitations, i.e., depth limita- 
tion of available drilling equipment. In-place explosive costs vary con- 
siderably, due primarily to the drilling characteristics of the geology. 
Recent experience with 70-ton DIHEST arrays in soil with drilling depths 
of 100 ft suggest in-place array costs of approximately $1000/ton of 
explosive. This approximation allows Figure 14 to be used directly as 
a DIHEST cost chart with the total explosive weight read as kilodollars. 
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Figure 14. Relationship Between DIHEST Variables 
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VI. CONCLUSIONS 


During the past few years, considerable gains have been made in 
understanding the performance of DIHEST arrays in dry soils and in the 
ability to design and predict their performance. In the analysis of 
previous data, the scaling relationships and generalized performance 
charts were developed. An empirical design procedure was also 
developed. It is therefore now possible to design and predict DIHEST 
array performance in dry alluvium soils with a reasonable level of 
confidence. 


Methods of analysis and prediction are now at a level at which 
design relationships can be developed for a new site and for different 
geologies with only a limited number of calibration shots. 


Ground motion experiments can be designed to have amplitudes and 
frequency contents comparable to the range of those expected in the 
ground-roll region of nuclear events and to those of very strong 
earthquakes. 
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ABSTRACT 


SIMULATION VS STIMULATION 


by 


Robert J. Port 
Henry F. Cooper, Jr. 


Testing plays a central role in all engineering and scientific 
enterprises, Tests often supply fuel to the creative process 
of conceiving and evolving new concepts; they refine and focus 
ideas in the engineering development and validation of designs, 
and they provide a basis for monitoring the quality of the 
products. In particular, testing has played an essential role 
in evolving U.S. land-based strategic systems, and in provid- 
ing confidence in their Survivability to nuclear attack. 

Since the atmospheric test ban conStrains, and in many cases 
prohibits, actual exposure of many system components to the 
expected nuclear environment, the U.S. has expended significant 
resources to develop and apply test procedures to subject sys- 
tem components, and entire systems, to Simulation environments. 
The fidelity of the simulation achieved in these tests has 
varied depending on test objectives, our knowledge of the 
environments being simulated, and our testing ability. 


This paper reviews the history of blast and shock testing 

for land-based systems, comments on the degree to which 
previous testing procedures has reflected our changing under- 
Standing of the nuclear environment and our ability to test, 
and discusses the requirements for high- vs low-fidelity 
testing, i.e., Simulation vs stimulation. Both design and 
targeting applications will be considered. 


BACKGROUND 


In any engineering endeavor, tests/experiments play key roles 


in: 
e understanding phenomena and conceiving systems. . 
developing basic principles. 
e evolving designs . . . identifying and correcting 
problems/improving concepts 


e proving designs in their final stage .. . proof/ 
acceptance testing 


e assuring design is implemented in production... 
quality control 


e assuring system is maintained as designed... 
quality asSurance. 


Nuclear weapon systems are no exception. 


However, nuclear weapon systems are distinctive in at least 
one sense. The final proof of the utility/success of a sys- 
tem/product is usually measured by its operational perform- 
ance following production . . . presumably under the condi- 


tions for which it was designed. In contrast, nuclear weapon 


systems have not been employed under anything approaching 
realistic wartime conditions, with the possible exception 


of the bombs dropped in Hiroshima and Nagasaki under a benign 


resisting environment. We are totally reliant on peacetime 
studies and testing for our understanding of the nuclear 
survivability and wartime reliability of our nuclear weapon 
systems, 


Everyone hopes the acid test will never occur. Indeed, the 
stated national purpose for such systems is to deter nuclear 


conflict, and in a major sense they will have failed in their 


purpose if they are ever used. Nevertheless, confidence 
that the systems will perform as designed in the event of 
a nuclear conflict is very important . . . otherwise their 
deterrent value would be significantly eroded. 


In particular, this confidence is crucial in assuring the 
existence of a second-strike posture, a keystone of U.S. 
Strategic policy. This policy requires that we have the 
capability to absorb a nuclear first strike before retaliat- 
ing with nuclear weapons. This second-strike policy places 
a premium on assuring the survivability of our strategic 
forces. If either we or our adversaries doubt that our 
forces are survivable, a potential instability could lead 


to nuclear war. If we have doubts, our adversary might per- 
ceive that we might not resist the temptation to preempt and 
minimize our losses under times of stress; hence, they could 
be motivated to preempt. To counter such unstable tendencies 
(which exist as part of Soviet doctrine and expectation of 
history), it is important that they believe our forces are 
Survivable and that preemptive attack is an unattractive 
option in any case--whether in times of stress or "“out-of~ 
the-blue". 


A key question is "what does it take to provide such confi- 
dence to ourselves and our adversaries?" Testing constitutes 
an important part of the answer to this question, and no 
Single test could be more impressive than a missile launch 
from an operational site following the detonation of a mega-~ 
ton surface burst a quarter of a mile away. For political, 
environmental, and other practical reasons, we cannot conduct 
such a test. Rather, we provide confidence in peacetime 
reliability by routinely flying missiles from Vandenberg AFB 
(under benign environments), and we have built some confidence 
in the wartime reliability of our silo systems by subjecting 
model and operational missile systems and components to test 
environments that simulate to some degree the nuclear envi- 
ronment expected under nuclear attack conditions. The inte- 
gration of results from these less than perfect tests provides 
the basis for our confidence in the hardness of our silo-based 
missile systems. 


Thus, it is important to consider the attributes of simulation 
testing and associated analysis procedures . . . including the 
key question: "How much fidelity in simulating the actual 
nuclear environment is necessary?" 


A purest desires a high-fidelity simulation as a matter of 
principle. On the other hand, a pragmatist recognizes that 
no Simulation is perfect; and furthermore, under some circum- 
stances either the conditions for a "perfect" simulation 
cannot be defined or the simulation test requirements are 
prohibitive. In addition, he notes that, if the system sur- 
vives "overtest" conditions, confidence would be assured and 
the entire question of test fidelity would be moot. 


Such an “overtest" is very desirable, and under some condi- 
tions is quite credibly possible. In other cases, the 


*It should be noted that whoever Ppreempts enjoys an advan- 
tage since peacetime reliability must be higher than wartime 
reliability . . . but it might not be good enough and this 
uncertainty adds to deterrence. U.S. survivability testing 
should have as a goal the maximizing of this uncertainty in 
the mind of our adversary. 
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uncertainties are so large that overtesting bounds provide 
very severe loading conditions, and the cost of assuring a 
sufficient hardness condition could be very high. Further- 
more, if the environment of primary concern is not known, or 
cannot be bounded with confidence, then it is not logically 
possible to define an overtest. 


For blast and shock testing, an input boundary condition is 
applied to the material within a control region, and meas- 
urements are made of both the loading and response parameters. 
A typical example is a structural response experiment where 
the loads are applied to the boundaries of the structure, and 
the prime output meaSurement is the stress and motion time 
history of the responding structure. Other examples involve 
experiments which load the air/ground interface over an earth 
half-space, loadings on soil and/or fluid islands around 
structures, loading on structures and portions of structures, 
and loadings on internal system and components. 


All the above experiments stimulate the test article in that 
they apply a loading to the boundary of a control region and 
excite a response of the material within the region. When the 
stimulation boundary condition duplicates or reproduces some 
aspect of a nuclear weapon effects input boundary condition, 
the experiment is said to simulate that aspect of the nuclear 
event. The simulation portion of the experiment must be accu- 
rately qualified to define the precise aspect of the nuclear 
environment that is being simulated, and the temporal and 
spatial extent of the simulation. 


Tests that provide high-fidelity simulation loadings are not 
necessarily a preferred means of understanding or qualifying 
the response of a test article. There are many problems 
which are better solved by developing the response function 
to a group of stimulation experiments. AS an example, the 
analytic response functions from step loading and impulsive 
loading (Green's function) can be combined to solve for the 
response of a complicated linear system to a complicated 
nuclear loading function. A similar procedure has been 
proposed to estimate the structural response from nuclear 
environments using a combination of the structural response 
measurements from a constant-pressure step input and an 
impulsive loading experiment. 


The fidelity of the test loading conditions refers to the 
exactness of the nuclear simulation and not the quality of 
the experiment. A low-fidelity stimulation test should pos- 
sess all of the qualities of predictability, repeatability, 
and measurability required in any worthwhile experimental 
research, 
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ENVIRONMENT FIDELITY ISSUES 


In attempting to conduct simulation experiments, we must know 
the environment well enough to define the simulation goals, 
and we must define and measure the test loading conditions 
well enough to meet these goals. This task is not easy. The 
first main problem is in specifying the nuclear environment 
to be simulated. As we will see in the next section, there 
are large uncertainties in defining nuclear environments, and 
many paSt programs have spent considerable resources gener- 
ating a high-fidelity simulation of someone's best estimate 
of a very uncertain nuclear environment. While this approach 
may be appropriate for some applications, we feel that an 
alternative approach of providing clear bounding overtests 
with lower fidelity experiments might provide higher con- 
fidence in test programs where demonstrating survivability 

is the objective. 


The large uncertainties in nuclear environments make it 
extremely difficult to address the actual fidelity issues 
within a given test program. When the environment is uncer- 
tain by factors of two or three, then any test that pro- 
vides an environment even close to the nominal design goals 
can realistically be considered to represent the salient 
features of a possible nuclear occurrence. 


One of the biggest problems in dealing with nuclear weapon 
effects uncertainties is differentiating between the random 
uncertainties and the systematic biases or ignorance uncer- 
tainties. The random uncertainties (such as variations in 
Structural material strengths about their mean values, or 
the azimuthal variation of weapon effects) can be evaluated 
with test environments that simulate nominal nuclear weapon 
effects. Such nominal environments can be used to predict 
the average results of a large number of weapons against a 
large number of targets, because the random uncertainties 
tend to integrate out to the nominal values. 


Systematic uncertainties, involving ignorance of either the 
nuclear weapon effects environment or the resulting system 
response analysis, cannot be addressed so easily in estab- 
lishing the fidelity requirement for a test program. Our 
present prediction techniques generally cannot resolve sys- 
tematic variations in the phenomenology of less than a fac- 
tor of two, and sometimes the possible bias errors are much 
larger. As an example of the difficulties resulting from 
these systematic uncertainties, consider the problem of 
designing a ground motion simulator to exercise the shock 
isolation system within a structure at a given site. What 
environment goal should be used? Since only one nominal 


environment is correct and we just do not know what it is, 
it seems obvious that the philosophy of developing a conser- 
vative overtest is mandatory. It also seems obvious that 
high fidelity for a conservative overtest should not be a 
mandatory requirement if achieving that fidelity is either 
difficult, time-consuming, or expensive. 


The question about fidelity is also strongly influenced by 
the ability to design and measure the test environment. 
There are some weapon effects, such as fireball tempera- 
tures, that cannot easily be simulated without rather large 
and very expensive test programs. Although the test program 
objectives may demand such tests, the available manpower, 
funding, and/or state-of-the-art simulation technology may 
not be available. It may therefore be necessary to settle 
for a low-fidelity stimulation test program and/or analysis 
effort to provide a low-confidence estimate of the required 
system response. 


The size of an experimental program and its required fidelity 
are strongly related. It may be very realistic to conduct 
high-fidelity experiments on statistically significant num- 
bers of component level tests, but prohibitively expensive 

to conduct a large number of full-scale experiments. Almost 
all large~scale tests consist of only one or two samples. 
Those who have studied the extrapolation of the results of 
experimental programs that have only a few data points 

will realize that it is not possible to make high-confidence 
Statements about system response issues from so little data. 
When the program plan is designed in such a way as to preclude 
high-confidence statements, one should seriously question any 
requirement for high-fidelity experiments. 


Any experiment that does not achieve a perfect simulation 
must address the significance of the nonsimulation response. 
There are two issues involved: those characteristics that 
are within the nuclear environment but are not part of the 
Simulation test, and those characteristics that are within 
the simulation test that are not part of the nuclear envi- 
ronment. Both issues are extremely important to the inter- 
pretation of the test results. 


It is also possible to have a high-fidelity simulation of 
one portion of the nuclear environment, such as vertical 
motion, and not attempt to simulate another portion such as 
the horizontal motion. The test could then be both a high- 
and low-fidelity simulation, and the judgment of the system 
response implications would depend on the critical failure 
mode analysis. We have had vertical simulation tests on 
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large~scale systems which turned out to be extremely vulner- 
able to the horizontal environment that was omitted from the 
test. One should now classify such a test as a low-fidelity 
system experiment. 


We also have performed tests that failed structures with 

some aspect of the test environment which did not simulate 

any nuclear weapons effects (such as debris impact from the 
HEST overburden). Such tests are difficult to interpret, and 
may be very misleading. It is therefore important to evaluate 
the fidelity requirements, and assess the impact of potential 
or planned low-fidelity aspects of the environment early 

in program formulation. 


ENVIRONMENT UNCERTAINTIES 


The first question one should ask when discussing simulation 
testing is what environment Should we be trying to simulate. 
It will become very apparent if one studies nuclear weapon 
prediction techniques in any detail that our understanding 
of weapon effects environments is very uncertain. Table l 
displays the approximate peak value uncertainties of several 
critical blast and shock environment parameters. Note that 
the blast and shock uncertainties are measured as factors 
expressing the ratio of upper to the lower bound estimates. 
Note also that the minimum entry in the table is the factor 
of two uncertainty in peak overpressure. This means we cannot 
estimate the peak overpressure at a given range from a 
nuclear detonation to better than a factor of two, and all 
other nuclear environments are more uncertain than the 
overpressure. 


The risetime to the peak overpressure is alSo very uncertain 
in some overpressure regions, and may be a strong function 
of the properties of the ground surface material. Figure l 
shows pressure-time data along three specially prepared 
radial lines from a single nuclear event: an asphalt line, 
a water line, and a desert line. Note that the asphalt and 
desert lines show a distinct precursor arrival before the 
main overpressure waveform, The water line is seen to 
approach the predictions for an ideal surface (also shown in 
Figure 1). It is assumed that the finite width and depth of 
the water line is the main reason for its less than ideal 
behavior, since a water surface should act like an ideal 
surface. 


There have been a number of airblast experiments that have 
attempted to provide high-fidelity simulation of the ideal 
waveform shown in Figure l. In particular, there are cur- 
rently such efforts in support of a system to be deployed 
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in a southwestern desert region. We maintain that to repro- 
duce the high-fidelity ideal curve actually would provide 

a low-fidelity stimulation of the expected nuclear pre- 
cursed overpressure waveforms for some overpressure levels. 
It should also be noted that the ideal waveform bounds the 
peak pressure, but need not bound the initial impulse and 
does not bound the peak dynamic pressure (not shown). One 
should be very careful when extrapolating such ideal wave- 
form tests to make survivability estimates for systems 
deployed in regions with nonideal surfaces. 


Ground motion experiments have shown that geologic materials 
are not homogeneous and isotropic and can introduce addi- 
tional uncertainties in ground motion waveform characteris- 
tics induced by the overpressure loading at the air-ground 
interface. Figure 2 shows the range of waveform variations 
observed on six azimuths at a given range from a 20-ton 
explosive source. Such azimuthal variation might also be 
expected in the ground motions from nuclear explosions. Peak 
ground motion is seen to vary by at least a factor of two for 
this single experiment. 


Note also in Figure 2 that the transverse motion is at least 
25 percent of the radial motion. A word of caution is 
directed to system designers who think the ground only moves 
radially and vertically in the vertical plane passing through 
the axis of symmetry. The amplitude of such transverse motion 
has often been observed to approach the amplitude of the ver- 
tical and radial motions. 


The vertical and horizontal waveforms shown on Figure 3 dis- 
play potential variation from different geologies. The data 
were measured at ~5-ft depth at ~100 psi from surface tangent 
high-explosive events that delivered approximately 100 tons 

of equivalent TNT energy. Unfortunately, the hard-rock experi- 
ments did not have data at this range, and the variation shown 
in Figure 3 is restricted to soil or soft rock geologies. 


The first arrivals in the top three materials (waveforms) of 
Figure 3 show an initial downward or outward airblast-induced 
ground shock, with the largest amplitudes in the highly com- 
pressible dry playa material. The bottom materials have near- 
surface water layers that cause the first motion to be an 
upward and outward outrunning ground shock. The direction 

of first motion is therefore dependent on the geology and 
range from the detonation. 


Figure 3 also shows considerable variation in the number of 


cycles and the frequency content of the overall waveforms. 
It is interesting to note that motion has not stopped after 
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one second from the 100-ton events, and less than 2-Hz fre- 
quencies are observed. If one is designing a system to be 
deployed in a range of geologies, it is necessary that the 
system performance be insensitive to the many and random 
details of the waveform characteristics. 


An example question central to the subject of this paper is 
how should one design a ground motion simulator to represent 
the input to a shock isolation system that will be located 
at the 100-psi range from even a 100-ton high-explosive 
event, much less a 1-MT nuclear event. After studying 
Figure 3, one should quickly come to the conclusion that 

a single high-fidelity stimulation is not meaningful. The 
best approach would involve test procedures that vary the 
Magnitude and direction of the initial motion, as well as 
the number of cycles and frequency content of the remaining 
wave trains. 


SIMULATION HISTORY 


The history of DNA simulation testing goes back to the Senate 
ratification of the Limited Test Ban Treaty in 1963. The 
chart on Figure 4 shows over 60 major experiments that have 
occurred over the past 14 years in support of several strate- 
gic weapon systems. We believe the chart includes all expe- 
riments that involved more than 20 tons of high-explosive and/ 
or tested structures of larger than 1/6 scale. The overpres- 
sure ranged from a few to several thousand psi, and the 
simulated yields ranged from 20 tons to a few megatons. 

The overall cost has not been compiled, but it is of the 
order of a few 100 million dollars, 


It is interesting to note there is almost no duplication. 
Only one or two experiments are repeats of previous experi- 
ments. In all other cases, sequential tests in the overall 
program have made significant changes in the yield, overpres- 
sure, scale size, geology, or structure being tested. In 
other words, there has not been a statistically significant 
Sample size in the entire program. The forces that mold 
Program plans tend to lead to only a few large-scale system 
tests for satisfying the stated program objectives. One 
should seriously think about the confidence level that can 
be achieved with only a few data points, before making deci- 
sions on developing high-fidelity simulators for proposed 
future test programs. 


The structural systems involved in these test programs 
include Minuteman, shelter and trench MX, deep underground 
basing, generic targets, bridges, POL, vehicles, tanks, and 
military equipment. Although several experiments raised 
more questions than they answered, they usually satisfied 
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the objectives or the specific system questions being 
addressed by the test program, 


The question of fidelity in past simulations has to be 
answered as a function of time, . Many previously purported 
high-fidelity experiments are now considered low-fidelity. 
In several instances, our understanding of the environment 
has changed, which would change our original interpretation 
of the fidelity of the experiment. We have also learned 
things about the system being tested which would indicate 
other environments were more critical than the design high- 
fidelity environment. In these cases, the experiment may 
have been a high-fidelity simulation of some portion of the 
environment but a low-~fidelity system response experiment. 


These experiments are unique test events that will probably 
never be repeated. The data, in some cases, provide the 
only useful means to check and extrapolate our prediction 
techniques to make estimates about the response of proposed 
systems. This extrapolation is usually difficult and uncer- 
tain because the test events usually did not include the 
correct structure and/or environment for direct application 
to future problems. This extrapolation does however repre- 
sent one of the significant uses of the results of an experi- 
mental program, and may in some cases be more useful than 
the original experimental objective. Several experiments 
have uncovered phenomenology gaps that have completely 
changed the direction of the major research programs, 


It is interesting to review the initial fidelity objectives 
of the past test programs. In general, the various simula- 
tion programs started off as low-fidelity experiments that 
gradually changed to produce better and more accurate simu- 
lations as the program continued. As an example, HEST 
testing was originally developed and applied because of 

a concern that low-level vulnerabilities might exist in 

a multi-silo attack scenario. The research community first 
attempted to provide a reasonably good simulation of the 
airblast loading. The emphasis on the overall fidelity of 
the HEST simulation came later. 


The DIHEST development was designed to simulate a best esti- 
mate for the direct-induced ground motion input in a hard 
rock geology. It was fully recognized that the state of the 
art for predicting ground motion was very poor, but it was 
not as bad as the state of the art for predicting the system 
response to typical values for the horizontal input. The 
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DIHEST is therefore a low-fidelity simulation, in the termi- 
nology of this paper, that can be used to evaluate/develop 
system models for predicting direct-induced ground shock 
response. 


The MINE THROW series was designed to provide a high-fidelity 
Simulation of the crater-related ground motions. The pro- 
gram was supposed to evaluate the simulation capability by 
comparing first with JOHNNY BOY, and then the CACTUS nuclear 
events. Unfortunately, the program was stopped for political 
reasons before we were able to perform the CACTUS simulation 
at Eniwetok. As will be shown on a later chart, MINE THROW I 
produced a very good fidelity simulation of the low-frequency 
Surface motion. : 


The BLEST was designed to provide a low-fidelity stimulation 
of the portion of the upstream airblast that would maintain 
the stress loading on the underground portions of a struc- 
ture for the required simulation time. The BLEST is not and 
never was a high-fidelity simulation of the upstream air- 
blast-induced ground motion. 


The DABS was an attempt to develop an inexpensive low-fidelity 
Simulator of high-pressure dynamic airblast loading. Alter- 
native techniques that produce higher fidelity simulation 

were much more expensive, but did eliminate the dust and 

high explosive product environment loading on the structure. 
The DABS experiments actually produced much better fidelity 
than was originally expected. 


The in-trench testing was clearly a low-fidelity test program 
Since it was not known, within rather large uncertainty bands, 
what to expect for the in-trench environment. 


The deep underground nuclear test events were phenomenolog- 
ical experiments that were designed to measure structural 
response and structural fragility under direct-induced 
ground shock loading. These experiments are obviously high- 
fidelity simulations of underground nuclear events, but of 
questionable fidelity for interpreting structural response 
to near-surface nuclear detonations, 


The high-explosive events involving stacks or spheres of high 
explosives do a good job of simulating the low-overpressure 
regions of a nuclear surface event. Comparison between the 
HE and NE will be shown in the next section. 


In summary, there are some high-fidelity aspects of all the 


past simulation experiments, but each test also has low- 
fidelity portions of the environments. The best simulations 


192 


have been performed on small-scale structures or at very low 
overpressures. All the high-overpressure experiments on full- 
scale structures have been, in our opinion, low-fidelity 
simulations, 


SIMULATION COMPARISONS 


The scope of this paper does not permit an effective review 
of the capabilities of all our simulation experiments, but 
we can comment briefly on some of the more important and 
most used techniques. 


Stacks or spheres of high explosives have been used exten- 
sively in the study of ground shock phenomenology, and the 
study of structural response in the low overpressure regions. 
Figure 5 compares free-air airblast data from high-explosive 
events with the peak pressure from an equivalent point source 
of energy. A nuclear event would produce results similar to 
the point source curve. It can be seen that high explosives 
grossly underestimate the peak pressure and impulse in the high 
overpressure region (greater than a few tens of bars). This 
low-fidelity simulation in the high-pressure region makes it 
very difficult, if not impossible, to extrapolate the upstream 
airblast, direct-induced ground shock, crater-related ground 
shock, and late-time ground roll from near-surface, high-explo- 
Sive events, Although the HE/NE extrapolation is questionable, 
the high-explosive tests are the only source of information in 
some geologies, and therefore are the only data to guide our 
predictions for nuclear weapons effects in these geologies. 


At low overpressures, less than ten bars, the high explosive 
is seen to do a very good job of simulating both the peak 
pressure and waveforms. The HE to NE extrapolation for low- 
pressure, air-induced ground shock can be performed with much 
higher confidence. The response of structures can also be 
evaluated, to the low-overpressure waveforms, and thus pro- 
vide confident predictions in the low-overpressure region. 


In most cases, large-scale high-overpressure simulations have 
been peformed with some version of the HEST technique. Fig- 
ure 6 shows six gauges, all located at the same range in the 
Hardpan I-3 HEST cavity. The heavy solid lines on the figure 
represent the desired nuclear overpressure waveform at the 
1200-psi range from a 125-kT detonation at a 210-ft HOB. 


Figure 6 shows that a HEST test only simulates the nuclear 
overpressure in a gross sense, since there are several 
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obvious disparities in the high-frequency region. The actual 
peak measured stresses occurred as very high frequency spikes 
with amplitudes that ranged from 1800 to 3600 psi. Large 
structural response experiments and ground motion phenomeno- 
logy experiments generally are not sensitive to these high- 
frequency spikes, and therefore can be adequately tested 

with a HEST environment. The HEST waveform also contains 
several oscillatory signals between 500 and 5000 Hz that 

are not represented in the ideal nuclear waveform, Caution 
must be exercised when using HEST techniques on model struc- 
tures where the characteristic frequency of higher order 
failure modes may approach these nonsimulation oscillatory 
levels. In ground motion experiments, the upper few feet 

of soil usually filters out the high-frequency oscillations, 
permitting a much better simulation of the predicted ideal 
nuclear waveform, 


Figure 7 shows a comparison between the ground motion pro- 
duced by the cratering effects of a nuclear event and the 
MINE THROW simulation technique. The JOHNNY BOY experi- 
ment was a 0.5-kT nuclear detonation at a 1.92-ft depth of 
burial. The MINE THROW I event was a buried stack of ANFO 
explosive that resembled a thick bowl configuration. The 
MINE THROW concept does not attempt to simulate the airblast, 
and therefore does not provide the initial airblast-induced 
Motion in the JOHNNY BOY waveform. The direct-induced, 
crater-related, and late-time ground roll seem to be very 
similar to those of the nuclear event. 


NEED FOR HIGH FIDELITY 


Having reviewed some of the past history and compared the 
nuclear events with some of the simulated waveforms, the 
question remains aS to what is the specific need for high- 
fidelity simulation experiments. Although the final answer 
to the question can only be evaluated with respect to the 
specific application and structural response problems under 
study, the following collection of observations are general 
to many applications. We will consider the applications in 
two very different categories: offense and defense. 


In offense applications, we have complete knowledge of our 
weapon systems (yield, HOB, etc.) but have incomplete knowl- 
edge of the target systems. Although we may know the gross 
details of structures, we seldom know specific details such 
as the surrounding geology, percentage of reinforcing steel, 
Strengths of structural materials, or amount of shock isola- 
tion on critical equipment. The target uncertainties there- 
fore dominate the target response problems. The targeteer 
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is therefore forced to use a conservative overkill environ- 
ment if he is to achieve high-confidence target kill. 
Consequently, the required test program objectives are to 
provide environments that will define the sure-kill levels 
of poorly defined structures. It is hard to make a case 

for a high-fidelity environment, when the objective is to 
develop a sure-kill level on an uncertain target. It would 
seem in general that a low-fidelity environment would be ade- 
quate for providing a least upper bound kill level for a 
generic class of targets. A high-fidelity environment may 
be useful to explore potential low-level vulnerabilities for 
very important systems in which the high-confidence, sure- 
kill level is larger than the available weapon arsenal can 
afford, 


Defense applications are almost opposite to the offense case. 
We have complete knowledge of the geology and structural 
configuration, but have little knowledge about the threat 

and resulting nuclear environments. The threat and environ- 
ment uncertainties therefore dominate the system response 
analysis and survivability calculations. The designer must 
therefore use a conservative overdesign approach to provide 

a high-confidence survivable system against an uncertain 
threat environment. It is hard to make a case for a high- 
fidelity test of an uncertain environment to provide confi- 
dence in a successful overdesign approach. It is possible 

in these situations to perform cost tradeoffs between 
confidence, fidelity, and system overdesign; but one has to 
be careful to adequately treat the environment uncertainties, 
and not perform detailed cost tradeoffs based on best-esti- 
mate design goals. 


Although it may be possible to design adequate low-fidelity 
overtests for those environments and response modes we know 
about, the major concern is for unknown vulnerabilities and 


response modes. There is a need for a high-fidelity simula- 
tion to check the final defensive system to insure that there 
are no unknown vulnerabilities to the high-fidelity aspects 


of the environment. The best case for a high-fidelity test 
of an uncertain environment is to provide confidence in a 
successful’ overdesign approach. 


SUMMARY ES Sh 3 


Since atmospheric nuclear testing has been outlawed, one must 
rely on simulation environments to provide the critical test 
levels for system development program objectives. Simulation 
testing is therefore a fallback position, and the question is 
"how much fidelity in simulating the actual nuclear environment 
is necessary?" 
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The question of required simulation fidelity is not an easy 
one to answer. Just to say simulation implies that we know 
the environment well enough to define the simulation goals, 
and believe that we can measure and provide an environment 
which meets these goals. There are, however, large uncer- 
tainties in predicting the nuclear environments to be simu- 
lated. These large uncertainties may imply that almost any 
environment close to the design goals may be realistically 
considered as a possible nuclear occurrence. It seems clear 
that to attempt to provide a high-fidelity approximation to 
someone's best estimate of a very uncertain environment is 
not a sound program objective, It may be more desirable to 
bound the environment and provide a clear overtest to insure 
a system survivability (or system vulnerability for offense 
application). Such a test may be of lower fidelity but of 
higher confidence. ° 


High-fidelity simulations are generally much more costly and 
require several calibration events to insure that the full 
system tests will meet the design test objectives. Low-fide- 
lity overtests are generally less expensive and require less 
proof of capability prior to full-scale testing. One cannot 
neglect questions about scaling and nonsimulation response 

in selecting the more desirable fidelity for implementing a 
test program. Nor is it possible to ignore the ability to 
design and implement a test program before establishing the 
fidelity objectives. Whether for high- or low-fidelity objec- 
tives, a test environment must be predictable, reproducible, 
and measurable, consistent with good experimental practices. 


The need for high-fidelity experiments depends strongly on 
how well we know the environment and the system response to 
the environment. When we know both the environment and 
response, then both high- and low-fidelity experiments can 
be designed and suitable cost tradeoff studies performed. 
When we know the environment but do not understand the sys- 
tem response, it is mandatory that we have a high-fidelity 
test environment. Otherwise, we would not be able to evalu- 
ate the deficiencies in a low-fidelity test program. 


The most common case is when the environment is not known accu- 
rately, and we are trying to develop an understanding of the 
system response. The less we understand of the system, the 
more we desire a fidelity simulation; but the less we under- 
stand of the environment, the more difficult it is to assure 

a high-fidelity simulation. In many cases, it would be more 
desirable to channel resources into performing several low- 
fidelity experiments that emphasize different aspects of the 
environment, than it would to perform only one or two higher 
fidelity simulations of very uncertain environments. 
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In the past, DNA has funded several large-scale simulation 
tests for application on Minuteman, Hard Rock Silo, Minuteman 
Upgrade, several MX concepts, Submarines, vehicles, tanks, 
and various military equipment. In many cases, the tests 
were advertised as high-fidelity simulations when in fact we 
learned more about the environment than we did about the sys- 
tem response. Many of the environment prediction techniques 
have changed over the years since the tests, which means the 
interpretation of the test results has changed from a high- 
fidelity to a low-fidelity application. Our understanding 

of the critical system response has also changed over the 
years, and some of the high-fidelity environments were not 
necessarily addressing the weakest links in the system design. 


The data from large-scale system tests provide unique informa- 
tion that probably will never be repeated, and thus form the 
primary data base for estimating the response for upgraded as 
well as new and proposed system response. Since the test data 
are being extrapolated to regions beyond those considered in 
designing the original test program, the extrapolation is nor- 
mally considered of low fidelity. Such low-fidelity extrap- 
olations have, however, been one of the most useful results 

of many of our past large-scale test programs, 


It seems clear that the pretest planning rationale has to 
seriously consider many issues before establishing the best 
fidelity for proposed test environments. Long-range consid- 
erations as well as immediate test objectives must be weighed 
before optimum decisions can be made. Stimulation versus 
simulation should therefore be a central issue of future DNA 
test program planning. 
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HE Blast Pressure Measurements 
by 


Capt David J. Ray 
Air Force Weapons Laboratory 
Kirtland AFB NM 


INTRODUCTION 


The Air Force Weapons Laboratory, under sponsorship of DNA 
and SAMSO, through its in-house resources and its contract with 
the University of New Mexico Civil Engineering Research Facility, 
is developing methods of simulating nuclear blast and shock effects 
with high explosives. Two of the methods used are the High Explo- 
sive Simulation Test (HEST) and the Dynamic Airblast Simulator (DABS) 
test. 


The measurement of transient pressure pulses in the range of 
1 MPa to 138 MPa is typically required for these high explosive 
Simulations carried on at AFWL and CERF. 


THE NATURE OF THE PRESSURE TIME HISTORY 


The transient pressure pulse in a nominal 5 MPa HEST resembles 
Figure 1, based on AFWL/CERF experience. This time history should 
be closely observed for several features. Figure 1 is an 80 ms 
plot. The initial peak appears to reach 15 MPa, but in actuality, 
was greater as the channel bandedge was 15 MPa. Therefore, 
the initial peak was three times or more the nominal 5 MPa effec~ 
tive cavity pressure. The rise time of the initial peak appears 
to be about 200 Usec, but the pressure rise was probably faster 
Since the transducer was shielded from the direct blast effects. 
This rise time corresponds to a bandwidth of 1750 Hz. 


Figure 2 is the 200 msec plot. The late time history of 
the pressure pulse must satisfy the boundary condition that the 
pressure returns to atmospheric pressure. If the correct final boundary 
conditions are met, the impulse should attain a final constant value. 
The pressure time history of Figure 2 appears to fail this test for 
validity, at least in this first 200 msec. 


A known zero initial boundary conditions may also be used to 
evaluate transducer performance. Figure 3 gives a pressure time 
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history which fails to meet the initial condition due to 
light sensitivity of the transducer. This is a measurement 
from the HAVE HOST S-1 Event which was a large DABS facility. 


THE PRESSURE-SENSING SCHEME 


The transduction scheme typically used by AFWL/CERF for 
pressure Measurements in HE events is the piezoresistive bridge 
diffused onto an edge-clamped disc of silicon. The pressure 
sensors are in steel mounts imbedded in concrete cylinders 
located in the floor of the HEST cavity, and the sensing surface 
is flush with the cavity floor. 


For the DABS events the transducers are imbedded in concrete 
on surface-flush mounts or on concrete structures and sometimes 
on blast wings as in Figure 4. The gages on the wings may be 
mounted on the sides of the probes for incident overpressure 
or in the probe tips for stagnation pressure. 


Suitable piezoresistive transducers are produced by Kulite 
Semiconductor Products, Inc., and ENDEVCO. Evaluations of the 
ENDEVCO 8511 series by AFWL/CERF are now scheduled with the 
objective of determining if it will perform as well as or better 
than the Kulite in a HEST environment. Figure 5 is a section 
view of a typical pressure transducer as is available from 
Kulite. 


SOME PROBLEMS WITH PRESSURE TRANSDUCERS IN BLAST ENVIRONMENTS 


One problem encountered in edge-clamped silicon discs is 
the resonance of the brittle disc which may crack if excited. 
Further, the resistance of the semiconductor bridge elements is 
quite temperature sensitive and, as mentioned earlier, light 
sensitive. Both the Kulite and the ENDEVCO transducers use a 
layer of GE TBS 758 silastic to protect the silicon diaphragm 
from temperature changes. 


In Figure 3 the rounding of the trace just before shock 
arrival was determined to be the sensitivity of the piezoresistive 
transducer to the light radiation from the burning explosive-air 
interface. Interestingly, the transducer silicon sensing disc 
was covered with about 1 millimeter of silastic and still indicated 
sensitivity to the light from the burning interface. This led to 
the introduction of a line of pressure transducers by manu- 
facturers which included a "photon barrier" on the sensing 
disc. This barrier is very thin and adds no significant mass 
to the disc but does shield the sensing elements from the 
offending light energy. 
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Experience at AFWL and CERF has shown that the further 
addition of materials such as RTV's silastics, and aluminum foil 
in front of the pressure-sensing surface may alter the response of 
such pressure transducers and may even cause failures at lower 
than rated pressure. Figure 6 is a shock tube test of a Kulite 
transducer with an additional layer of RTV (1 mm) and aluminum 
foil over that applied by Kulite. The overshoot is visible at 
about three times the pressure level seen after the gage 
oscillations have damped out. Figure 7 shows the overshoot with 
two layers of 1 mm RTV and a piece of aluminum foil. Again, 
the overshoot is about three times the late time pressure. Also, 
the resonant frequency has been lowered. AFWL concluded that 
additional silastics used to prevent drift due to temperature or 
light effects may cause gage failure when used in a HEST or DABS 
environment. 


The construction materials used in HEST pressure cavities 
are varied and have, in the past, included metal rings, plywood, 
nails, screws and styrofoam. At detonation, these materials are 
propelled at high speed throughout the cavity. The potential damage 
to instrumentation from flying metal and wood is obvious, but 
burning pieces of styrofoam can also cause instrumentation problems. 
In order to protect the pressure gages, two shield schemes are used 
by AFWL and CERF. The blast diffusion plate (carport) and the debris 
shield are the two arrangements typically used. Figure 8 is the 
AFWL "carport." It is a 3/8-inch steel plate with a hole pattern 
drilled into it and mounted over the pressure transducer on 1-1/4 
inch spacers. This plate also has the effect of slowing the rise 
time of the pressure spikes seen by the gage. The exact amount of 
rise time slow-down is not known, but test experience seems to 
indicate a slow-down of about 20 usec. This would have the 
effect of lowering the initial peaks and enhancing gage survival. 


Figure 9 is a side view of the AFWL debris shield. This 
shield was originally developed for head-on pressure measure- 
ments in the debris-laden DABS blast flow environment. It isa 
small cap with holes, as shown in Figure 10. It threads directly 
over the pressure gage with a 0.039-inch gap between the center of 
the shield and the surface of the pressure gage. Figure 1l is a 
photo of a shock tube test on a Kulite pressure transducer with 
the AFWL Mod VIII debris shield. The upper trace has the debris 
shield and shows a rise time slow-down of about 10 usec. 
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An extensive shock tube investigation of debris shields was 
done by CERF at AFWL's request. This investigation used tested 
a cap with eight holes (Figure 12). During the test series, the 
eight holes were progressively enlarged and also tested with 
various spacing between the cap and the sensor face. Reference 1 
gives the complete test results and CERF evaluation. The resultant 
eight-hole shields have been given numbers which correlate to 
their test number in the original investigation; therefore, eight- 
hole debris shield #10 is the configuration used in Test #10 of 
the attached original CERF report. Figures 13 and 14 give the 
test results for two of the more popular configurations; they are 
the #10 and #18 eight-hole debris shields. The #10 shield has a 
rise time of about 20 Usec, while the #18 rise time is about 150 
usec. 


GAGE PLACEMENT IN THE TEST BED 


The AFWL and CERF placement of pressure gages in HEST test 
beds has been guided in the past on the intuitive feeling that, 
by spreading the explosive cords at the positions of the pressure 
gages, maximum distance from the gage to the cords could be 
obtained, thereby increasing the survivability of the transducers. 
This is now felt to be in considerable error in that CERF calcula- 
tions of pressure levels near multiple line explosive sources 
show an enhancement of the pressure midway between two such line 
sources of pressure. According to these CERF calculations, the 
position of minimum initial pressure is directly under the explo- 
sive cord, as shown in Figure 15. 


Some test evidence of pressure enhancement has been seen on 
open detonating cord tests. Figure 16 shows the result of such a 
test and indicates a region of high pressure between the explo- 
sive cord as evidenced by the plate indentation. 


CONCLUSION 


A review of these experiences indicates that, for HE 
simulations, instrumentation survival cannot be left to chance 
or intuition. The fundamental problem in ensuring survivability 
is satisfying the sometimes mutually-exclusive conditions of 
choosing maximum gage range to ensure survivability during high 
initial spikes and minimum gage range to obtain the best signal- 
to-noise ratio in the data amplitudes for the resolution required 
for the particular test. ‘ 


The following rules guide the AFWL in pressure measurements 
in HEST and DABS events and are current AFWL practice. 
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1. No silastics, RTV’s, or other materials should be 
placed in contact with transducer face. The face should be left 
as it comes from the manufacturer with only the manufacturer's 
material covering the silicon sensing disc. 


2. The initial detonation pressure peaks in the HEST 
depend on the location with respect to the explosive cord, but can 
be counted on to be greater than three times the effective cavity 
pressure. Therefore, gage ranges must be selected accordingly. 


3. The Kulite HKS-375 series transducers have a history 
of suriving in HE events. The ENDEVCO 8511 transducers are 
new and after suitable evaluation may also be included in HE 
events in order to gather data as to their survivability and 
performance in such environments. 


4. The AFWL carport should be used whenever possible to 
diffuse early-time spikes and to protect from flying debris in 
the HEST cavity. Carports should not be used in DABS events. 


5. The APWL Mod VIII or CERF eight-hole debris shield 


should be used to slow the rise time and further protect from flying 
debris. 


6. The most successful combination for HEST pressure mea- 
surements in the AFWL's experience has been both the carport and 
the debris shield. 


7. ‘The pressure gages should be placed directly under the 
explosive cord and should not be placed symmetrically between 
two cords or at a constant radius from three or more equally 
spaced cords. As much distance as possible should separate the 
cord from the gage. 
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MACH 130 AIR SHOCK 
ATTENUATION STUDIES* 


H. David Glenn, Lawrence Livermore Laboratory 
HOWARD R. KRATZ, Systems, Science & Software 
DOUGLAS D. KEOUGH, SRI International 
ROBERT P. SWIFT, Lawrence Livermore Laboratory 


ABSTRACT 

YVoitenko compressors were used to generate ~ 4.5 cm/usec air shocks 
in a steel and a grout outlet pipe. Diaphragm burst times and time-of- 
arrival data of shock-front luminosity along the 2-cm-i.d. exit pipes are 
provided by fiber-optic ports. Pressure profiles were obtained for this 
high enthalpy shock propagation for the first time and at many locations 
in both experiments. The Voitenko-grout experiment represents the first 
laboratory attempt to study shock propagation with this type of compressible 
wall material. 

The primary purpose of these two experiments was to examine the 
effect of wall material on high-energy shock propagation. In the inter- 
val between 10- and 250-cm from the diaphragm the velocity and peak pressure 
of the shock front attenuated from = 4.5 to 0.5 cm/usec and * 21 to = 0.2 kb, 
respectively. Over this distance the shock propagation gave nearly identical 
results for both experiments within the range of experimental accuracy. 
This rapid attenuation of the shock front velocity and pressure is attri- 2. 
buted to ablation and entrainment of wall material. An interesting feature 
that was observed, confirmed by multiple measurement techniques, was 
rapid oscillations of the pressure profiles. The results indicate that 
the shock propagation is independent of wall composition to the extent 
of materials considered over the 250-cm distance of the outlet pipes. 


*Work performed under the auspices of the U.S. Department of Energy under 
contract No. W-7405-ENG-48; and DNA Subtask J24AAXIX955. 
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A primary threat to the MX vehicle in the buried trench basing mode 
is air shock propagation in the shallowly buried tunnel (Figure 1). Evalua- 
tion of this threat is based on numerical codes that simulate high-energy 
air shock propagation in open pipes. 


For more than a decade experimental/theoretical studies have 
indicated that ablation and entrainment of wall material was the principal 
attenuation mechanism for high-velocity (> 1 cm/usec) air shock propaga- 
tion. '-9 The above conclusion was primarily based on agreement between 
calculational predictions and experimental results for time-of-arrival 
(TOA) of the shock front at specific locations along open pipes containing 
air at ambient conditions. Improvements have been made 5 in the early! 
theoretical description of the ablation processes contained in these 
hydrodynamic codes. However, an expanded experimental data base is defi- 
nitely needed to evaluate and/or improve the dynamic description of the 
ablation processes behind the shock front. For that purpose, a modified 
version?*° of the Voitenko compressor’ was used to generate » 4.5 cm/usec 
air shocks in a steel and a grout outlet pipe. 


Figure 2 and Table 1 illustrate the experimental setup and diagnosti 
coverage used in the steel shock tube experiment. Details of the Voitenko 
compressor operation are given in reference 8 of this symposium proceedings.* 
Following diaphragm break the air shock propagates down an 0.02 mi.d. bore 


*The diaphragm broke at a time of 58.1 usec following detonation of the 
high-explosive. 
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for a distance of 6 m. Fiber optics located in the steel wall transmit 
Tuminosity associated with the air shock to two display boards that are 
scanned by streaking cameras. A framing camera focused on the first 
0.20 m of outlet pipe indicated that no venting of gases occurred behind 
the shock front for > 40 usec after the diaphragm broke. High explosive 
gases from the driver section obscured the field of view after this time. 
Pressure gages, at locations given in Table 1, provided pressure profiles 
and TOA information for the shock front. A summary of air-shock front 
TOA's for optical and electronic diagnostics, with respect to diaphragm 
burst time, is given in Table 2. These preliminary results indicate 

a close correlation for the various diagnostics data. Some of the pres- 
sure results will be presented below when comparisons with corresponding 
grout experimental results are shown. 


Figure 3 and Table 3 show the experimental setup and diagnostic 
coverage used in the grout shock tube experiment. The Voitenko compressor 
and first 15-cm of the outlet pipe are identical with the steel experi- 
ment. Thus, experimental results for both experiments can be compared 
for similarity of flow conditions over the initial segments of the outlet 
pipes. The diaphragm broke at a time of 57.6 usec following detonation 
of the high-explosive. This is a 0.5 usec earlier than in the steel 
experiment and may be attributed to the slightly higher chamber pressure 
of 11 atmospheres compared to 10.66 atmospheres for the steel experiment. 
Following diaphragm break the air shock propagates down an 0.02 m i.d. bore 
for a distance of 3m. The walls of the bore for the final = 2.85 m 
are composed of DR-1 grout.” A summary of air-shock front TOA's for 
optical and electronic diagnostics, with respect to diaphragm burst time, 
is given in Table 4. Comparison indicates a very close correlation of 
the various diagnostics data. 


A comparison of shock front TOA for both experiments are presented 
in Figure 4. These results indicate surprisingly similar propagation charac- 
teristics over the first 2.0 m with slightly greater attenuation seen 
between 2.0 and 2.5 m for the grout experiment. 
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A further confirmation of propagation similarities can be obtained 
by comparison of a few corresponding pressure profiles from bar gages. 
Figure 5 shows the pressure profiles at 0.10 m from the diaphragm.* Not 
only are the pressure peaks in relatively good agreement but the general 
shape appears reproducible for « 20 usec behind the shock front. Framing 
camera coverage of this area indicated a small amount of venting around 
the bar gage occurred at * 25 usec after diaphragm burst and continued 
until high explosive gases obscured the field of view approximately 20 usec 
later. This may explain the reduction in the pressure profile at later 
times for the grout experiment. However, at the time venting is first 
detected the shock front has reached = 0.80 m and consequently venting 
can expect to have little effect on conditions in and for some distance 
behind the front. The main purpose of this comparison is to illustrate 
the very similar initial conditions of the flow for both experiments. 


Figures 6 and 7 show the pressure profiles obtained at 1.0 and 2.0-m 
from the diaphragm, respectively. Again we note the close agreement in 
peak values and general wave shape. A common characteristic of all the 
profiles shown is the large oscillations in pressure that appears to occur. 
The flat pack (SRI) gages in the grout experiments showed similar results 
for the first two or three oscillations both in peak magnitudes and spacing 
between the peaks. A piezoelectric (LLL) gage at 2 m in the steel experi- 
ment agreed generally with the bar gages in peak values and relative 
spacing between the peaks. However, both the flat packs and piezoelectric 
gage indicate a much smoother profile behind the first two or three 
oscillations than observed with the bar gages. Thus, although the first 
few oscillations are felt to be realistic, the subsequent oscillations 
are felt to be a problem related to the design or installation of the bar 
gages. At this time only peak values for the flat pack gages have been 
analyzed and thus no corresponding pressure profiles are available for 
comparison. 


*The time axis for the following pressure profile comparisons refer to 
detonation of the high explosive and thus have not been corrected for 
diaphragm times. 
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Figure 8 illustrates the decay of the peak pressure in the shock 
front as a function of axial distance from the diaphragm. Although the 
peak pressures in the steel experiment are consistently lower than the 
grout experiment, except for the results at 0.10 m, it is not known whether 
these differences are really significant at this time. 


Examination of the TOA and pressure profiles indicates a slightly 
larger attenuation for the air shock at 2.0 m and beyond for propagation 
in grout pipe compared to the steel pipe. To the extent of the experi- 
mental comparison, the rather surprising conclusion is that air shock 
propagation is rather insensitive to whether the wall material is com- 
posed of incompressible (steel) or relatively compressible (grout) material. 


NOTICE 
“This report was prepared as an account of work 
Reference to a company or product sponsored by the United States Government. 


name does not imply approval or Neither the cited States nor ie Oniled Stats 
r Department of Energy, nor any of their employees, 
recommendation of the product by nor any of their contractors, subcontractors, or 
the University of California or the their employees, makes any bdatecle express or 
implied, or assumes any legal liability or respon- 

US. Department of Energy to the sibility for the accuracy, completeness or 
exclusion of others that may be usefulness of any information, apparatus, product 
suitable or process disclosed, or represents that its use 

i would not infringe privately-owned rights.” 
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AXIAL 
DISTANCE 
(cM) 


0 


150 


VOTTENKO 


AIR SHOCK TIME-OF-ARRIVAL 


FIBER OPT 
TOA 
(us) 
(58.1) 0.0 
1.2 
2a 
555 
5,9 
8.5 
11,6 
14,4 
22.6 
S17 
43,3 
57.8 
F755 
103.8 
183, 
295, 
448, 
611. 


TABLE 2 


STEEL EXPERIMENTS 


Ics 53 BAR 
GAGE TOA 


(us) 


0.0 


SS a 


102.6 
(223) 187.9 
(495) 
(595,-/00,) 
(737.) 
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TABLE 3 


GROUT VOITENKO EXPERIMENTAL DIAGNOSTICS 


DISTANCE 
(cm) _____CdDIAGNOSTICS 
0 F (Fiper Optics) 
2 F 
5 F 
10 F, B, (Bar Gace) 
20 Py Bee got 
30 F, B, Py W 
40 F 
50 Pye Box anu i 
15 F 
100 Ey Bs rah 
125 F 
150 F, B, P, W 
1/5 F 
200 Pp sBo 1 Pa 
250 Fy Bo oe 
300 B 
*P = PRESSURE SENSOR (PIEZORESISTANT) 
**W = WALL MOTION SENSOR (PARTICLE 


VELOCITY LOOP) 
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TABLE 4 


VOITENKO GROUT EXPERIMENT 
AIR SHOCK TIME-OF-ARRIVAL 


AXIAL FIBER OPTICS S? BAR SRI Fiat SRI 
DISTANCE TOA Gace TOA Pack TOA = VeLocity TOA 
(cm) (us) (us) (us) (is) 
0 (57,6) 0.0 
2 13 
5 fad 
10 3,8 3,9 
20 ows 6.6 6,8 55 
30 8.6 9.6 9,4 9,4 
40 11.4 
50 14.6 15.9 15.0 13,4 
75 22.9 
100 31.8 55.1 32.4 32.6 
125 43,0 
150 57.4 58,5 58.6 56,8 
175 76,9 
200 107.7 IZ 108,4 107.4 


250 208, (218) 211.9 193.4 
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Figure 4 Air shock propagation in a steel and a grout pipe. 
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Figure 5 Pressure profiles obtained with bar gages at a distance of 0.10 m from the diaphragms 
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Figure 6 Pressure profiles obtained with bar gages at a distance of 1.00 m from the diaphragms. 
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Figure 7 Pressure profiles obtained with bar gages at a distance of 2.00 m from the diaphragms. 
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SHOCK FRONT PRESSURE VS RANGE 
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Figure &. Peak shock front pressure versus axial distance from the diaphraagm. 
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HE SIMULATION OF GROUND SHOCK-INDUCED MX TRENCH 
COLLAPSE 


J. C. Baker, H. R. Kratz, F. I. Peterson, 
K.. D.: Pyatt, J%.9.Gy Ee Bailey, 
and J. L. Waddell 


The work sponsored by the Defense Nuclear Agency 
ABSTRACT 


In order to establish the feasibility of HE simulation 
of ground shock-induced MX trench collapse, a calculation of 
the close-in ground shock from a 1 MT surface burst over de- 
sert alluvium was performed. In this way the pressure wave- 
form arriving at a full-size trench buried 1.9 meters and 
Situated 12.5 meters from the burst point was determined. 

The calculated direct-induced ground shock arrives at the 
trench 400 usec after detonation, with a peak pressure of 

90 kilobars that decays to half value in 200 usec. This 
highly impulsive load initiates the collapse of the trench 
‘wall at the spall velocity (approximately 0.2 cm/usec). Due 
to the extremely transient nature of the direct-induced 
loading, little wall displacement actually takes place during 
the passage of the waveform. Since the direct~induced ground 
shock also decays very rapidly with slant range, only about 
12 meters of the trench are subjected to the very high initial 
closure velocities. However, airblast-dominated closure 
continues to pinch off the trench at much lower velocities 
out to several hundred meters. In defining the HE simulation 
experiment, only the simulation of the direct-induced ground- 
shock was considered. 


The simulation experiment consists of a rectangular 
slab of HE positioned over a buried trench. Calculations 
indicate that this source is capable of inducing a 90 kilobar 
shock at the trench wall. Moreover, the correct pressure-time 
history can be achieved by adjusting the thickness of the HE 
slab. The simulation experiments are conducted with 1/16th- 
scale trenches. The HE slab is 120-cm square by 6.93-cm thick 
and the scaled trench is a concrete pipe 26.5-cm ID with a 
2-cm wall. The primary instrumentation is high-speed photo- 
graphy to record details of the collapse process. Pressure 
in the soil and inside the rapid closure region is measured. 
Pressure measurements are also made to assess the magnitude 
of induced flow in the pipe outside the rapid closure region. 
However, the emphasis of the simulation experiment is on energy 
coupling and not on pipe flow. Results from a pre-shot 
simulation experiment are presented. 
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The 30-day study. group of MX trench testing requirements 


concluded that strong, direct-induced ground shock resulting 

from a near-miss nuclear burst is an energy source for trench 
collapse that is highly uncertain and has a significant impact 

on the environment definition. Present uncertainties in the 
modeling of direct-induced ground shock collapse are manifest 

as order of magnitude uncertainties in the plug pressure environ- 
ment. The study group further concluded that the data requirement 
for strong ground shock induced collapse was not being addressed 
in the then current MX experimental program, but that this re- 
quirement could be met by static HE tests. This paper will pre- 
sent the definition of the nuclear environment resulting from a 

l MT surface burst and the experiment designed tc simulate the 
effect of that environment on an MX trench. 


Systems, Science and Software under separate contract 
had performed a detailed calculation of a l MT surface burst 
dubbed source 3/5. This calculation used STREAK which is a 
multi-material Eulerian radiation-hydrodynamic finite difference 
code. The calculation included a detailed description of the 
source, the nuclear fireball, airblast coupling, and detailed 
zoning in the shock region - both direct-induced and airblast- 
induced. This code was used quite successfully in the previous 
experiments on the underground tests Ming Blade and Husky Pup, 
which were designed to validate the predictive codes used in 
energy coupling calculations. Source 3/5 was carried to a final 


time of approximately 800 microseconds at which point the airblast 


1, Review of MX Trench Testing Requirements, A Report of the 
Findings of the Joint DNA, SAMSO, and AFWL 30-Day Study Group, 
November, 1977, 
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is at a radius of approximately 85 meters and has a shock pressure 
of about 1.3 Giga Pascals (13 kbars). Figure 1 illustrates 
airblast values for a 1 MT surface burst. The curve marked 

"Max Surface Airblast" is the peak shock value of the blast wave 
at that time. The curve labeled "Min Surface Airblast" is the 
relatively constant (in space) value the pressure drops to be- 
hind the shock. Note that at times earlier than .1 msec these 
two curves come together. The well-formed blast wave does not 
appear until about this time following a 1 MT surface burst. 

The only data available lie at pressures below 100 MPa (1 kbar). 


The calculation was monitored at points corresponding to 
the top of the MX trench (approximately 1.5 meters below the 
original ground surface) at various radii. Two of these curves 
are shown as dashed lines in Figure 1. The close-in values 
(4 5 m) are not of interest to this problem since the question 
of direct coupling of the nuclear device to a trench was addressed 
in the underground test Hybla Gold. The 25 m curve on the down- 
side is mostly an extrapolation from the end of the calculation. 
This range reaches pressures of approximately 1 GPa (10 kbars) 
which is lower than values of interest. The 12.5 m range shows 
a peak pressure of approximately 9 GPa (90 kbars) which is above 
levels usually available through airblast coupling using HE. 

Yet this value is in the interesting range for trying to under- 
stand how pipe crush will couple energy down the length of the 


trench. 


Turning to Figure 2, in (a) is plotted the 100 GPa (1 mbar) 
and 10 GPa (100 kbar) isobars. This stress level does not occur 
at the same time at all locations so this plot is an amalgama- 
tion of several times. These curves are displaced downward at 
radii less than 5-10 m because the calculation included a cy- 
lindrical room roughly approximating another trench configura- 
tion directly under the explosion point. In Figure 2 (b) 
plotted as a solid line is the pressure at point A in 2 (a) as 
a function of time. The calculated ground shock arrives at the 
trench 400 microseconds after detonation, with a peak pressure 
of 9 GPa. This wave decays to half-value in 200 microseconds. 
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The highly impulsive load represented by this pressure-time 


history initiates the collapse of the trench wall at a spall 
velocity of approximately 0.2 cm/usec. Due to the transient 
nature of the loading, most of the wall displacement actually 
takes place after the passage of the waveform. Since the direct- 
induced ground shock also decays very rapidly with slant range 
{cf: Figure 2 (a)}, only about 12 meters of the trench are sub- 
jected to the very high initial closure velocities. However, 
airblast dominated closure continues to pinch off the trench 

at much lower velocities out to several hundred meters. Admit- 
tedly, there are uncertainties in this calculation cf the pre- 
dicted stress field following a 1 MT surface burst, but the 
calculation does suggest the environment to be simulated if one 
wants to study the energy untimately coupled to the pipe by the 
ground shock-induced collapse. 


Figure 2 (b) also illustrates the type of scaled pressure- 
time history (dashed curve) obtainable from HE coupled directly 
into soil. The simulation experiment consists of a rectangular 
slab of HE positioned over a buried trench. Calculations in- 
dicate that this source is capable of inducing a 9 GPa shock 
at the trench wall. Moreover, the correct pressure-time history 
at this point can be achieved by adjusting the thickness of the 
HE slab. (The decay of pressure with depth will probably not 
be the same in the HE case as in the nuclear, but the nuclear 
calculation has not been carried to the point where the shock 
has totally traversed the trench location. From earlier data, 
it appears that the nuclear value might reach about 2.5 - 3.0 
GPa at the trench bottom and the HE case about 1.5 GPa.) The 
simulation experiments are conducted with 1/16€th - scale trenches. 
Table 1 lists some parameters of the HE. The slab is approx- 
imately 1.2 meter square by 6.93 cm thick (Comp C4 with a total 
mass of 166 kg. (This is about the maximum surface charge 
which can be fired at the 5? Green Farm test site.) The areal 
size of the slab is dictated by efforts to avoid any edge 
effects affecting the ground shock until the wave has passed 
below the bottom of the trench. 
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-Table ie 


1/16 scale experiment 


un slab thickness 6255S 


Depth of pipe burial 7.62 cm + 2 cm 
+ liner thickness 
HE mass (1.2 m x 1.2 m) 166 kg 
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The scaled trench is a concrete pipe 26.5 em ID with a 
2 cm wall, 10 meters long. (The scaled depth of burial is 4 
feet, rather than the 5 proposed for an MX trench, but since we 
are interested in energy coupling to the remainder of the pipe 
by the collapse, this is not a serious difference and it does 
allow one to reduce the amount of HE required.) Since concrete 
pipes with walls this thin are not available commercially, the 
pipe was cast by s? in sections about 72 cm long in specially 
built forms. The pipe is mounted in a trench which is back- 
filled with Overton sand, the top of the pipe being 7.62 cm 
below the top surface of the sand. (See Figure 3 for the experi- 
mental set-up.) Overton sand was chosen as the soil material 
rather than alluvium because of the necessity to pack the material 
around the pipe with a relatively uniform density. In back- 
filling the trench, the sand is "rained" into the trench to a 
density of approximately 1.8 gm/cm?, with samples being measured to 
confirm the consistency. The HE charge is positioned on the 
surface of the sand directly above the center of the pipe as 
shown in Figure 3. (Note that the horizontal and vertical 
scales are different in this figure.) The HE charge is deton- 
ated simultaneously (to within + 2 usec) over its entire sur- 


face with two flyinc—plate-type detonators. 


The main diagnostic to measure the rate and shape of the 
pipe collapse is high-speed photography, using a framing camera 
running at about microseconds per frame. The camera looks 
directly down the axis of the pipe via the turning mirrors 
indicated in Figure 3. The pipe is back-lighted from the oppo- 
site end with a high-explosively~-driven argon "candle", the 
length of which is great enough that the shock wave does not 
arrive at the end of the candle until after the collapse of the 
pipe is complete. What is actually photographed is a shadow- 
graph of the collapsing pipe directly under the HE charge. In 
this way, both the rate of collapse and the cross-sectional 
shape of the collapsing pipe are determined. Using an opal 
glass window at the end of the argon candle was found to produce 
much more uniform backlighting than when a transparent window 


was used, 
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As indicated in Figure 3, various gauges are used to 
measure pressure and time of arrival both in the pipe and on 
the ground surface above the pipe. The three bar gauges shown 
are to measure the air pressure directly under the collapsing 
portion of the pipe. In addition to these gauges, two bar 
gauges are also located in the center plane perpendicular to 
the pipe axis, but at an angle of 20 degrees to the vertical. 
The purpose of these gauges is to measure the air pressure in 
the lobes which may form on either side of the axis during 
pipe collapse. The manganin flat packs on the bottom of the 
pipe under the HE charge are to measure the pressure produced 
by the impact on the bottom off the pipe of the collapsing upper 
portion of the pipe wall. In addition to these gauges, manganin 
gauges are also distributed in the sand at various distances 
below the HE charge. These gauges measure the ground shock 


incident on the pipe. 


TOA and pressure gauges located along the inside of the 
pipe but outside the region of the HE charge measure the time- 
of-arrival and magnitude of the pressure pulse driven down the 
Pipe by the collapse. Additional gauges on the ground surface 
above the pipe will measure the air blast incident on the ground 
outside the HE charge. Bundles of light pipes are mounted at 
various positions along the top inside of the pipe. These light 
pipes lead to flash tubes which are located beyond the end of 
the concrete pipe. The ends of the light pipes in the concrete 
pipe are pointed towards the turning mirrors and are photographed 
by a remotely-located HYCAM camera. The purpose of this arrange- 
ment is to measure the rate of collapse of the pipe outside the 


rapid closure region under the HE charge. 


In order to prove the concept and get ideas for timing 
settings on scopes in the main experiment, a set of two pre- 
shot tests were performed. These tests were mainly to check 
out the optics and argon candle performance. Commercial 6 inch 
concrete pipe and scaled HE charges were used. The first 
test indicated that careful alignment of all the optics was 


absolutely essential for obtaining pictures of high enough 


250 


quality to analyze. The second pre~shot test used the same 
scaling as the first, but the inside of the pipe was coated 
with a black, non-reflective material. Both opal glass anda 
negative density filter were used to get uniform illumination 
of the pipe cross-section and more than a day was spent in 


getting good alignment of the mirrors and the pipe. 


Figure 4 has four photographs relating to this 6 inch 
pre-shot test. Figure 4 (a) shows the pipe and turning mirror 
frame. The argon candle was later installed at the near end 
of the pipe. The bunker is beyond the truck visible in the 
background. Figure 4 (b) shows the sand bed over and around 
the pipe. A plastic tarp has been erected to keep direct sun- 
light off the HE after it is emplaced. The HE charge (approx- 
imately 3.5 cm thick in this case) and the flying plate detonator 
are visible in Figure 4 (c). A post-shot view of the test bed 
is shown in Figure 4 (d), with a technician standing in the 
main charge crater for scale. The large piece of plywood in 
the right foreground is part of the turning mirror framework 


seen in Figures 4 (a) and 4 (b). 


A STREAK prediction calculation for the full scale ex- 
periment was performed to determine if there were any potential 
problems with the design. This was done in X-Y plane geometry, 
taking a cross-section of the pipe,soil, and charge. We were 
interested in pressures in the air as the pipe squeezed off and 
stagnation pressures at the bottom of the pipe. This calculation 
did not include either radiation or strength, since neither 
was significant for this regime. The code uses tracer particles 
to delineate material boundaries, and this is what will be shown 
in figures to come. Both the inside and outside boundaries of 
the pipe will be presented in the tracer plots. They will be 
compared with the photographs taken during the second pre-shot 
test described above. The calculated times will be scaled to 
the times for this particular event. Comp B JWL equation of 
state coefficients were used in the calculation. 
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Figure 5 is the first of a series of plots comparing the 
event and the code. On the left [Figure 5 (a)] is a series of 
photos taken down the pipe while the collapse is in progress. 
There are 7.63 usec between frames and the collapse started at 
about Frame 9, which is not shown. To the right [Figure 5 (b)] 
is an artist's rendition (dashed line) of the pipe shape at the 
central frame of the three presented. The heading indicates 
it is for the 6 inch experiment, it is Frame 11, and the time 
is 84 usec from detonation of the main HE charge. The quality 
of reproduction of the photos is not good enough to show 
the sharp boundary which is apparent between the air and the 
Pipe wall. Figure 5 (b) also has the calculated configuration 
over-laid on the photo rendition. The shape is remarkably 
similar and the location is in reasonable agreement. (The cal- 
culation used an equation of state of alluvium at density approx- 
imately 2.1 rather than sand at a density of approximately 1.6- 


1.7 as in this experiment.) 


Figures 6 (a) and 6 (b) make the same comparison at 107 
usec, Figure 7 (a) and 7 (b) at 130 usec, Figures 8 (a) and 8 
(b) at 152 usec, Figures 9 (a) and 9 (b) at 176 usec, and Figures 
10 (a) and 10 (b) at 199 usec. In all of the photos there is 
no indication of a break-up of the pipe-wall/air interface. 
The calculation is in qualitative agreement with the observed 
behavior, although it is protruding more on axis. This is 
thought to be due simply to coarse zoning problems in the regions 
of the developing lobes. The lobe effect was not expected when 
the calculation was performed and thus was not adequately 
accounted for in the problem set-up. The air in these lobes 
reaches a calculated pressure of only about 100 MPa (1 kbar), 
which may not be effective in driving a strong shock down the 
length of the tunnel. The gauges on the main experiment will 
give data relating to the efficiency of the air compression and 
the validity of the numerical calculation. The experiment is 
currently set for detonation sometime during the week of December 
A. ASTSs 
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B) 


(A) 


Experimental/Calculational Comparison 


Figure 10. 


In summary, then, it is possible to simulate with HE the 
very close-in calculated nuclear blast-induced environments. 
To a certain extent, it is immaterial in this case how accur- 
ate (at least to within a factor of 2) the numerical prediction 
is. Given that the calculation is reasonably correct, the sim- 
ulated environment can be produced to allow a test of the 
effective energy coupling to a pipe as a result of collapse 
initiated by the strong, direct-induced ground shock. The 
photographic technique for measuring pipe collapse has been 
refined and confirmed. Finally, the evidence from the 6 
inch pre-shot test is that the high pressure, impulsive loading 
of the pipe does not cause a break-up at the pipe-wall/air 
interface as the collapse continues. Rather, a distinct separa- 
tion is maintained and lobes of trapped air form to either side 
of the vertical axis of collapse. The effect of this geometry 
on energy coupling down the pipe will be investigated in detail 


in the impending experiment. 
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COMPUTATIONAL STUDIES OF A VOITENKO COMPRESSOR™ 


P. S. Brown and M. L. Lohmann 
Lawrence Livermore Laboratory 
Livermore, California 94550 


ABSTRACT 


Two successful shock tube experiments’! were recently performed 
at the Lawrence Livermore Laboratory. The experiments simulated the 
high enthalpy flow conditions that might arise in the nuclear blast 
encounter of a tunnel-based missile. The experiments, using a modified 
Voitenko high explosive generator, produced megabar pressure conditions 
at the inlet of the 2 cm-diameter air-filled shock tubes. 

Computer calculations have been done to model the 10°-fold in- 
crease in air pressure in the Voitenko generator. These calculations are 
necessary to describe the source for the experiments and to permit better 
interpretation of the experimental data. A detailed understanding of 
the phenomenology of the compressor is also essential for optimization 
of future experiments, especially those involving scaling up to larger 
size shock tubes. 

An arbitrary Lagrange-Eulerian computer code was selected to calcu- 
late the late-time high explosive driven motion of a steel plate into the 
bell shaped compressor section. We found that ordinary Lagrangian and 
Eulerian codes had difficulty in properly treating the steel-air inter- 
faces, whereas an arbitrary Lagrange-Eulerian code was capable of modeling 
the complex flow of air past the steel interfaces in the compressor region. 
Using an improved equation of state for air, excellent agreement has been 
obtained in the source timing for the shock tube experiments, and we now 
have a much more detailed understanding of the Voitenko generator operation. 


*Work performed under the auspices of the U.S. Department of Energy under 
contract No. W-7405-ENG-48; and DNA Subtask J24AAXIX955. 
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INTRODUCTION 

In May and June of 1978 two shock tube experiments |!) 
at Lawrence Livermore Laboratory in support of the MX trench program. The 
purpose of the experiments was to simulate the high enthalpy flow condi- 
tions that might be encountered in the nuclear attack of a tunnel-based MX 
missile. The air-filled shock tubes, one made of steel and the other of 
grout, each had a 2-cm inside diameter and were driven by a high explosive 
Voi tenko 2,3) compressor. The aaa’ design was identical to that 
employed in a similar experiment 4 performed ten years earlier. 

It is desirable to understand how the Voitenko compressor works for 
purposes of designing a more efficient device, especially if the compressor 
is to be scaled to larger sizes. Accordingly, calculations were done to 
model its operation. An additional motivation for doing the calculations 
was to provide a better understanding of experimentally observed pressure 
pulse shapes and pressure oscillations. Calculations of the Voitenko 
generator are difficult because of the complex air flow in the late stages 
of the air compression, which necessitates an Eulerian hydrodynamic treat- 
ment. The difficulties are increased when dissimilar materials like air 
and steel are mixed in an Eulerian cell. The steel mass in the cell over- 
whelms the air, preventing the air from moving past the steel surface, 
and in effect, the air sticks to the steel. We used an arbitrary Lagrange- 
Eulerian code to overcome these difficulties. 


were performed 


Calculational Methods 

The Voitenko compressor is depicted in Figure 1. The air chamber 
pressure is 1] atmospheres, and the air in the exit shock tube is at 
1.0 atmospheres. The high explosive burn and the first 20 pseconds of 
the compressor operation were done with the Vector Hemp( 5-7) 2D Lagrangian 
code. The HEMP configuration at zero time is shown in Figure 2. The 
high explosive was initiated as a planar source along the left boundary 
of the problem, and the abscissa is an axis of symmetry. The HEMP confi- 
guration at 20 yseconds is shown in Figure 3. The initial motion of the 
steel plate into the chamber is planar except where the driver plate meets 
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the steel housing. This lack of planarity is not considered serious since 
this part of the plate soon impacts the hemispherical portion of the 
chamber. Only that portion of the HEMP problem between an axial position 
of 12 cm and 24 cm, and a radial position of 0 cm to 12.06 cm was linked 
to the Eulerian code. 

The calculation was continued using the arbitrary Lagrange-Eulerian 
code in which problems are described by a quadrilateral Lagrangian grid. 
The code basically does a Lagrangian treatment using the TENSOR 8,9 differ- 
ence equations. Individual nodes become Eulerian as needed, when the grid 
begins to distort or when relative zone thicknesses become too thin. 

At that time, the code does a second order accurate multifluid Eulerian 
treatment. The present calculations employed polynomial equations of 
state for the various materials. An improved equation of state!) was 
used for the air and was modeled in the code by a 28 term ratio of poly- 
nomials. The code has recently acquired a strength of materials treatment, 
an option that was not available at the start of these calculations. 


Results 


The configuration and zoning at 20 usecs are shown in Figure 4. 
The abscissa is an axis of symmetry and the upper boundary is a rigid wall. 
Figure 5 shows the configuration and velocity field at 23.2 psecs. The 
beginning of a radial air jet is visible at a radius of 9 cm; and 
Figure 6 shows an expanded view of this jet. We had originally tried an 
ordinary Eulerian code to do this calculation and were unable to observe 
a radial jet, even when extremely fine zoning was used. 


Figure 7 shows the configuration and velocity field at 27.25 usecs. 
For economy of calculation, part of the problem from a radial position 
of about 7.2 cm up to 12.06 cm has been discarded. Figures 8-13 show the 
further motion of the plate into the chamber and the resulting air flow 
field. Note the formation of a thin bladder of air between the steel 
plate and the compressor wall. The plate slows down at 29 usecs, and by 
29.69 usecs the bladder of air is expanding. 
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Figure 14 gives a closeup view of the air pouring into the central 
chamber region just prior to rupture of the diaphragm that separates the 
11 atm air in the chamber from the 1 atm air in the shock tube. We 
estimate that the diaphragm breaks at about 29.75 usecs which is in excellent 
agreement with the breakout times of the current experiments of 28.6 usecs 
and 29.1 usecs, and a breakout time of 30.1] usecs as previously repopted’ 
for one of the earlier experiments. The breakout time for the earlier 
experiment has been reduced by 0.6 usec to account for the shock transit 
time across a steel pad located between the plane wave lens and the high 
explosive. The breakout time is probably sensitive to the air chamber 
pressure which was somewhat different for the three experiments. 

Figures 15, 16 and 17 depict the initial flow of the air down the 
shock tube. There is a significant amount of distortion in the steel 
plate and the steel housing near the throat of the shock tube at late 
times and this appears to be restricting the flow from the bladder into 
the shock tube. We believe that incorporation of strength of materials 
in the calculation might lessen this distortion and lead to enhanced flow 
of air into the shock tube. At times from about 30.0 to 30.7 useconds, 
the downward radial air velocity in the bladder is about 0.7 cm/usec at 
a radial position of 2 cm, the velocity through the constricted region 
between the plate and housing is about 1 cm/usec, and the velocity down 
the shock tube is about 2 to 4 cm/usec. At 30.7 usecs we estimate that 
about 3 gms out of a total initial amount of 29 gms of air has entered 
the shock tube, and that about 2/3 of the total initial air is still 
available for flow into the tube (i.e., has not yet been trapped by the 
steel). 

Figures 18, 19 and 20 show the pressure profiles of the air shock 
developing at late times in the shock tube region. These calculations 
are now being continued in order to follow the air flow down to the first 
pressure gage at a position of 30.7 cm (10 cm from the diaphragm). We 
hope eventually to get a comparison of the calculated pressure pulse with | 
that measured at the 30.7 cm location. 
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Figure 21 gives the calculated energy delivered to the air as a 
function of time. The discontinuities in the curves at 27 yseconds result 
from discarding part of the problem above a radius of 7.2 cm. Only a 


fraction of the energy is lost in that process. 


Figure 22 compares the 


calculated plate and shock trajectories with those measured in the earlier 
experiments. The experimental points were taken from reference 4. 


Conclusions 


These calculations have given us a much better understanding of how 
the Voitenko compressor works. Its operation is considerably more compli- 
cated than we had originally conceived (we had not expected formation of 
an air bladder), and must be analyzed with relatively sophisticated 
computational methods. It appears that treatment of strength of materials 
in the steel sections may improve the computational results and allow the 
air to flow more readily into the shock tube; and it is planned to repeat 
the calculations including strength of materials. 

The computations will be continued until the shock wave passes 


the location of the first pressure gage. 
of the convergent wave interaction in the inlet region to the pipe in 
order to explain some of the experimentally observed pressure pulse 


shapes and oscillations. 
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ANFO FOR NUCLEAR WEAPONS EFFECTS SIMULATION 


by 
J. Petes 


Naval Surface Weapons Center 


LII-a 


ANFO for Nuclear Weapons Effects Simulation 


The history of ANFO (ammonium nitrate/fuel oil) as an explosive 
for use in nuclear weapons blast and shock effects simulation is a 
short one, hardly twelve years long (reference 1). In the mid-60's, 
the DNA (Defense Nuclear Agency) was seeking alternative explosion 
sources for use on its large scale military hardware and structures 
response operations. TNT (Fig. 1) which had been in use since the 
nuclear test moratorium had evidenced a number of significant draw- 
backs. For one, TNT was becoming shorter and shorter in supply 
because of the military demands brought on by the Vietnam conflict. 
For another, the cost of preparing large charges such as required 
for many of the test operations was becoming greater and greater. 
It was estimated that this total cost for material, preparation, and 
charge emplacement was one dollar a pound, or $1,000,000 for a 
500-ton charge. And on a technical basis, a serious drawback of the 
block built TNT charges was the unpredictable and severe blast anom- 
alies (reference 2) that were present upon detonation of the charge 
(Fig. 2). These anomalies, with jetting and ahead running blast 
waves, extended out to 1,000-ft from GZ and covered as much as 25-30% 
of the test area. They usually imposed much larger than expected 
blast forces on the test structures and made analysis of results 
difficult in some cases and, in others, invalidated expensive and 
important response studies. 


Aware of these problems, in 1966, the Naval Surface Weapons 
Center (formerly Naval Ordnance Laboratory) suggested the use of 
ANFO for DNA operations. Enough information about ANFO was 
available, primarily from the mining industry, to indicate the 
merits of this explosive. It was readily available; dozens of 
plants throughout the country made the material in huge quantities- 
approximately 500 million tons per year. It was cheap; in 1966 the 
cost of ANFO in 50-lb bags or in tank trucks delivered to a test 
site was about five cents a pound, (today, it is about ten to fifteen 
cents a pound). It is safe, relatively insensitive. AN and FO 
separately are not classed as explosives; mixed as ANFO, it is less 
sensitive to all the usual tests such as rifle bullet and impact, 
than TNT. And as deduced from its large use. by the mining industry, 
it is effective, effective in breaking and moving rock and earth. 
But in 1966 the questions were: Is ANFO effective as an airblast 
generator? Would it produce anomalies? Could it be made into 
500-ton charges? Would the output be predictable and reproducible? 
These and other questions not faced by the mining industry had to 
be answered before DNA and the military testing community would 
accept ANFO as a TNT substitute. 


Starting in 1966 with bootleg 8-lb ANFO charges and working 
up through 1976 with a large multi-agency effort, pre-DICE THROW II, 
a 120-ton ANFO charge, was fired (Fig. 3); most of the questions 
posed were satisfactorily answered. So, in October 1976, DICE THROW 
was conducted ~ the first full scale DNA sponsored operation using 
ANFO (Fig. 4). The charge contained 600-tons of ANFO and the results 
were good - predictable blast and ground shock, and anomalies of less 
extent and severity than on TNT shots (reference 3). This year in 
August, ANFO again was the explosive of choice when on operation 
MISERS BLUFF II six 120-ton charges were fired simultaneously. And 
on the schedule for 1980, there is planned a 600-ton shot for the 


289 


MISTY CASTLE series. ANFO has come of age. 


What about the questions asked? And perhaps more important, 
what are the answers, what are the characteristics of ANFO partic- 
ularly as it is used in large, multi-ton masses? As indicated 
previously, much of the early information on ANFO was developed by 
and for the mining industry for their particular needs. In most 
mining and earth moving operations, boreholes are drilled into the 
native rock and earth and the ANFO is poured into the hole. The 
hole diameters are relatively small, in the neighborhood of 6 to 
10 inches and because of the insensitivity of the ANFO, dynamite 
boosters are used at 10-ft intervals along the length of the bore- 
holes. Hence, much of the information available in 1966 was for 
small, highly confined charges. For nuclear weapon simulation blast 
effects tests, large quantities of explosive are required - in the 
hundreds of tons range - and confinement is to be avoided lest the 
confining case or structure produce unwanted and deleterious frag- 
ments. 


The mining industry had shown that unconfined charges had low 
output as indicated by detonation velocity (Fig. 5). We found 
Similar low output for small charges as measured both in detonation 
velocity and airblast. But as the charges became larger and larger, 
apparently self-confinement took place and both detonation velocity 
and airblast increased. In terms of TNT equivalence (Fig. 6), we 
found that as the unconfined charge size was increased from 8-lb 
sizes to 1,000-lbs and larger, the TNT equivalence for airblast 
increased from 0.47 to 0.82; beyond about 1,000-lbs, ‘the TNT 
equivalence remains as 0.82. Incidentally, this TNT equivalence 
is the reason why we now use 120-ton and 600-ton ANFO charges to 
reproduce the blast effects (beyond a few charge radii) of 100- 
and 500-ton TNT charges. (Because the detonation velocity and 
pressures of ANFO are less than those of TNT, we cannot get as high 
initial airblast pressures as generated by TNT; for most test 
operations this is of little consequence.) 


It is interesting to note that while TNT equivalence of ANFO 
has stabilized around 0.82 for large charges, the detonation velocity 
appears to be creeping up as the charges get larger (Fig. 7). 
Comparing charges with similar geometries, for hemispheres the 
velocity goes from 4,200 m/s to 4,600 m/s as the charge increases 
from 1,000-lbs to 100-tons. And for pre-DICE THROW II, a 120-ton 
capped cylinder, the measured average velocity was 4,790 m/s; for 
DICE THROW, a charge of similar geometry but weighing 620-tons, the 
detonation velocity was 5,020 m/s. 


Why this continuous increase in average velocity? Again, early 
data provide a plausible answer. Detonation velocity increases with 
increases in density (Fig. 8). Translating these data into terms 
appropriate to the large simulation charges, it seems reasonable to 
assume that as the charge gets larger, i.e., higher, the prilled ANFO 
in the lower layers of the charge become more tightly packed, more 
dense. Hence, it could be expected that higher velocities would be 
reached in the lower layers. On DICE THROW, where rather extensive 
detonation velocity measurements were made with rate sticks and by 
photo interpretation, the detonation velocity, indeed, was found to 
be a function of the height of the layers (Fig. 9). The lowest 
layer, which has to be the most densely packed because of all the ANFO 
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on top of it, had a detonation velocity of 5,600 m/s; the highest 
layer at which measurements were made - opposite booster #7 - had 
the lowest detonation velocity (reference 4). It would be nice to 
make in situ density measurements on future large ANFO shots to 
check this reasoning. 


Now, an obvious question can be asked: how does this chang- 
ing detonation velocity affect output - airblast, for instance? 
The evidence is that there is little practical effect. Needham 
(reference 5) calculated the overpressure and impulse as a function 
of range for a 600-ton ANFO sphere with two different densities, 
one with 0.8 grams/cc, and the other with 0.9 grams/cc. Por pressure 
(Fig. 10) and impulse (Fig. 11), you cannot tell the difference. And 
as shown in the previous figure, the fact that the TNT equivalence 
for airblast holds steady at 0.82 for charges larger than 1,000-lbs 
also indicates that output is not affected. It should be noted, 
however, that airblast is not a particularly sensitive measure of 
output for any charge ~ ANFO, TNT, or any other explosive. Ten 
percent differences in charge weight translate into only 2.1% 
differences in distances for a given pressure; this difference cannot 
be ascertained in the field where we are lucky to get + 10% answers. 
This is why, incidentally, we use the TNT equivalence of ANFO as 0.8 
instead of 0.82 for building large charges; it is easier to multiply 
by 0.8 and in terms of effects, we cannot tell the difference. The 
validity of the 0.8 equivalence and the insensitivity of airblast 
to density variations is summed up in Figure 12 showing the pre-shot 
predictions and the measured values on DICE THROW; excellent agree- 
ment and little scatter (reference 5). 


There is at least one other parameter of ANFO that merits 
attention - the fuel oil content. Commercial ANFO, the kind we use 
on test operations, is a stoichiometric mixture of ammonium nitrate 
and diesel fuel oil in a 94.5 to 5.5 weight ratio. This mixture, in 
which all the reactants are consumed without the need for atmospheric 
oxygen, gives maximum output as indicated by detonation velocity for 
a given density and degree of confinement (Fig. 13). Note that quite 
a spread in FO content can be tolerated ~ from about 3% to 10% - 
without severe reduction in detonation velocity. This is fortunate 
for simulation work because we have found that ANFO manufacturers 
apparently have some difficulty in maintaining a 5.5% fuel oil content. 
On DICE THROW the FO averaged about 6% for the whole charge but 
individual bags ranged between 4.9% to 7%. On MISERS BLUFF II, while 
the average FO content for all six charges was again 6%, one maverick 
bag was as low as 2% and another as high as 12.7%. 


Apparently these variations in fuel oil content do not affect 
the overall charge output; as shown before, pressure data fall on 
prediction curves. But these variations in a given charge do cause 
concern. On one hand, the FO content influences the sensitivity to 
initiation of the ANFO; on the other hand, the variations in fuel oil 
content introduce inhomogeneities in the charge. The lower the FO 
percentage, the higher the sensitivity (Fig. 14). The material even 
at 2% FO is still pretty insensitive but we would like to keep it 
at the 4-6% level where we get the best of everything.. 


The maverick bags with too high or too low FO percentages in- 
validate one of the possible and claimed advantages of ANFO charges - 
homogeneity of charge material throughout the charge. It was 
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inhomogeneities in TNT block built charges that were considered 
to be a source of anomalies (Ref. 2). It may be that the few 
anomalies still witnessed on ANFO shots, while less severe and 
less extensive than those seen on TNT shots, originate in these 
FO discrepancies. Better production control at the ANFO plants 
is called for. This may be nit-picking considering the overall 
advantages of ANFO over TNT, but if it is easily and cheaply 
attained, efforts to get uniform fuel oil content would be 
desirable. 


So much for some of the hydro- and thermodynamic properties 
of ANFO and how they impact on simulation tests where combined 
airblast and ground shock effects are of concern as on DICE THROW 
and MISERS BLUFF. ANFO because it comes in convenient 50-1b bags 
can be used in many configurations other than spheres, hemispheres, 
and cylinders. A few years back, ANFO was used for simulating 
eraters and direct-induced ground shock from a nuclear weapon surface 
burst (Ref. 6). In CHEST (Cratering High Explosive Simulation 
Technique) (Fig. 15), a preformed crater is lined with bags of 
ANFO and, using multipoint initiation, detonated. The radius of 
the preformed crater is selected such that the detonation pressure 
of the ANFO matches some iso-pressure contour of the nuclear burst. 
The ensuing crater formation and ground shock then is like that from 
the nuclear burst. This type of arrangement could be used for 
testing underground structures and facilities - missile silos, command 
and personnel shelters, utilities, underground stores and other 
underground structures. Another interesting technical point 
determined during the CHEST development is that the detonation 
pressure and velocity could be lowered to meet requirements by 
adding inert polystyrene beads to the ANFO to lower its density. 
(And, of course, detonation velocity and pressure can be increased 
by compacting the ANFO or by using smaller sized or groundup prills 
thus increasing the bulk density.) 


ANFO can be used at sea for testing the vulnerability of Navy 
ships and their equipment to nuclear weapon proportioned airblast. 
In one concept (Fig. 16), a large barge loaded with 500-tons, or 
so, of ANFO can be part of a towed or anchored array with test 
ships deployed around the charge at predetermined pressure levels. 
Many ships could participate in a test with realistic, underway 
conditions. 


ANFO could also be used by the Navy to simulate underwater 
shock waves. A floatable, compliant tank holding up to hundreds of 
tons of ANFO can be towed to sea, made to flood and sink to some 
predetermined depth where it would be detonated. Studies (Ref. 7) have 
shown that a system such as this would work down to a depth of 
5,000-ft, that the ANFO, which obviously would be compressed to 
high density, will reliably detonate. Thus, underwater nuclear 
weapon bursts could be simulated and with ships underway, realistic 
operational conditions could be achieved. In sea operations the 
known hygroscopicity of ANFO can be adequately countered by the use 
of tarpaulins and waterproof plastic bags. 
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In summary, ANFO is here for airblast and ground shock use 
on large target response tests; it can be used to simulate many 
hydrodynamic effects of military interest in most media*. The 
properties of the material are as well established as for many 
other explosives - and probably better established than for most 
so called "non-ideal" explosives. The material is still relatively 
cheaper compared to other explosives. It is easily available, easy 
to handle, and safe. As I once said, only the imagination of 
military and civilian engineers and scientists can limit the 
application of large masses of exploding ANFO to its problems; 
fortunately, this imagination is limitless. 


J. Petes 
Naval Surface Weapons Center 
Silver Spring, Maryland 


*Note that ANFO is not the cure-all for all nuclear weapons effects 
simulation techniques. Its limitations, e.g., large size required 
to get full output, must be observed. Thus, at present ANFO is 
eminently suited for large charge configurations; its limitations 
remain yet to be exploited. 
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FIGURE 3. PRE-DICE THROW I! — COMPLETED CHARGE 
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CONVENTIONAL HIGH EXPLOSIVES FOR NUCLEAR SIMULATION 


by 
Gordon F. Lederman, Jr., Capt., USAF 


Air Force Weapons Laboratory 


SO. 


Conventional High Explosives for Nuclear Simulation 


Abstract 


Conventional high explosives have been used as energy 
sources for nuclear simulators for quite some time. However, 
little thought has been given to developing a balanced methodology 
to address the technical and economic factors involved in optimizing 
the explosive used for a particular simulation technique. In this 
paper, an attempt is made to discuss the relevance of technical 
factors such as the ideality of the explosive, the detonation 
products, initiation characteristics, detonation velocity, critical 
or failure diameter, and whether or not the material is charac- 
terizable by standard tests. In addition, the often more complex 
economic factors of availability, production methods, field opera- 
tions, support facilities, and, of course, cost are also addressed. 
Many factors for discussion, such as safety and handling, include 
consideration of both technical and economic data. Specific field 
examples, as well as suggestions for future work, are discussed and 
several high explosive systems are detailed. This paper is meant, 
not as a complete technical evaluation, but rather as a managerial 
primer to bring forth topics for consideration. 
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As Mr. J. Petes of the NSWC points out in his paper, the 
non-ideal explosive, ANFO, has been and will continue to be 
widely used to simulate nuclear explosions. However, it is 
quite important to recognize the simulation limitations of 
ANFO and similar explosives. 


ANFO has been widely used as a surface-detonated explosive 
to simulate a nuclear event. In these experiments, no attempt 
has been made to significantly enhance the yield of the detonation 
by confinement or directing of the blast wave. Thus, this may be 
described as an unsophisticated simulation. That is, the scaled 
energy input of a nuclear device is input into the environment 
(ground and air) by producing that energy at ground zero with 
a non-nuclear (high explosive) energy source (Figure 1). For 
example, a tenth scale-simulation of a 1 MT nuclear device would 
require the detonation of 500 tons of TNT or approximately 600 
tons of ANFO. The explosive would be detonated at Ground Zero 
(GZ) and the surrounding environment would interact with the 
results of the detonation to produce the appropriate shock and. 
other interaction with the test articles. 


_ Obviously, such a simulation can easily escalate greatly in 
scope and cost, particularly if high yield weapons are being 
simulated at large scale. For example, a half-scale simulation 
of a 5 MT weapon would require approximately 375,000 tons of ANFO. 
At $200.00 per ton, the explosive cost alone would be $75 million. 
Such a test would require explosive quantities not readily avail- 
able, a massive construction project, and a test area encompassing 
several counties. This obviously absurd example is meant only to 
illustrate the necessity for nuclear simulators which focus energy 
and shock in order to reduce the amount of explosives and the cost 
and scope of the project. 


There are numerous ways of increasing the sophistication 
of a simulation and thus decreasing the scope and cost. The most 
obvious is by limiting the environment "exercised" to that per- 
tinent to the test structure (Figure 2). This limitation, as well 
as other methods of increasing the sophistication of the simulation, 
assume some degree of knowledge of the way in which nuclear detona- 
tions interact with the environment before the test structure is 
affected. The Dynamic Airblast Simulator (DABS) (Figure 3) and 
the High Explosive Simulation Technique (HEST) (Figure 4), dis- 
cussed elsewhere at this symposium, significantly limit the amount 
of environment exercised (References 1 and 2). Among other things, 
they move the energy source explosives closer to the target than 
the distance to GZ. The energy input is reduced further by con- 
fining ("tamping") the explosive in some way and by directing the 
shock. 
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Thus, to reduce the scope and cost of large scale high 
explosive nuclear simulations, we have developed techniques that 
require energy sources significantly smaller than multi-ton 
constructs of bagged ANFO. The larger of these simulators, 
because of the scale, yield being simulated, duration of simula- 
tion, and other factors, may still use many thousands of kilograms 
of explosives. Thus, budgetary constraints such as constructability 
of the explosive charge, may require the use of explosives other 
than the ideal explosives such as TNT. 


Non-ideal explosives, such as ANFO, suffer from several 
drawbacks that must be considered in their use. Ideal explosives 
have relatively small failure diameters (measured in several milli- 
meters) and properties such as detonation pressures (hundreds of 
kilobars) and velocities (5-9 mm/us) that are invariant with 
charge size above 10-20 mm diameter. These explosives are well 
characterized and their energy deposition is well characterized 
and reproducible. Total energy release from the explosive is 
within a few (less than 5) microseconds after passage of the 
detonation front. Non-ideal explosives often have failure dia- 
Meters in the range of several tens to several hundred milli- 
meters and their detonation pressures (usually in the tens of 
kilobars) and velocities (usually less than 6 mm/yis) vary signi- 
ficantly with charge diameter up to diameters of hundreds of 
millimeters. These explosives are not well characterized and, 
unlike ideal explosives, adequate equations of state (E.0O.S) do 
not exist for these explosives except under extremely limited 
circumstances. Total energy release from these explosives, if 
it takes place at all, may take up to several hundred microseconds 
after passage of the detonation front. Figure 5 summarizes these 
considerations. 


Thus, when using non-ideal explosives in simulation appli- 
cations, a number of considerations must be made that are per- 
tinent to the particular application. These factors are reviewed 
below. They may be grouped into the following categories: explo- 
sive properties, budgetary considerations, constructability, and 
safety. 


Explosive Properties 
The important explosive properties may be divided into two 
groups: detonation properties and physical properties. In the 


former category, I will discuss the properties important during 
energy release, such as detonation velocity, detonation pressure, 
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Figure 5 


COMPARISON OF IDEAL & NON-IDEAL EXPLOSIVES - 


IDEAL NON-IDEAL 

HOMOGENEOUS INHOMOGENEOUS 

SMALL FAILURE DIAMETER LARGE FAILURE DIAMETER 
MINIMAL DIAMETER EFFECTS DIAMETER EFFECTS 

SHORT RUN-UP LONG RUN-UP 


REACTION AT DETONATION FRONT LENGTHENED REACTION ZONE 
STEADY VELOCITY (7 - 9 mm/s) VARIABLE VELOCITY (2 - 6 mm/ps) 
HIGH ENERGY, PRESSURE LOWER ENERGY, PRESSURE 


KNOWN EQUATION OF STATE, PRE- EMPIRICAL EQUATION OF STATE 
DICTABLE 
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detonation chemical products, and explosive E.0.S. In the latter 
category, I will briefly touch on such topics as failure diameter, 
density, strength, and formulation. 


Detonation Properties 


Depending on the type of simulation, a precise knowledge 
of the detonation velocity of the explosive may be important. 
For example, in the HEST simulations, a cord explosive is normally 
used as the energy source. To simulate the nuclear shock front 
propagation velocity, the cord is angled to the direction of 
propagation. The component of the detonation velocity in the 
direction of shock front propagation becomes the shock front 
propagation velocity (Figure 6). In an application such as Dice 
Throw, a very accurate knowledge of the detonation velocity was 
neither attainable nor necessary. The detonation pressure refers 
to the pressure generated at the detonation front by the explosive. 
This will be reflected in the pressure-time history seen at a 
later time at a measuring point some distance from the explosive. 
In some instances, it may be desirable to input a given amount of 
energy as rapidly as possible--hence, an explosive with a very 
high detonation pressure may be desired. For purposes of sustain- 
ing a load, it may be desirable to avoid a high, initial pressure 
spike and instead, use an explosive with a lower detonation pressure 
and longer duration pulse. The chemical products from an explosive 
May be important for a particular application. In recent DABS 
work, the production of air-like detonation products was desired. 
This immediately eliminated from consideration explosives with 
very negative oxygen balances (high particulate carbon output) 
and explosives containing metals such as aluminum. All the above 
properties as well as others may be lumped together by con- 
sidering the E.0.S. of the explosive. Ideal explosives, as 
mentioned, have well-defined E.0.S.'s that allow a complete 
description of the explosive output given certain easily measure- 
able data such as density and chemical composition. Non-ideal 
explosives may have E.0.S.'s developed for them from experimental 
data that are applicable within a very narrow range of circum- 
stances. For example, an empirical E.0.S. was developed by 
Lawr ence Livermore Laboratory for the ANFO used in Dice Throw. 
Such £.0.S.'s, however, must be used very cautiously because of 
their limitations. The E.0.S. for a non-ideal explosive may even 
vary on a macroscopic scale within a single batch of such an 
explosive, because of compositional variations in the batch. 


Physical Properties 
The most important property in this category is the failure 
diameter and the related topic of the variation of other 


properties with charge size. The failure diameter may be defined, 
for cylindrical geometries, as that minimum diameter of explosive 
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which will sustain a high order detonation. The surface-to-volume 
ratio and the surface energy losses of an explosive charge deter- 
mine the failure diameter. This property can become important 

in two ways. If a charge is small enough in a particular dimen- 
sion, deflagration or low-order detonation may occur causing the 
explosive to release only a small fraction of its energy within a 
useful time period. Further, if a non-ideal charge is within 2 

or 3 multiples of its failure diameter, its detonation pro- 
perties may not be predictable. Thus, its usefulness for pre- 
dictive models may be limited. ANFO, for example, has an uncon- 
fined failure diameter of 100 mm. Further, its detonation proper- 
ties do not reach equilibrium values even at 300 mm (References 

3, 4, and 5). Thus, while it has reliable detonation properties 
in multi-ton piles, its properties may vary with thickness in many 
applications. The density of an explosive is an important property 
for two reasons. First, the density determines many secondary 
explosive properties, such as detonation velocity. Second, and 
particularly important if there are volume constraints, a denser 
material has a higher energy density than the same material ina 
less dense cast. The physical strength of the material relates 
primarily to field construction and will be deferred until the 
section on Constructability. 


Budgetary Considerations 


The first and foremost budgetary consideration is the avail- 
ability of the explosive. It does little good to determine that 
a particular explosive is the best for your simulation purposes if 
the only way it can be had in quantity is by building your own 
explosives' plant. A few examples will suffice. Amatol 80/20 
is, in terms of detonation velocity and pressure, predictability, 
physical strength, and detonation products, a very attractive 
explosive for use in a DABS facility. Because of its viscosity as 
a melt, however, no manufacturer is remotely interested in producing 
even small, experimental batches of it--hence, it is not readily 
available. ANFO, at the other extreme, is the most widely used 
explosive in the world. It is available from dozens of different 
manufacturers at prices from 10% down to 1% that of ideal explo- 
Sives. Related to availability is, of course, cost of the explo- 
Sive. This consideration is rather obvious and requires no dis- 
cussion. However, a hidden explosive cost should be mentioned: the 
manufacturer/customer relationship. If this relationship is 
friendly, it can mean a willingness on the part of the manufacturer 
to make slight, but important technical improvements without add-on 
costs or excessive contractual red tape. A related and often 
unconsidered cost of the explosive is the cost of the support 
facilities. These run all the way from cost of transport and 
cost of storage to the cost of installation and maintenance. 
Canadian experiments using detonable gas balloons have required 
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relatively high costs for compressed gas storage tanks, pumping 
facilities, and maintenance of secure and safe operations (Ref. 6). 
Liquid or slurry explosives might require on-site mixing, pumping 
facilities, and appropriate molds. 


The bottom line in all discussions of budgetary constraints 
is the "bang per buck." That is, for the energy delivered to the 
system desired, what is the total unit cost of the energy. This 
figure would naturally include all the factors discussed above. 
For simulations such as DICE THROW, ANFO is not only low in cost, 
but is also inexpensive to transport in bulk, easy to handle, 
store, and assemble in 50-1b sacks, and requires no special support 
facilities--it gives by far the most "bang per buck" for this 
application. 


Constructability 

Under the topic Constructability, I first include the manu- 
facturability of the explosive itself. As previously mentioned 
for Amatol 80/20 and also for some of the Amatex explosives, ; 
Manufacture may be difficult because of the viscosity of the melt. é 
This and other contributing factors, such as insufficient agitation 7 
and too rapid or too slow cooling of the melt may lead to poor * 
quality control of the explosive composition. Hence, an explosive ~ 
IMay perform other than expected simply because it literally is not 
the same explosive as it was thought to be. Additional quality 
control may be prohibitively expensive. 


A second subtopic under Constructability is the actual mounting 
of the explosive for use. If the explosive is gaseous, liquid, 
slurry, or a low-strength solid, it will require some sort of 
form. For a gas or liquid, this form would have to he a high- 
strength air- or liquid-tight container. In addition to required 
additional construction and loading equipment, the form material 
might very well become unwanted high velocity debris. A high- 
strength solid, such as the Iretol 30T22-C used in DABS facilities 
May provide their own support if they are suitably designed. 
Iretol is cast around dowels that are glued into holes drilled 
in a plywood mount. The mount is then attached directly to the 
facility wall (Figure 7). 


The initiation system must, of course, be an important con- 
Sideration in the construction of the simulator. The explosive 
Must be able to be initiated in a fairly simple, repeatable 
Manner. That is, the explosive cannot be insensitive to simple 
initiation, for example, by detonating cord or blasting cap. The 
initiation must be reliable and must perform in a repeatable, 
timely manner. For example, in a DABS facility, it is important 
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PHOTOGRAPH OF THE EXPLOSIVE DRIVER WALL 


FIGURE 7. 


OF A DABS FACILITY 
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to achieve near simultaneous initiation of the entire explosive 
wall. To do this, individual 300 mm x 300 mm square blocks of 

the explosive are all reliably initiated simultaneously with a 

detonating cord initiation train that is easily constructed 


‘beforehand and attached to the driver wall (Figures 8a.and 8b). 


Safety 


I have intentionally left safety for last, not to slight it, 
but rather to stress its importance. Utilizing a commercially 
available non-ideal explosive in itself helps to ensure that the 
explosive has been safely used in numerous, widespread applica- 
tions. This is not to say that safety can be taken for granted. 
Rather, especially during transport and storage, well-established 
safety procedures have been formulated and tested and should be 
adhered to. However, each new simulator use presents new 
problems in safety both in the handling of the explosives during 
mounting and in the mounted explosives themselves. For example, 
certain explosives, although drop insensitive, may be very highly 
skid sensitive. Skid sensitivity becomes extremely important 
where explosives are being moved horizontally over a rough surface. 
Sensitivity may increase dangerously at elevated temperatures, 
such as those found in Southern Arizona. Another safety issue that 
becomes important for larger quantities of explosives is the cri- 
tical mass--that mass above which an explosive's normal breakdown 
reactions may generate enough heat to initiate a sustained chain 
reaction. This mass may decrease at elevated temperatures, under 
confinement conditions, etc. 


Summary 


I have attempted to briefly discuss the uses of non-ideal 
explosives in nuclear weapons' simulation. We in the simulation 
community have turned to these explosives for numerous reasons. 
ANFO, which can be used in the larger, inefficient simulations, has 
several drawbacks for use in the smaller, more sophisticated 
simulations. In selecting an explosive for simulator use, factors 
of explosive properties, budgetary restraints, construction con- 
siderations, and safety must be addressed. These factors may 
seem obvious and trivial. If this were so, however, the simu- 
lation community would not have the history it has of difficulties 
with agreeing on the proper use and adequate prediction of the 
output of explosive drivers. Even at this writing, new explosives 
are being considered and tested as the energy sources in present 
simulators and for new simulation systems. 
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FIGURE 8b. PHOTOGRAPH OF THE EXPLOSIVE DRIVER 
IN A DABS FACILITY WITH THE DETONATING CORD INITIATION TRAIN 
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The list of considerations presented here is far from all- 
inclusive. To apply these and other criteria properly, a parti- 
cular simulator application must first be determined. Particu- 
larly important are the scale and yield and the sophistication 
of the simulation. All pertinent considerations must then be 
listed and weighted appropriately. Then, candidate non-ideal 
explosives may be compared. Figure 9 shows an evaluation chart 
used for DABS explosive evaluation. The final test, of course, 
is when the chosen explosive adequately performs the simulation 
task. The explosive energy sources for nuclear simulators have, 
for too long, been taken for granted. The explosive must be 
investigated early in the simulation design process, not after 
simulation difficulties begin to appear. 
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